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ABSTRACT 

Analysis  and  design  of  linear  feedback  control  systems  by  application 
of  the  basic  Mitrovic  method  is  well  established.   However   Li  tie  has  bf- - 
done  in  applying  the  basic  Mitrovic  method  to  nonlinear  systems.   The  objec 
tive  of  this  work  was  twofold.   First,  Mitrovic" s  method  as  applied  to 
systems  with  a  single  gain-variable  nonlinearity  predicts  limit  cycles  with 
results  comparable  to  existing  analysis  methods.   Secondly 9  air.  analog  :c  npv 
er  study  was  made  on  a  system  with  two  gain- variable  nonlinearities  and  8 
results  interpreted  by  Mitrovic" s  method. 

The  analog  computer  results  were  obtained  by  a  Donner  Scientific  Corp  i 
tion  analog  computer.,  Model  3100,  at  the  U.  S.  Naval  Postgraduate  Schc  -  . 
Monterey j,  California.   The  writer  wishes  to  express  his  appreciation  for  t  :  - 
assistance  and  encouragement  given  him  in  this  endeavor  by  Professor  Geos 
Thaler  of  the  U.  S.  Naval  Postgraduate  School. 
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PART  I 
PREDICTION  OF  LIMIT  CYCLES  BY  MITROVIC ' S  METHOD  FOR 
CONTROL  SYSTEMS  WITH  ONE  GAIN  VARIABLE  NONLINEARITY 
Introduction. 

Mitrovic's  method  is  easily  applied  to  nonlinear  control  systems  if 
the  nonlinearity  can  be  described  by  an  input-output  characteristic  in 
which  the  output  is  an  odd  single-valued  function  of  the  input.   Such  a 
nonlinearity  can  then  be  expressed  as  an  instantaneous  variable  gain 
with  zero  phase  shift.   This  property  allows  the  nonlinear  effect  to  be 
expressed  as  a  variable  coefficient  of  the  system  characteristic  equation. 
For  the  nonlinearity  in  the  error  channel,  such  a  system  can  only  have  limit 
cycles  if  the  linear  system  is  unstable  since  the  nonlinearily  contributes 
no  phase  shift. 

Mitrovic's  method  determines  absolute  stability  for  a  linear  system,, 
which  has  fixed  coefficients  in  the  characteristic  equation,  by  assuming 
two  of  the  coefficients  to  be  variable.   The  method  derives  equations 
for  two  coefficients  with  S   and  COfy  of  the  S-plane  as  parameters. 
With  the  JU)  axis  (3  =  0)  as  the  S-plane  mapping  contour,,  the  curve  of 
one  coefficient  versus  the  other  coefficient  is  plotted  using  different 
values  of  O^f)   .   If  this  curve  encircles  a  point  defined  as  the  working 
point  or  M-point  in  a  clockwise  direction  for  increasing  LOj^   ,  the 
system  is  stable.   The  working  point  is  determined  by  the  fixed  values 
of  the  two  coefficients  which  were  assumed  variable. 

For  a  nonlinearity,  as  described  above,  one  or  more  coefficients  of 
the  characteristic  equation  are  variable.   This  manifests  itself  on  the 
variable  coefficient  plane  as  a  locus  of  possible  working  point  locations. 


Thus  a  limit  cycle  for  the  nonlinear  system  exists  for  the  conditions 
at  the  intersection  of  the  working  point  locus  and  the  stability  (^  ■  0) 
curve.   A  typical  Mitrovic  plot  is  shown  in  Figure  I.   The  frequency  of 
the  limit  cycle  is  read  from  the  stability  curve  at  the  intersection. 
T,he  amplitude  is  determined  approximately  from  the  gain  of  the  non- 
linearity  obtained  from  the  working  point  locus. 
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Figure  I 
Coefficient  Plane  of  Mitrovic' s  Method 
The  amplitude  determined  in  this  manner  will  necessarily  be  less 
than  the  actual  amplitude.  '  The  gain  of  the  nonlinearity  is  defined  as 
the  ratio  of  the  instantaneous  output  to  the  instantaneous  input  and  the 
actual  amplitude  could  be  determined  accurately  if  during  a  cycle  the 


minimum/maximum  instantaneous  gain  could  be  determined.   However.,  the 
non-linear  gain  determined  from  the  Mitrovic  plot  is  a  "psuedo  -  average" 
gain  over  a  complete  cycle,  since  this  "psuedo  -  average"  gain  is  deter- 
mined from  an  "average"  M-point  position  over  a  cycle. 

The  shape  of  the  stability  curve  on  a  coefficient  plane  depends  on 
the  two  coefficients  chosen  and  the  order  of  the  characteristic  equation. 
If  the  coefficients  of  the  S°  and  S  are  chosen  variable.,  then  the  stability 
curve  for  a  third  order  system  is  a  straight  line  on  the  Bo   versus  q. 
plane,  for  a  fourth  order  system  the  stability  curve  is  a  parabola,  and 
so  on.   If  different  coefficients  are  chosen  variable,  then  the  stability 
curve  has  a  different  shape  on  the  corresponding  coefficient  plane. 

The  working  point  locus  depends  also  on  the  linear  system  and  speci- 
fically the  number  of  zeros  in  the  loop  transfer  function.   For  no  zeros 
in  the  loop  transfer  function  and  the  nonlinearity  in  the  forward  path, 
the  working  point  locus  is  as  straight  line  segment  parallel  to  the  3o 
axis  since  the  variable  nonlinear  gain  appears  only  in  the  constant  co- 
efficient of  the  characteristic  equation.   For  a  linear  system  with  one 
zero  in  the  loop  transfer  function,  the  working  point  locus  is  again  a 
straight  line  segment  with  a  specific  slope  on  the  3e>  versus  q,   plane. 
Other  systems  can  be  conceived  with  different  placements  of  the  nonlinear- 
ity but  the  basic  working  point  locus  is  always  a  straight  line  segment 
on  the  appropriate  coefficients  plane. 

The  use  of  Mitrovic' s  method  for  a  nonlinear  system  for  which  the 
linear  loop  transfer  function  has  two  or  more  zeros  becomes  more  compli- 
cated.  In  this  case  the  nonlinear  gain  appears  in  three  or  more  co- 
efficients of  the  characteristic  equation.   With  only  two  assumed  variable 


coef f icientSj  the  stability  curve  becomes  a  family  of  curves  with  the 
nonlinear  gain  a  parameter.   To  determine  the  limit  cycles  for  such  a 
system  the  intersection  of  the  straight  line  working  point  locus  and 
the  family  of  stability  curves  must  be  found. 

The  various  systems  chosen  to  analyze  in  the  chapters  of  Part  One 
are  completely  arbitrary.   The  author  had  no  specific  system  or  specific 
device  in  mind  when  the  block  diagram  was  drawn.   Consequent ly9  all  signals 
in  the  various  systems  are  referred  to  in  volts.   The  following  guidelines 
were  followed  in  order  to  simplify  the  analysis  without  disrupting  the 
arbitrary  nature  of  the  system : 

1.  The  poles  and  zeros  of  the  linear  loop  transfer  function 
were  in  general  chosen  one  octave  apart.   This  was  done  in  order  to 
sketch  the  actual  Bode  magnitude  curve  in  the  describing  function  anal- 
ysis from  the  asymptotic  plot  by  applying  the  usual  corrections  at  the 
corner  frequency  and  two  octaves  either  side  of  the  corner  frequency. 

2.  The  magnitudes  of  the  poless  zeros,,  and  gain  constants 
were  chosen  to  eliminate  magnitude  and  time  scaling  procedures  in  the 
analog  computer  simulations  of  the  system. 


CHAPTER  I 
SATURATION  IN  ERROR  CHANNEL 
I-A  Loop  Transfer  Function  with  No  Zeros. 
%        The  block  diagram  of  the  system  chosen  to  analyze  is  given  in  Figure 

I-A. 


Figure  I-A 

Saturation  in  Error  Channel 
Loop  Transfer  Function  with  No  Zeros 


I-A-l.   Mitrovic's  Method. 

Assigning  the  variable  gain  N  to  the  nonlinearity,  the  character- 
istic equation  becomes, 


s3-r3s2+  as  t  /aw  =o  . 


(I-A-l.l) 


Assuming  the  last  two  coefficients  variable  the  characteristic 
equation  is  rewritten, 


53-+  35%  6,S  +  8d  =0 


(I-A-1.2) 


From  Appendix  A,  Mitrovic's  equations  for  B  and  B  are, 

e0  =  -feuM^)  +  ojo30z(^)]  d-A-1.3) 


and 


B,=  3  00n$2^)  +  uV<j53(S)  .  (I-A-1.4) 

Substitution  of  the  values  of  the  0  functions  Appendix  A  for  ^  =  0 
gives  the  parametric  equations  of  the  stability  curve  as 

8<>  -  3  6LV,2-  (I-A-1.5) 

and 

6,=  (X>r?     •  (I-A-1.6) 

From  comparison  of  equations  I-A-l.l  and  I-A-1.2,  it  is  seen  that 
the  equations  describing  the  M- point  locus  are 

8$  =  \0  N  (I-A-1.7) 

and 

B,  ^  2.  ,  (I-A-1.8) 

Thus  the  set  of  four  parametric  equations  I-A-1.5  through  I-A-1.8 
describe  the  two  curves  whose  intersection  is  desired  for  limit  cycle 
determination. 

At  the  intersection,  Equation  I-A-1.5  equals  Equation  I-A-1.7  and 
Equation  I-A-1.6  equals  Equation  I-A-1.8.   Thus, 

LOf\       =2-  (I-A-1.9) 


and 


0  Nave  =  3  lo0  . 


(I-A-1.10) 


Solution  of  these  two  simultaneous  equations  yields, 

U^=   1,414    r<xd/s<sc  .      , 


(I-A-l.ll) 


and 


Nave  -  0  .  h  • 


(I-A-1.12) 


Since  the  nonlinear  gain  N  is  defined  as  the  ratio  of  the  output 
to  the  input, 


rr\&*  -  n 


<LV<2 


(I-A-1.13) 


Thus  the  frequency  of  the  limit  cycle  is  given  by  equation  I-A-l.ll 
and  the  approximate  amplitude  by  I-A-1.13  or 


6m*y  =3.33  Volts  . 


(I-A-1.14) 


The  above  solution. is  presented  graphically  in  Figure  I-A-l.l. 
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Figure   I-A-l.l 


B, 


Graphical  Solution  for  Saturation  in  Error  Channel,  Loop  Transfer 

Function  with  no  Zeros 


The  limit  cycle  is  a  stable  limit  cycle  as  determined  by  the  follow- 
ing reasoning.   If  the  input  signal ,  Q^  ,  is  large  enough  such  that  N 
is  small  enough  to  place  the  instantaneous  M-point  at  say  f/\  »   in  Figure 
I-A-l.l.   The  system  is  instantaneously  stable  and  the  oscillations  de- 
crease.  As  the  oscillations  decrease9  N  increases  moving  the  instantan- 
eous M-point  above  the  stability  curve.   Here  the  system  is  unstable  and 
the  oscillations  increase  thereby  decreasing  N  and  moving  the  M=point  back 
into  the  stable  region.   Eventually  a  dynamic  equilibrium  will  occur  in 
which  the  "pseudo  average"  M-point  lies  on  the  intersection  of  the  two 
curves.   If  the  input  signal  is  small  enoughs  N  will  be  large  enough 
such  that  the  instantaneous  M-point  lies  at  say  fVL  at  which  point  the 
system  is  unstable.   By  similar  reasoning  as  aboves  the  "pseudo  -  average" 
M-point  will  again  eventually  end  up  at  the  intersection. 


I-A-2.   Root  Locus  Method. 

The  root  locus  method  of  analyzing  a  system  with  such  a  non-linearity 
consists  of  construction  of  the  linear  root  locus  and  determining  the 
point  where  the  root  locus  crosses  the  J60  axis. 

The  system  characteristic  equation  is 

5*  +  3s2  t  2  s  +  ION  =  0  .  (i-a-2.1) 

Substitution  of  S=-i^->  an<3  simplifying 

-Jw3-  3  to*  +j2tu  +  I0M  =0  .  (i-a-2.2) 

By  requiring  that  both  the  real  and  imaginary  part  of  equation  I-A- 
2.2  go  to  zero  independently , 

-JtO3  +  J  2  60  =  0  (I-A-2. 3) 


and 


-OOJ   +  I0M  =  0  •  (I-A-2. 4) 


Dividing  equation  I-A-2. 3  by   J^O  the  frequency  is  obtained  as 

LO  -   1.414  rad/sec. 

By  substituting  (jq     into  equation  I-A-2. 4  the  "pseudo  -  average" 
gain  over  a  complete  cycle  is  obtained  as 

Mai/S  =0.6 

from  which  the  approximate  amplitude  of  the  limit  cycle  is  determined 
as 


Ill   1 1  I II  I     ll  I  I  .. 


A  sketch  of  the  root  locus  solution  Is  given  in  Figure  I-A-2.1. 
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Figure  I-A-2.1 

t 
Root  Locus  Solution  for  Saturation  in  Error  Channel 

Loop  Transfer  Function  with  No  Zeros 
The  root  locus  method  also  gives  a  stable  limit  cycle  by  the  follow- 
ing reasoning.   If  the  input  signal  is  large  enough,  N  will  be  small 
enough  such  that  the  instantaneous  root  location  is  at  f^     in  Figure 
I-A-2.1.   The  system  is  stable  and  the  oscillations  decrease.   As  the 
oscillations  decrease  N  increases  moving  the  instantaneous  root  into  the 
right  half  plane.   The  system  is  now  unstable,  the  oscillations  increase 
causing  N  to  decrease,  and  the  instantaneous  root  nioves  back  into  left 
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half  plane.   Eventually  a  dynamic  equilibrium  is  established  in  which  the 
"average"  root  location  is  at  the  imaginary  axis  cross-over.   If  the 
input  signal  is  small,  N  is  large  placing  the  instantaneous  root  in  the 
right  half  plane  at  ITq,.  .   By  a  similar  reasoning  the  "average"  root 
location  over  a  cycle  will  again  lie  at  the  imaginary  axis  cross-over. 
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T-A-3.   Describing  Function  Method. 

Application  of  describing  function  theory  to  the  given  system  entails 
plotting  the  linear  loop  transfer  function  and  the  negative  reciprocal  of 
the  describing  function  on  the  gain-phase  plane  (Nichols  plot).   Limit 
cycles  of  the  system  are  then  given  by  the  intersections  of  the  two  curves 
on  the  gain-phase  plane. 

The  describing  function  for  a  saturating  element  with  unity  gain  in 
the  linear  region  is 

GD  =  #■[$!*"' R  +  RVi-R*J  Zo°  (i-a-3.1) 


where 


R  =  4s**"  ^A-3-< 

t  max  • 

The  linear  loop  transfer  function  in  frequency  response  form  is 

Since  G  has  no  phase  shift  thenl/G  has  an  associated  phase  angle 
of  -180° j  and  the  1/G  curve  lies  entirely  on  the  -180°  axis  of  the  gain- 
phase  plane.   Consequently,  the  only  region  of  the  linear  loop  transfer 
function  which  is  of  interest  is  the  region  in  which  the  phase  shift  is 
-180°.   The  values  of  gain  in  db  and  phase  of  the  linear  loop  transfer 
function  in  the  region  of  -180°  which  were  used  to  plot  on  the  gain- 
phase  plane  are  given  in  Table  I-A-3.1.   These  values  were  obtained  from 
a  Bode  diagram  plot  of  the  linear  loop  transfer  function. 
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TABLE 

I-A-3.1 

CO 

IGMUb 

&QUJ   d&3 

1.0 

9.8 

-162 

1.1 

8.5 

-167.5 

1.2 

7.1 

-172 

1.3 

5.8 

-176 

1.4 

4.5 

-179.5 

1.5 

3.3 

-183 

1.6 

2.2 

-186 

The  values  of  1/G  in  db  were  calculated  for  various  values  of  R 
according  to  equation  I-A-3.1  and  appear  in  Table  I-A-3.2.  The  gain- 
phase  plane  of  the  system  is  shown  in  Figure  I-A-3.1. 


TABLE  I-A-3.2 


R 

So 

l/So 

\/SD     elfe 

0.0 

0 

OO 

©C 

0.1 

0.127 

7.88 

17.92 

0.2 

0.253 

3.95 

11.94 

0.3 

0.376 

2.66 

8.5 

0.4 

0.496 

2.02 

6.1 

0.5 

0.609 

1.642 

4.32 

0.6 

0.716 

1.396 

2.9 

0.49 

0.598 

1.673 

4.48 
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By  interpolation  of  Figure  I-A-3.1,  the  intersection  occurs  at 

(JO  =  1.41  rad/sec 
and 

R  =  0.495    0 
Thus  from  Equation  I-A-3.2  the  amplitude  of  the  limit  cycle  is 

6*ax=  10,1  volts   • 

Assuming  the  1/G  curve  to  be  the  locus  of  the  critical  point  of 
the  system,  the  describing  function  method  predicts  a  stable  limit  cycle 
by  the  following  reasoning.   If  the  input  signal,,   @>J.   is  large  R  will 
be  small  and  the  instantaneous  critical  point  lies  above  the  GC^)  curve 
say  at  point  R  in  I-A-3.1.   The  system  is  momentarily  stable  due  to  a 
positive  phase  margin  and  the  oscillations  will  decrease.   As  the  oscil- 
lations decrease,  R  increases  and  the  critical  point  moves  below  the 
0  0^0)  curve.   The  system  is  now  unstable  due  to  a  negative  phase 
margin  and  the  oscillations  increase  moving  the  critical  point  above  the 
GO^i)  curve.   Eventually  a  dynamic  equilibrium  is  reached  in  which  the 
"average"  critical  point  over  a  cycle  lies  at  the  intersection  of  the  two 
curves.   If  the  input  signal  is  small,  R  is  large  and  the  critical  point 
is  momentarily  below  the  (90u))curve  at  $k   ,   The  system  is  unstable 
and  the  oscillations  increase.   Similar  reasoning  shows  again  that  the 
"average"  critical  point  over  a  cycle  will  eventually  lie  at  the  inter- 
section of  the  two  curves. 
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I-A-4.   Analog  Computer  Simulation. 

The  simulation  of  the  given  nonlinear  system  is  given  in  Figure 
I-A-4.1.  The  coefficient  pot  settings  with  associated  resistances  and 
capacitances  for  real  time  and  real  magnitude  scaling  are  given  in  Table 
I-A-4.1. 

TABLE  I-A-4.1 


Pot  Setting 


al 

100 

a2 

-  1.0 

a3 

=  1.0 

a4 

=  1.0 

a5 

-  1.0 

a6 

=  1.0 

a7 

=  1.0 

a8 

-  1.0 

a  = 

=  0.05 

Initial 

Set 

ting 

Associated  Elements 


For  Input  Step  of  E.  volts  and 


Rx  =  1  Meg, 


R, 


2=1  Meg 


R3  =  1  Meg 

R.  =  1  Meg 


R,.  =  1  Meg 


R,  =  0.1  Meg 
R,  =  0.5  Meg 


Rg  =  1  Meg 


For  a  5 

Voltage 


fl 


^fl 


m  1  Meg 


Meg 


Rf2  a  1  M^g 


Cc,  -  1  uf 

Cfl  -  1  uf 

cf2  =  I«f 

Cf2  -  1  uf 

C-„  =  1  uf 
f3 


To  accurately  set  the  back-biasing  potss  "a"s  a  sine  wave  e     1 :r 

was  connected  to  the  input  of  the  saturation  simulator.   The  two  pots 
were  then  adjusted  until  the  saturated  output  voltage  was  exactly  5  volts 
on  a  recording  device.   The  frequency  of  the  generator  was  set  at  0.01  cps 
and  the  static  characteristic  of  the  saturation  simulation  was  obt 
by  an  X-Y  recorder.   This  characteristic  appears  as  Figure  I-A-4.2. 
The  system  was  first  excited  with  a  2  volt  step  input  signal. 
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system  was  found  to  go  in  to  the  limit  cycle  shown  in  recorder  trace 
I-A-4.1.   The  system  was  then  excited  with  a  20  volt  step  input  and  was 
found  to  go  into  the  same  limit  cycle  (recorder  trace  X-A-^.2).   From 
the  recorder  traces  the  limit  cycle  is  a  stable  limit  cycle  with  an 
angular  frequency  of  1.395  rad/sec  and  10.0  volts  amplitude.   The  phase 
portrait  for  the  system  appears  in  Figure  I-A-4.3. 
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CHART  NO.  RA2921  30  BIIUSH      IXSTKU.\IE\TS 


Recorder  Trace  I-A-4.3 

Limit  Cycle  Measurement s  Error  Channel  Saturations,  Loop 
Transfer  Function  with  No  Zeros.  • 


«\ 


Scales:     £  "  1  Volt/Line    X  -  0.5  Volts/  Line 

Time  -  0.2  Second/Divisior 
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I-A-5.   Comparison  of  the  Four  Analytic  Methods  in  Limit  Cycle  Pre- 
diction for  the  System  of  Figure  I-A. 


Method 

Frequency 

Amplitude 

(rad/sec) 

(volts) 

Mitrovic*  s 

1.414 

8.33 

Root  Locus 

1.414 

8.33 

Describing  Function 

1.41 

10.1 

Analog  Simulation 

1.395 

■■■  -      ■■■  -   — ■ 

10.0 

24 


I-B.   Loop  Transfer  Function  with  One  Zero. 

The  block  diagram  of  the  system  chosen  to  analyze  is  shown  in 
Figure  I-B. 


s»-*K8 


$+8) 


Figure  I-B 
Saturation  in  Error  Channel  Loop 
Transfer  Function  with  One  Zero 
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I-B-l.   Mitrovic's  Method. 

If  N  is  the  instantaneous  gain  of  the  nonlinearity  the  character- 
istic equation  for  the  system  is 

S4+  l3S3  +  44s2+(3Z-f30ON)s  +  C>00N-0  .  (i-b-i.i) 

Assuming  the  last  two  coefficients  to  be  variables,  the  characteris- 
tic equation  becomes, 

S^+  l353  +  44s1  +  6,  s  ■+  Bo  =  0  (i-b-1.2) 

where 

Bj  =  32  +  300  N  (I-B-l. 3) 

and 

30  =  600  N    ,  (I-B-l. 4) 


From  Appendix  A,  Mitrovic's  equations  for  B  and  B-  are 


8o=-[^6O^(«)+i3^302(^)f  Wh4£«0j  a 


-B-1.5) 


and 


6,  =  W^foCt)*  l3w^$/g)  +  w^tt)  . 


(I-B-l. 6) 


Substitution  of  the  values   of  the   0    functions    from  Appendix  A   for 
Q     =  0  gives, 

80  -44u>n2-UGn4  (I-B-l. 7) 

and 

8,   -    13  60n  ^         .  (I-B-l. 8) 

Thus   the   sets   of  parametric  equations   I-B-l.  3,    I-B-l. 4,    I-B-l.  7, 
and   I-B-l. 8   describe   the   stability  curve   and   the  working  point    locus   on 
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the  B  versus  B.  plane.   To  find  the  intersection  of  two  curves  solve 
o         1 

Equation  I-B-1.8  for  CO^     and  substitute  in  Equation  I-B-1.7  to  obtain 

A  A  B,2- 

Solving  Equation  I-B-1.4  for  N  and  substituting  in  Equation  X-B-1.3 
gives, 

6,  =  32+  £>.£"  B0  (i-B-i.io) 

At  the  intersect  ion* equations  I-B-1.9  and  I-B-1.10  must  be  satis- 
fied simultaneously.  Solving  Equation  I'E-1.  10  for  Bn  and  substitut- 
ing in  Equation  I-B-1.9  yields, 

g^-234  8,  "  10816  =  0  (I-B-l.ll) 

Solving  for  B..  gives 

B.  =  273.54,  and  B.  =  -39.54 
1  1 

The  negative  value  is  meaningless^  thus  substitution  of  the  posi- 
tive value  into  equations  I-B-1.8  and  I-E-1.3  gives  the  frequency  and  the 
"average"  gain  respectively,  as 

COn   =  4.587  rad/sec 
and 

Na,e=  °'8051 
The  approximate  amplitude  of  the  limit  cycle  is  given  by 

E 

C.  =        sat 

max     N 
or 

£     =  TTT^i      =  6.21  Volts 
^  max    0.8051 

The  graphical  representation  of  the  solution  is  shown  in  Figure 
I-B-l.l. 
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Figure  I-B-l.l 

Graphical  Solution  for  Saturation  in  Error  Channel, 

Loop  Transfer  Function  with  One  Zero 

Mitrovic's  method  predicts  a  stable  limit  cycle  by  following 
reasoning  to  that  given  in  Seeder.  I-A-l  of  this  chapter. 
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I-B-2.   Root  Locus  Method. 

The  system  characteristic  Equation  is 

s4  +  13s3  +  44s2  +  (32  +  300  N)s  +  600  N  =  0.  (I-B-2. 1) 

Substitution  of  s  =  J  LO  and  simplifyings  yields 

u/-Jl3u^-44toa  +  >)(32-/-300N)a)+  hOOH  -  O   *  (i-b-2.2) 


By  requiring  that  both  the   real   and   imaginary  parts   of  equation 
I-B-2. 1   go   to   zero   simultaneously;,    then 

uf-/\4U)2+   600  N  =  0  (I-B-2 .3) 

and  -tteai>32«)+"  J3006ON       =  0.  (I-B-2. 4) 

Dividing  Equation   I-B-2. 4  by    J(/J      and   solving   for  N  gives 

N  ■     300 W*  "    300  •  (1-1-2.5) 

Substituting   this    into  Equation   I-B-2. 3  gives   the    following  equa- 
tion  for    (JO 

A  2. 

CO i    -l$£0  -   64  =  0     «  (I-B-2. 6> 

2 
Solution  of  Equation  I-B-2. 6  for  UJ       gives 

6U2  =  21.042, 

from  which 

CO  =   4.587  rad/sec. 

2. 

Substitution  of  CU   into  Equation  I-B-2. 5  gives  the  "average""  non- 


linear gain  over  a  cycle  as 


N    =  0.8051 
ave 


From  this  the  approximate  amplitude  is  obtained  as 

£•   »   =  n  BL    =  6-21  Volts 
max    0.8051 

A  pictorial  representation  of  the  root  locus  solution  appears  in 

Figure  I-B-2. 1. 
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The  root  locus  method  also  predicts  a  stable  limit  cycle  by  follow- 
ing similar  reasoning  to  that  given  in  Section  I-A-l  of  this  Chapter. 


LINEAR    SYSTEM   CLOSED- 
LOOP   RCOT 


-<r 


■*■ 


•8 


LIMIT  CYCLE 
6J-4.f87  rod/sec. 

iY*„«=o.8af 


Figure  I-b-2.1 
Root  Locus  Solution  for  Saturation  in  Error  Channel 
Loop  Transfer  Function  with  One  Zero 
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I-B-3.   Describing  Function  Method. 

The  describing  function  is  the  same  as  that  given  in  Equation  I-A-3.1 
of  this  Chapter. 

The  linear  loop  transfer  function  in  frequency  response  form  is,, 


JuXJ(jU  +  lXO,25-JW+0(0,l2rjO>+/) 


(T-B-3.1) 


Table  I-B-3. 1  gives  the  magnitude  and  phase  of  Equation  I-B-3. 1 
which  were  used  to  plot  the  linear  loop  transfer  function  on  the  gain- 
phase  plane  of  Figure  I-B-3. 1.   These  values  were  obtained  from  a  Bode 
diagram  plot  of  the  linear  loop  transfer  function. 

TABLE  I-B-3. 1 


(JO 

|SOco)/Jb 

/60*>)  ^25    ' 

4.05 

4 

-174.5 

4.3 

3 

-177 

4.55 

2 

-180 

4.85 

1 

-183 

5.15 

0 

-190.5 

Table  I-B-3. 2  gives  the  values  of  1/G   in  the  area  of  interest 
These  values  were  obtained  from  Equation  I-A-3.1. 


TABLE  I-B-3. 2 


R 

GD 

I/Gd 

l/So   <**> 

0.6 

0.716 

1.396 

2.9 

0.68 

0.794 

1.26 

2.01 

0.69 

0.803 

1.246 

1.92 

0.7 

0.812 

1.232 

1.83 
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By  interpolation  of  Figure  I-B-3.1,  the  intersection  and  thus  the 

limit  cycle  occurs  at 

(jj    =  4.55  rad/sec. 
and 

R  =  0.682 
From  Equation  I-A-3.2,  the  amplitude  of  the  limit  cycle  is  then 


a 


=   7.33  Volts 


max    0.682 

The  describing  function  predicts  a  stable  limit  cycle  by  following 
similar  reasoning  to  that  given  in  Section  I-A-3  of  this  Chapter. 
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I-B-4.   Analog  Computer  Simulation. 

The  analog  simulation  of  the  given  nonlinear  system  is  sh 
Figure  I-B-4.1.   The  coefficient  pot  settings  with  associated  resistances 
and  capacitances  for  real  time  and  real  magnitude  scaling  are  given  in 
Table  I-B-4.1. 

TABLE  I-B-4.1 


Pot  Set tins 

Associated  Elements 

Ei 
ai  =  100 

For  Input.  Step  of  E,  Volts 

and  R  =  1  Meg,   R   ~   1  Meg 

a2  =  1.0 

R2  -  1  Meg,    Rfl  -  1  Meg 

a3  =  1.0 

R3  =  1  Meg,     Rf2  =  1  Meg 

a.  =  1.0 

4 

R,  =  0.5  Meg,   Cf,  =  1  uf 

a  =  0.4 

R_  =0.1  Meg,   Cr1  =  1  uf 

5        &9    fl 

a,  =  1.0 

6 

R,  =  0.1  Meg,   R   =  1  Meg 

a?  -  1.0 

R?  =  0.1  Meg,   R   =  1  Meg 

ag  =  0.3 

Rg  =0.1  Meg,   Cf2  =  1  uf 

ag  =  1.0 

R9  »   1  Meg,     Cf2  -  1  uf 

a1Q  .  1.0 

R   =0.1  Meg,  Cf3  -  1  uf 

ail  "  °'8 

R   =  0.1  Meg,  G   =  1  uf 

al2  =  1.0 

R.,  -  1  Meg,    C-.  =1  uf 
12       °J           r4 

a^  .  1.0 

R.  ,  =  1  Meg,    R„.  =  1  Meg 
13            f4 

a  =  0.05 

For  a  5  volt 

initial  Setting 

Saturation  Voltage 

To  accurately  set  the  back=bias  pots  "a"  a  sine  wave  generator  was 
connected  to  the  input  of  the  saturation  simulator.   The  two  pots  were 
then  adjusted  until  the  saturated  output  voltage  was  5  volts  on  a  re- 
cording device. 
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The  static  characteristic  for  the  saturation  simulator  appears  as 
Figure  I-A-4.2  in  Section  I-A-4  of  this  Chapter, 

The  system  was  first  excited  with  a  4  volt. input  signal.   The 
system  was  found  to  go  into  the  limit  cycle  shown  in  recorder  trace 
I-B-4.1.   The  system  was  then  excited  with  a  10  volt  step  input  signal. 
The  system  was  found  to  go  into  the  same  limit  cycle  (recorder  trace 
I-B-4.2).   From  these  two  recorder  tracess  the  limit  cycle  is  seen  to  be 
stable  with  an  angular  frequency  of  4.58  rad/sec  and  7.2  volts  amplitude, 
The  phase  portrait  for  the  system  appears  in  Figure  I-B-4.2. 
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I-B-5.   Comparison  of  the  Four  Analytic  Methods  In  Linit  Cycle  Predict 
for  the  System  Given  in  Figure  I-B. 


Method 

Frequency 

Amplitude 

rad/sec 

(volts) 

Mitrovic 

4.587 

6.: 

Root  Locus 

4.587 

6.21 

Describing  Function 

4.55 

7.33 

Analog  Simulation 

4.58 

7   2 
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CHAPTER  II 
DEAD  ZONE  IN  ERROR  CHANNEL 
II-A  Loop  Transfer  Function  with  No  Zeros. 

The  block  diagram  of  the  arbitrarily  chosen  system  which  is  to  be 

•«• 
analyzed  is  given  in 'Figure  II-A. 


e± 


Gc 


Figure  h-a 

Dead  Zone  in  Error  Channel  Loop 
Transfer  Function  with  No  Zeros 
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II-A-1.   Mitrovic's  Method. 

The  characteristic  Equation  of  the  given  system  is j, 

S4+lls3  +  26sN   I6S-*    fcOAl  =.  0   .  (ii-a-i.I) 

Assuming  the  last  two  coefficients  variable,  the  characteristic  equa- 
tion becomes, 

S44  IIS3  -H  Z^S2"  ■+    8,  5  -+60  =  O  (II-A-1. 2) 

where  the  parametric  equations  of  the  working  point  locus  are, 


and 


B  =  60  N  (II-A-1. 3) 

o 


B  =  16  ,  (II-A-1. 4) 

From  Appendix  A,  Mitrovic's  equations  for  B  and  B.  are, 

i30  =  -fcfcw^rtf)-*  \\u>f0t(<S)  +  uJ^fet&J  (11-A-1. 5) 


and 


6,  =  26^n^a(a)+  iitoffaltf)  +  w«3^4f») 


:i-a-i.6v 


Substitution  of  the  values  of  the  0  functions  from  Appendix  A  for 
j£   =0  gives  the  parametric  equations  of  the  stability  curve  as 

80-  2£>(J0r>2  -60n4  (II-A-1. 7) 

and 

6,  =  II  OO^2*     ,  (II-A-1. 8> 

At  the  intersection  of  the  M-point  locus  and  the  stability  curve> 
equations  II-A-1. 3  and  II-A-1. 4  and  Equations  II-A-1. 7  and  II-A-1. 8 
must  be  satisfied  simultaneously.   Thus  at  the  intersection, 

Jj6cvf-  16  (II-A-1. 9) 
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from  which  the  angular  frequency  of  the  limit  cycle  is  determined  to  be, 


l()     =  1.206  rad/sec. 


Also  at  the  intersection 

*'.,.%  4  _ 


26iun*- C*>n   -  £>0U      - 


(II-A-1.10) 


I       Substitution  of  ,  CO*2"  from  Equation  II-A-1.9  into  Equation  II-A-1.10 
gives  the  "average"  gain  of  the  nonlinearity  over  the  cycle. as, 


N    =0.595  , 
ave 


To  obtain  the  approximate  amplitude  of  the  limit  cycle,  reference  is 
made  to  Figure  II-A-1.1  which  is  a  sketch  o£  the  nonlinearity  static 
characteristic. 


Figure  II-A-1.1 
Limit  Cycle  Amplitude  Calculation  from  the 


Dead  Zone  Characteristic  Knowing  N 


ave 
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From  the  figure  the  following  equations  are  obtained, 

£     =  5  +  c  (II-A-1.11) 

**"  max 

c  =  a  cos  45°  =  0.707a  (II-A-1.12) 

a  «=  5  fr™  (II-A-1.13) 

Sin  B 

A  =  Tan"1  N  (II-A-1.14) 

ave 

and 

B  =  180  -  A  -  (180  -  45)  ,  (II-A-1.15) 

Entering  equation  II-A-1.14  with  the  "average"  gain  over  a  cycles 

A  =  Tan"1  0.595  =  30.8°  . 
Substituting  this  in  Equation  II-A-1.15  gives 

B  =  45°  -  30.8°  =  14.2°  . 
Substituting  for  A  and  B  in  Equation  II-A-1.13  gives, 


a  .  5  Sin  ??-*  =  10.469  . 
sin  14.2 

Substituting  for  "a"  in  Equation  II-A-1.12  gives 

c  =  0.707  (10.469)  =  7.402. 

And  finally  Equation  II-A-1.11  gives 

&  =   5  +  7.402  =  12.402  Volts 

max 

The  graphical  solution  of  the  intersection  is  sketched  in  Figure 
II-A-1.2. 

Mitrovic's  method  predicts  an  unstable  limit  cycle  by  the  following 
reasoning.   If  the  input  signal  is  small,  N  will  be  small  and  the  in- 
stantaneous M-point  will  be  in  the  stable  region  such  as  M  b   in  Figure 
II-A-1.2.   The  system  is  stable  and  the  oscillations  will  decrease  in  ampli- 
tude.  As  the  oscillations  decrease,  N  becomes  even  smaller  driving  the 
M-point  even  further  into  the  stable  region.   Consequently  the  M-point 
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STABILITY  CURVE  ISO) 


Figure  II-A-1.2 
Graphical  Solution  for  Dead  Zone  in  Error  Channel, 
Loop  Transfer  Function  with  No  Zeros 
can  never  go  into  the  unstable  region  and  all  oscillations  will  eventual- 
ly die  out.   If  the  input  signal  is  large,  N  is  relatively  large  and  the 
instantaneous  M-point  lies  in  the  unstable  region  such  as  M  .  The  system 

3 

is  unstable  and  the  oscillations  increase,  N  increases  toward  N  =  1.0, 
and  the  M-point  moves  farther  into  the  unstable  region.  Consequently, 
the  oscillations  increase  indefinitely  in  magnitude. 
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II-A-2.   Root  Locus  Method. 

Substitution  of  s  = ->  uj  into  the  characteristic  equation  given  in 
Equation  II-A-1.1  and  evaluating  the  higher  powers  of  J   gives,, 

604-Jilio3-26a>2  +jlkbo  +  (oON  -  0  .  (ii-a-2.1) 

By  requiring  that  both  the  real  and  imaginary  parts  of  Equation  II-A- 
2.1  go  to  zero  independently, 

4        2 
LU    -  2^60  +  60  N  -  0  (II-A-2. 2) 


and 


3 
-\\LQ   +  J/6o;=  0   •  (II-A-2. 3) 


Dividing  Equation  II-A-2. 3  by  J  (JO  s  gives 

11  UJ2"    -  16  =  0  (II-A-2. 4) 

from  which  the  angular  frequency  of  the  limit  cycle  is  obtained  as 

U->     =  1.206  rad/sec. 

2. 
Substitution  of   0J   from  Equation  II-A-2. 4  into  Equation  II-A-2.  k: 

and  solving  for  N    yields 
&      ave  J 

N    =  0.595  . 
ave 

To  determine  the  approximate  amplitude  of  the  limit  cycle s    Equations 
II-A-1.11  through  II-A-1.15  are  again  used.   Thuss 

£     =  12.402  Volts 

max 

A  graphical  presentation  of  the  root  locus  solution  is  given  in 
Figure  II-A-2.1. 

The  root  locus  method  predicts  an  unstable  limit  cycle  by  the  follow- 
ing reasoning.   If  the  input  signal  is  smalls  N  is  small  and  the  loop  gain 
is  small  placing  the  instantaneous  root  at  K\   in  Figure  II-A-2.1.   The 
system  is  stable  and  the  oscillations  decrease  in  magnitude.   As  the 
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Figure  II-A-2.1 
Root  Locus  Solution  for  Dead  Zone  in  Error  Channel, 
Loop  Transfer  Function  with  No  Zeros 
oscillations  decrease,  N  also  decreases,  driving  the  root  location  farther 
into  the  left  half  plane.   Consequently,  the  system  is  always  stable  and 
the  oscillations  eventually  die  out.   If  the  input  signal  is  large,  N  is 
also  relatively  large  making  the  loop  gain  high.  The  instantaneous  root 
location  is  at  f^   .  The  system  is  unstable  and  the  oscillations  increase. 
As  the  oscillations  increase,  N  increases  toward  N  =  1.  The  root  location 

moves  farther  into  the  right  half  plane.   Thus  the  system  is  always  un- 

i 

stable  and  the  oscillations  will  increase  indefinitely. 
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II-A-3.   Describing  Function  Method. 

The  describing  function  for  a  dead  zone  nonlinear! ty  is 


with 


R  = 


(II-A-3.1) 


(II-A-3.  2) 


where  b  is  the  total  minus  to  plus  dead  zone. 

The  linear  loop  transfer  function  for  the  given  system  in  frequency 
response  form  is 

3.7r 


G   Qoj)  =- 


Table  II-A-3.1  gives  the  magnitude  and  phase  of  Equation  II-A-3. 3 
for  the  different  frequencies  which  were  used  to  plot  the  linear  loop 
transfer  function  on  the  gain-phase  plane  of  Figure  II-A-3.1.   These 
values  are  obtained  from  a  Bode  diagram  plot  of  Equation  II-A-3. 3. 

TABLE    II-A-3.1 


(ii~a-3.  y. 


60 


60(0)  ^     1  j&Uul  a^ 


1.1 
1.2 
1.3 
1.4 


7.7 

6.2 
4.7 
3.5 
2.0 


-168 
-174 

-179.5 
-184 

-189.5 


TABLE 

II-A-3. 2 

R 

Go 

I/Gd 

l/Godb 

0 

1 

1 

0 

0.2 

0.746 

1.34 

2.55 

0.3 

0.624 

1.60 

4.10 

0.32 

0.61 

1.64 

4.29 

0.35 

0.577 

1.73 

4.78 

0.4 

0.504 

1.98 

5.94 

Table  II-A-3. 2  gives  the  values  of  1G/D,  in  the  area  of  interest, 

a 

These  values  were  obtained  from  equation  II-A-3.1. 
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By  interpolation  of  Figure  II-A-3.1,  the  intersection  and  thus  the 
limit  cycle  occurs  at 

(X)  =  1.21  rad/sec 
and 

R  =  0.34  # 
From  Equation  II-A-3.2,  the  amplitude  of  the  limit  cycle  is  then 

C  max    2R 
or 

£     =  14.71  Volts 

w  max 

The  describing  function  method  predicts  an  unstable  limit  cycle  by 
the  following  reasoning.   For  small  input  signals,  R  is  relatively  large. 
Assuming  the  -1/G  curve  to  be  the  locus  of  the  system  critical  point,,  then 
the  instantaneous  critical  point  is  at  R  in  Figure  II-A-3.1  and  the  sys- 

ex 

tern  is  stable  due  to  a  positive  phase  margin.   The  oscillations  will  de- 
crease causing  R  to  decrease,  thus  moving  the  critical  point  farther  up 
the  -180°  axis.   This  motion  of  the  critical  point  increases  the  positive 
phase  margin.   Thus  the  system  is  always  stable  for  small  signals  and  all 
oscillations  will  eventually  die  out.   For  large  input  signals  R  is  small 
placing  the  instantaneous  critical  point  at  R,  .   The  system  is  unstable  due 
to  a  negative  phase  margin  and  the  oscillations  increase.   As  the  oscilla- 
tions increase  R  becomes  even  smaller  moving  the  critical  point  down  the 
-180°  axis  and  increasing  the  negative  phase  margin.   Consequently ,  for 
large  signals  the  system  is  always  unstable  and  the  oscillations  increase 
indefinitely. 
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II-A-4  Analog  Computer  Simulation. 

The  analog  simulation  of  the  given  nonlinear  system  is  shown  in 
Figure  II-A-1.   The  coefficient  pot  settings  with  associated  resistances 
and  capacitances  for  real  time  and  real  magnitude  scaling  are  given  in 
Table  II-A-4. 1. 

TABLE  II-A-4. 1 


Pot  Setting 

Associated  Elements 

E. 

ai  =  100 

For  Input  Step  of  E,  Volts 
Rx  =  1  Meg,     Rfl  =  1  Meg 

a2  =  1.0 

R2  =  1  Megs     R_.  =  1  Meg 

a3  =  1.0 

R3  =  1  Meg,     Rf2  =  1  Meg 

a.  =  1.0 

4 

R,  =  1  Megs     C  .  =  1  uf 

a5  =  1.0 

R5  =  1  Megs    C   =  1  uf 

a,  =  0.6 
o 

R6  =  0.1  Megs   Cf2  =  1  uf 

a  =  1.0 

R?  =  0.5  Megs   Cf2  =  1  uf 

ag  =  1.0 

Rg  =  0.1  Megs   C   =  1  uf 

a9  =  0.8 

R.  =  0.1  Meg,   Cf  =  1  uf 

a1()  =  1.0 

R1Q  =  1  Megs    Cf4  =  1  uf 

an  =  1.0 

Ru   =   1  Megs    R   =  1  Meg 

a  =  0.5 

For  a  5  Volt  Dead  Zone  Each 

Initial  Setting 

Side  of  Zero 

To  accurately  set  the  back-bias  pots  "a"s  a  sine  wave  generator  was 
connected  to  the  input  of  the  dead  zone  simulator.   An  X=Y  recorder  was 
then  calibrated  for  5  volts  per  inch  on  each  axis.   The  input  and  output 
of  the  simulator  was  then  connected  to  the  X-Y  recorder  and  pots  "a"  were 
adjusted  until  both  break  voltages  of  the  static  characteristic  curve  were 
5  volts.   The  static  characteristic  curve  appears  as  Figure  II-A-4. 2. 
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The  system  was  first  excited  with  a  15  volt  step  input  signal  and 
the  system  was  found  to  be  stable  as  shown  in  recorder  trace  II-A-4.1. 
The  system  was  then  excited  with  a  17  volt  step  input  and  was  found  to  be 
unstable  as  shown  in  recorder  trace  II-A-4.2.   The  unstable  limit  cycle 
was  found  to  occur  for  an  input  step  of  17.07  and  is  shown  in  recorder 
trace  II-A-4.3.   From  recorder  trace  II-A-4.3  the  limit  cycle  is  seen  to 
have  an  angular  frequency  of  1.207  rad/sec  and  an  amplitude  of  15  volts. 

The  phase  protrait  of  the  system  is  shown  in  Figure  II-A-4.3. 
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CLEVELAND,  OHIO 


PHINTEO  IN  U.S.A. 


Recorder  Trace  II-A-4.1 

Dead  Zone  in  Error  Channel,  Loop  Transfer  Function  with  No  Zero's,  13  Vol; 

Step  Input 

SCALE:  £,  -  1  Volt/Line      £ '-  1  Volt/Line         X  -  1  Volt/1  - 

Time  -  1  Second/Division 
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vTHUIVIENTS  DIVISION  OF  CLEVITE  CORPORATION  CLEVELAND.  OHIO 
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BRUSH     INSTRUMENTS  division  of  clevite  corporatio 


Recorder  Trace   XX-A-4,2 

Dead  Zone   in  Error  Channel s    Loop  Transit  with 

No  Zeros,    17  Volt   Step   Input 


SCALE:    £-2  Volts/Line 

X-2.V6  Its/ Line 


£.  —  2  Volts/Lire   ' 
.  Ti^e  -  1  Second/Division 
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C  II  CHART  NO.  RA2921   32  BRUSH      INSTRUMEN 


Recorder  Tr&ce   II-A-4.3 

Limit  Cycle-,  Dead  Zone  in  Error  Channel,,  Loop  rransfeir  Function 
with  No  Zeros,  16.07  Volt  Step  Input 


SCALES: 


£  -  1  Volt/Line      £  -  1  Volt/Line        X  -  1  Volt/i. 

Time  -  1  Second/Division 
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II-A-5  Comparison  of  the  four  analytic  methods  in  limit  cycle  prediction 
for  the  system  of  Figure  II-A. 


Method 

Frequency 

Amplitude 

(rad/sec) 

(volts) 

Mitrovic 

1.206 

12.402 

Root  Locus 

1.206 

12,402 

Describing  Function 

1.21 

14.71 

Analog  Simulation 

1.207 

15.0 
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II-B  Loop  Transfer  Function  with  One  Zero. 


The  block  diagram  of  the  system  chosen  to  analyze  is  shown  in  Fig, 


II-B. 
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bW'  5(5+0.5)C5-MXSt^(548) 

• 

Figure  II-B 
Dead  Zone  in  Error  Channel  Loop 
Transfer  Function  with  One  Zero 
II-B-1.   Mitrovic's  Method.  , 

If  N  is  the  instantaneous  gain  of  the  nonlinearity  the  character- 
istic equation  for  the  system  is    t 


5rtl/.r54t3l.i*S3+^S1  +  (8+£'M)s  -r  20H  =0  . 


(II-B-1.1) 


Assuming  the  last  two  coefficients  to  be  variable,  the  characteris- 
tic equation  becomes, 


srt  j/.rs4  +  Z\SsZ  +  2ci  s1  +  6,  s  +  So  =-  0 


(II-B-1. 2) 


where 


8,  =    8  +£*N 


(II-B-1. 3) 


and 


G*  =      20N       . 
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(II-B-1. 4) 


From  Appendix  As  Mitrovie c  s  equations  £cx    B  and  B  are;, 

o  I 


80=  -[Wt»t?ftL&WSm\U8)  tll.W^W)  +  ^rffaCV)] 


(II-B-1, 


and 


6,  =  29w^i(K)+3/.r60n2^«)  +  l/.5iW«3^)+^V^W)  •    (II-B-1. 6] 

Substitution  of  the  values  of  the  0  functions  from  Appendix  A  for 

^"0   gives, 

60=  29  60,2-  JJ.iTu^  (II-B-1. 

and 

6,  =  3/,r60^-dO,4   ♦  (II-B-1. 8) 

The  sets  of  parametric  equations  II-B-1. 39  II-B-1. 4S  II-B-1. 7S    and 
II-B-1. 8  describe  the  M-point  locus  and  the  stability  curve  on  the  B 
versus  B1  plane.   Direct  solution  of  these  equations  for  the  intersection 
becomes  fairly  complicated;  consequently,,  a  straight  graphical  soluti  n  is 
performed. 

Table  II-B-1. 1  gives  the  values  of  B  and  B„  for  values  of   (jOa 
plotted  in  Figure  II-B-1. 1.   Elimination  of  N  from  Equations  II-B-1. 3  and 
II-B-1. 4  yields 

8,  =  8  +  -Js-  (11-B-1. 9 

which  is  a  straight  line  with  a  B  intercept  of  +  8  and  a  slope  of  +  4. 
The  intersection  of  the  two  curves  is  seen  to  be  approximately  at  ((L~   0.56 
rad/sec. 

To  accurately  determine  the  intersection  more  values  of  stability  curve 
are  calculated  and  the  value  of  B,  thus  determined  is  substituted  in  Equa- 
tion II-B-1. 9. 
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TABLE  II-B-1.1 


lq^ 

B 
o 

Bl 

0.5 

6.53 

7.81 

0.6 

8.95 

11.21 

0.8 

13.85 

19.75 

1.0 

17.5 

30.50 

1.2 

17.9 

43.23 

1.4 

12.6 

57.86 

1.588 

0 

73.0 

The  intersection  occurs  when  the  two  values  of  B  are  equal.   These  addi- 
tional calculations  appear  in  Table  II-B-1.2 


TABLE  II-B-1.2 


<On 

Bl 

B 
o 

B1 

o 

0.55 

9.43724 

8.72017 

5.74897 

0.567 

10.02355 

8.13459 

8.09419 

0.568 

10.05857 

8.15910 

8.23428 

0.57    | 10. 12679 

8.20816 

8.50716 

and 


By  interpolation  the  intersection  occurs  at 

6U0  =  0.5673  rad/sec 


B  =8.140 
o 

From  Equation  II-B-1.4,  the  average  nonlinear  gain  over  a  cycle  is 

N    =  0.407 
ave 
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Since  the  nonlinearity  is  the  same  as  the  one  of  Section  A-l  of 
this  chapter,  Equations  II-A-1.11  through  II-A-1.15  are  used  to  obtain 
the  approximate  amplitude  of  the  limit  cycle.   Substituting  N    into 
Equation  II-A-1.14  yields 

A  =  Tan"1  0.407  =  22.15°  , 
Substituting  A  into  Equation  II-A-1.15  yields 

B  =  22.85°  , 
Substituting  A  and  B  in  Equation  II-A-1.13  gives 

_   sin  22.15    ,  oc 
a  =  5  — : — 00  Q,  =  4.85  . 
sin  22.85 

Substituting  "a"  into  Equation  II-A-1.12  gives 

c  =  0.707  (4.85)  =  3.42  . 
Finally  Equation  II-A-1.11  gives 

£     =  5  +  3.42  =  8.42  volts  . 

'•"  max 

Applying  the  same  reasoning  as  given  in  Section  II-A-1  to  Figure 
II-B-1.2,  Mitrovic's  method  predicts  an  unstable  limit  cycle. 
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II-B-2.   Root  Locus  Method. 

Substitution  of   S-^W  into  Equation  IX-B-1.1  £;nd  Bimplii 
higher  order  J  terms  gives 

By  requiring  that  both  the.  real  and  imaginary  parts  :; 
II-B-2, 1  go  to  zero  independently  gives 

v-S"   k«  .  r-,,3 


0  0J -  -J  3  /.vTCO  <r)(foSU)(X)  -  0  ci-B-2.2; 


and 


ll.)Tio'-29io2+20N  -  o  •       (ii  b-: 

Dividing  Equation  II-B-2. 2  by  Ju)  and  solving  for  5N  gi^ 

A 
5  N  =  -  U)   +  31.5  U)    -8   ,  C-B-2.4) 

Substitution  of  Equation   II-B-2. 4   into   Eqwati        I  1.3   and   sisnplif 

ing  gives 

7.5  (JO    4     +     97  IX)   2        -    32   =  Q      , 

Solving  this  equation  for  WJ   gives 

\jj2     *   0.3219 

from  which 

60  =  0.567  rad/sec  » 

Substituting  the  value  of  tO '  back  into  Eq  a  .  :  II-B-2. 4  gi 

"average"  nonlinear  gain  over  a  cycle  as, 

N    -  0.4072  . 
ave 

To  find  the  approximate  amplitude  of  the  limitt  cycles  Equations 
II-A-1.11  through  II-A-1.15  are  again  used.   These  equations  yield 

£      =8.41  volts 
max 
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A  graphical  presentation  of  the  solution  is  sketched  in  Figure  II-B- 


2.1. 


Applying  the  same  reasoning  as  given  in  Section  II-A-2  to  Figure  II- 


B-2.1,  the  root  locus  method  predicts  an  unstable  limit  cycle. 


\ 


\        LINEAR  SYSTEM 
\  CLOSED  LOOP    ROOT 
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II-B-3.   Describing  Function  Method. 

The  describing  function  for  dead  zone  is  given  in   Equation  II-A-3.1. 

The  linear  loop  transfer  function  for  the  given  system  in  frequency 
response  form  is 


)(0.i2fi*jaj+i) 


(II-B-3.1) 


Table  II-B-3.1  gives  the  magnitude  and  phase  of  Equati  n  II-B-3.1 
for  the  different  frequencies  which  were  used  to  plot  the  linear  loop 
transfer  function  on  the  gain-phase  plane  of  Figure  II-B-3. 3.   These 

values  are  obtained  from  the  Bode  diagram  plot  of  Equation  U-E-3.1. 

TABLE  II-B-3.1 


<D 

I60")U*> 

/Gi^Ldeq] 

0.495 

+   10 

-172 

0.53 

+     9 

-176 

0.56 

+     8 

-180 

0.60 

+     7 

-184 

0.638 

+     6 

■=188 

Table  II-B-3. 2  gives  the  values  of  1/G  plotted  in  Figure  II-B-3.1. 

D 

These  values  are  obtained   from  Equation   II-A-3.1    for  various    values  of  R. 

TABLE    II-B-3. 2 


R 

Go 

l/Go 

1/6/3  cib 

0.4 

0.5044 

1.983 

5  .  ;-'  4 

0.45 

0.447 

2.237 

6.99 

0.48 

0.4132 

2.420 

7.48 

0.49 

0.4021 

2.4817 

7.91 

0.5 

0.3910 

2.558 

5.16 

By  interpolation  of  Figure  II-B-3.1  the  intersection  and  thus  the 
limit  cycle  occurs  at 
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CO    -   0.56  rad/sec 
and 

R  -  0.495  # 
From  Equation  II-A-3.2S  the  amplitude  of  the  limit  cycle  is  then 

£     =  om1?^ —  =  10.1  volts. 
max    2(0.495) 

Applying  the  same  reasoning  given  in  Section  II-A-3  to  Figure  XX-B- 

3.1,  the  describing  function  predicts  an  unstable  limit  cycle. 
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II-B-4  Analog  Computer  Simulation. 

The  analog  simulation  of  the  given  nonlinear  system  is  shown  in 
Figure  IX-A-4.1.   The  coefficient  pot  settings  with  associated  resist" 
ances  and  capacitances  for  real  time  and  real  magnitude  scaling  aire  gl 
in  Table  XX-A-4.1. 

TABLE  XX-A-4.1 


Pot    Setting 

Associated  Elements 

E. 

i 

ai  r      100 

For  Input   Step  of  E     volts 
R}   *   1  Meg,    Rfl    -   1  Meg 

a2   =   1.0 

R2   -   1  Meg,    R       «  1  Meg 

a3  -   1.0 

R.   =   1  Meg,    R  ,   =   1   Keg 

a.    =   0.4 
4 

R.    =0.4  Meg,   Cr,    ■  1   uf 
4                    OJ      rl 

a_   =  0.8 

j 

Rc   -  0.1  Meg,    Cfl   sluf 

a.,   =  1,0 

6 

R6  -  1  Megs          Rf3  ^   1  Meg 

a7   =   1.0 

R?    -   1   Meg,           R        »   1   Meg 

a8   .  0.5 

R8   -  0.1  Meg,    Cf2   ■  1  uf 

a9   *   1.0 

R9   =  0.5  Meg,    Cf2   =   1  uf 

a„n   =   1.0 
10 

R1Q   =   1   Meg,    Gf3   .    1   uf 

an  ■  i-o 

Rll   "  l  Me8>    Cf3  -  l  uf 

a12   =   1.0 

Rl2   -   1  Meg,    C£4  «   1   uf 

a^  =  0.5 

R13  =   1  Meg,    Cf4  -   1   u£ 

*14  e  !-0 

R14  =   1  Meg,    Cf5   -   1   uf 

a  =  0.05 

Initial   setting 

For  a  5  volt  Dead   Zone 
Each  Side  of  Zero 

The  back-bias  pots  "a"  were  accurately  set " In  the  same 
described  in  Section  IX-A-4.   The  static  characteristic  of  i 
linearity  Is  shown  in  Figure  XX-A-4.1. 
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The  system  was  first  excited  with  a  10.2  volt  step  input  and  v 
found  to  be  stable  as  shown  in  recorder  tra«e<  tI-B-4.1. 

was  then  excited  with  a  10.6  volt  step  input  and  foui 

as  shown  in  recorder  tr<=ce  XI->B-4.2a   The  unstable  U  It  cycle  - 

difficult  to  obtain  but  was  found  to  lie  between  step  inputs  of  10.3 

and  10.4  volts.   These  two  responses  are  shown  in  recorder  traces  It-E-4.3 

and  II-B-4.4.   The  average  of  the  two  traces  show  the  angular  frequen  , 

be  0.552  rad/sec  and  the  amplitude,,  using  the  first  two  peaks  of  es 

traces  to  be  9.9  volts. 

The  phase  portrait  of  the  system  is  shown  in  Figure  IX-B-402. 
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Recorder  Trace  II-B-4.1 

Dead  Zone  In  Error  Channel, 
Loop  Transfer  Function  with 
One  Zero,  10.2  Volt  Step 
Input.  i 

SCALES: 

t-    1  Volt/Line 
-   0.5  Volts/Line 

Time   r    5   Seconds/Division 


JSTRUIVIE1SJTS  division  ( 
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Recorder  Trace  II-B-4.2 

Dead  Zone  in  Error  Channel, 
Loop  Transfer  Function  with 
One  Zero,  10.6  Volt  Step  Input 

SCALES : 

%■    -  1  Volt/Line 
£  -  0.5  Volts/Line 

Time  -  5  Seconds/Division 


<\ 


RECORDER   MARF 
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Recorder  Trace  II-E-4.3 

LIMIT  CYCLE,  Dead  Zone  in  Error  Channel 9  Loop  Transfer  Function  with  One  Zero 

SCALES:   <£,-.!  volt /line        X  -  0.5  Volt /line    '  Time  -  1  Second/divia 

74  ~" 
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II-B-5  Comparison  of  the  Four  Analytic  Methods  in  Predicting  Limit 
Cycles  for  the  System  Shown  in  Figure  II-B. 


Method 


Angular  Frequency 
rad/sec 


T 


Mitrovic's  0.5673 

0.567 
Describing  Function  |  0.560 
Analog  Simulation      0.552 


Amplitude 

Volts 


8.42 

8.41 
10.1 
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CHAPTER  III 

IDEAL  RELAY  IN  ERROR  CHANNEL 

III-A  Loop  Transfer  Function  With  No  Zeros 

j 
The  system  arbitrarily  chosen  to  analyze  is  shown  in  the  block 

^ 

diagram  of  Figure  III^A. 


t    » 


N     * 

X 

-5* 

&k>  =  S(StO,5JfS^2j(5 +3X5+4) 


€b 


Figure  III-A 

Ideal  Relay  in  Error  Channel 

Loop  Transfer  Function  With  No  Zeros 

III-A-1  Mitrovic's  Method 

Assigning  N  as  the  instantaneous  gain  of  the  relay,  the  character- 
istic equation  can  be  written  as 

S^.vf  S4  +30.5"s3  +3752  -HZS  +\ZH   =*0...    '  (III-A-1. 1) 


Assuming  the  last  two  coefficients  variable,  the  characteristic 
equation  becomes, 

A    —  ^    3 


S^9.s-s>30.r£ +37s2-+  6,5+Bo-O, 


where 


B  =  12  N 
o 


(III-A-1. 2) 


(III-A-1. 3) 
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and 

B-  =  12  .  (III-A-1. 

1 

From  Appendix  A,  Mitrcvic''s  equation  for  B  and  B,,  a«9 

o      i 


and 


6,  -  ^7oi,^^J+3o1^/^^J^/6'^V/^f^V^  •       (iii-a-i 

Substitution  of  the  values  of  the  0  functions  from  Appendix  A  for 
«   =0  gives  the  parametric  equations  of  the  stability  eurve  as, 

Bo  =  37  6l\,     -   9.5  UJ^         (III-A-1. 7) 

and 

B.  -  30.5  UJn2    "  ^n    «       (III-A-1. 8) 

To  find  the  intersection  between  the  R- point  locus  (Equations  III-A- 
1.3  and  IIX-A-1.4)  and  the  stability  curve  (Equations  III-A-1. 7  and  III- 
A-1.  8) ,  the  graphical  approach  is  chosen.   Table  III-A-1. 1  gives  the 

values  of  B  and  B,  for  various  (jQn  that  are  plotted  in  Figure  III-A-1 ,1, 
o      1 

The  M-point  locus  is  seen  to  be  a  straight  line  parallel  to  the  B  axis 

o 

at  B,  =  12.   From  Figure  III-A-1. 1,  the  intersection 'is  seer:  to  occur  at 
1 

d0^=  0.62. 


TABLE  III-A-1.1 


£0n 

B 
o 

h 

0.5 

8.65 

7.56 

0.6 

12.08 

10.85 

0.7 

15.85 

15.70 

0.8 

19.8 

19.2 

1.0 

27.5 

29.5 

1.2 

33.6 

41.8 

1.4 

36.0 

56.0 

1.6 

32.5 

71.4 

1.8 

20.3 

88.2 

1.974 

0 

103.6 

To  accurately  determine  the  intersection  more  values  of  B  and 

o 

Bj  are  calculated  for  60^  close  to  0.62.  These  values  appear  in  Table 
III-A-1.2.  By  interpolation  from  Table  III-A-1.2  the  intersection  and 
thus  the  limit  cycle  occurs  at 

^Or)  =  0.6315  rad/sec 
and 


B  =  13.24 
o 


TABLE    III-A-1.2 

000 

B 
o 

B 
1 

0.62 

12.8191 

11.5764 

0.63 

13.1891 

11.9480 

0.631 

13.2259 

11.9854 

0.632 

13.2631 

12,0229 

The  "average"  nonlinear  gain  over  a  cycle  is  determined  from 


79 


Equation  III-A-1.3  as, 


N    =  1.103 
ave 


From  this  the  approximate  amplitude  of  the  limit  cycle  is 

P      =  4,53  volts 
^  max 

Note  that  the  M-point  of  the  linear  system  in  Figure  Ijfl-A-I.l  shows 

that  the  linear  system  is  stable.   It  is  well  known  that  an  ideal  relay 

gives  a  limit  cycle  for  a  linear  system  of  order  three  or  above  and  the 

occurance  of  a  limit  cycle  is  independent  of  linear  system  stability. 

This  is' easily  explainable  on  the  B  versus  B  plane.   The  ideal  relay 

o         1 

has  an  infinite  maximum  instantaneous  gain  whereas  the  nonlinear ity  used 

in  Chapters  I  and  II  have  a  maximum  instantaneous  gain  of  one.   Thus  on 

the  B  versus  B,  plane  the  M-point  locus  is  unlimited  for  the  ideal  re- 
o         1 

lay  but  does  have  a  maximum  value  for  the  other  non linearities.   Conse- 
quently, an  infinite  straight  line  segment  must  necessarily  intersect  the 
stability  curve  at  some  point  producing  a  limit  cycle. 

Mitrovic's  method  predicts  a  stable  limit  by  the  same  reasoning  used 
in  Chapter  I,  Section  I-A-l.   This  same  argument  can  be  applied  to  the 
ideal  relay  since  this  nonlinearity  can  be  envisioned  as  saturation  with 
an  infinite  gain  in  the  linear  region. 
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III-A-2   Root  Locus  Method. 

Substituting  $  -Ju)  into  the  characteristic  Equasti 
evaluating  the  higher  order  powers  of  J  yiel 

Requiring  that  the  real  And   imaginary  parts 

go  to  zero  independently  gives s 

and 

Dividing  Equation   IIX-A-2.2   by   ^(x   and   ?    1  ring     i    i  U)       g: 

2 
U)        m  0.39E 6, 

from  which 

60  =  0.6314     rad/sec     , 

2 
Substituting  ^0        into  Equation   III-A-2. 3   and   so]         =    for  N   ; 

N  =   1.1032 

ave 

From  this  "average"  nonlinear ity  gain  over  a  cyeles  the  approxi  .  ' i 
amplitude  is  determined  as, 

£     =4.533  volts 

max 

The  graphical  presentation  of  the  solution  by  i 
is  shown  in  Figure  III-A-2. 1„ 

The  root  locus  method  predicts  a  stable  limit  cycle  t         ;  the 
same  arguments  to  Figure  III-A-2. I  that  are  gi1  a  i  is  CI  £p5  as 
I-A-2. 
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LIMIT  CYCLIf 
W  =  0.63)4   rbA/Scc 


Root  Locus   Solution   foi    Ideal    Re  Ln  Em    _-    _.-- 

Loop  Transfer   Function  wi 
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IXT-A-3  Describing  Function  Method. 

The  describing  function  for  an  ideal  relay  is 

.0 


/J* 


(III-A-3.1) 


ff\o./ 


The  linear  loop  transfer  function  for  the  given  system  In  frequency 
response  form  is, 


G  (Jrt)  - 


^(2JuH|)(O,£Juj+|)(O.333^+/X0'^^  +  0 


(III-A-3.2) 


Table  III-A-3.1  gives  the  magnitude  and  phase  of  Equation  III-A-3.2 
for  the  different  frequencies  which  were  used  to  plot  the  linear  loop 
transfer  function  on  the  gain-phase  plane  of  Figure  III-A-3.1.   These 
values  are  obtained  from  a  Bode  diagram  plot  of  Equation  III-A-3.2. 


IABLE    III-A-3 

1 

CO 

\G(^)\db 

/GfrU  <i«* 

0.527 

2.0 

-167.5 

0.565 

1.0 

-172 

0.605 

0 

-  176.5 

0.625 

-0.5 

-173. 5 

0.648 

-1.0 

-181 

0.69 

-2.0 

•185.5 

0.735 

-3.0 

-190 

Table  III-A-3.2  gives  the  values  of  1/G  ir.  the  area  cf  Interest  as 
obtained  from  Equation  III-A-3.1  for  various  values  of  error. 

TABLE  III-A-3.2 


C.  r<><3  y 

Gn 

//An 

//fto    <Lb 

1.273 

' 0.7854 

-2.10 

6 

1.061 

0.942 

-0.52 

5.8 

1.098 

0.911 



-0.80 
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By  interpolation  of  Figure  XXI-A-3.1S  the  intersection  and  thus 

the  limit  cycle  occurs  at 

60=  0.635  rad/sec 

and 

C,     =  5«b  volts  o 
max 

The  describing  function  method  predicts  a  stable  limit  cycle  by 
applying  the  same  arguments  in  Section  I-A-3  of  Chapter  I  to  Figure  XII- 
A-3.1 
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I1I-A-4  Analog  Computer  Simulation. 

The  analog  simulation  of  the  given  nonlinear  system  Is  sh  . 
Figure  III-A-4.1.   The  coefficient  pot  settings  with  assc 
ances  and  capacitances  for  real  time  and  real  magnitt  cde  scaling  -are  g. 
in  Table  III-A-4.1. 

TABLE  III-A-4.1 


Pot  Setting 

Associated  Elements 

Ei 

a  *  

1    100 

For  Input  Step  of  Ev  Volt 
R1  -  I  Meg,   Rf]  .■=  ^  1  Meg 

a?  =  1.0 

R2  =  1  Meg9   R   «  1  Meg 

a   -  (See  discussion) 

R.3  -  0.1  Meg,  Rf2  =   1  Meg 

a.  =  0.4 
4 

R4  =  0.1  Meg,   Cfl  -  I  of 

a^  .  0.5 

3 

R5  «  1  Meg,   Cf.  ■  1  uf 

a,  -  0.3 
6 

R.  =0.1  Meg,   C  -  -  1  uf 

a_  =  1.0 
1 

R„  =  0.5  Meg,   C,,  =  1  uf 
i                                          1  & 

ag  =  1.0 

Rg  -  1  Meg,    Cf3  -  1  u 

a9  =  0.3 

R9  =  0.1  Meg,  Cf3  =  1  uf 

aio  *  1-° 

R1Q  -  1  Meg,  Cf4  =  1  uf 

aH  ~  °'4 

Rn  -  0.1  Meg,  Gf4  -  1  uf 

ai2  s  l 

R--  =  1  Meg,   Crc  >  1  uf 
12       os    f5 

a  s  0.05  (Initial 

For  5  Volt  Voltage 

setting) 
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The  simulation  of  the  ideal  relay  is   seen  to  be  a  simple  saturation 
simulation.   However,,  the  gain  of  the  linear  portion  is  wade  very  hi 
by  the  insertion  of  coefficient  pot  a  .   To  set  the  back-bias  pots 
and  pot  a. j,  a  sine  wave  generator  was  connected  to  the  input  of  the  r« 
lay  simulator.   The  output  was  monitored  and  the  back-bias  pots  "a"  were 
then  adjusted  until  the  output  was  limited  at  +  5  volts.   The  output  and 
input  was  then  connected  to  an  X  -  Y  recorder  to  obtain  a  static  chars 
istic  curve.   Pot  a  was  then  adjusted  until  the  slope  p 

tion  was  as  nearly  vertical  as  possible  without  the  amplifier  overload 
The  static  characteristic  curve  appears  as  Figure  IIX-A-4.2. 

The  system  was  first  excited  with  a  1  volt  input  signal  and  was  found 
to  build  up  to  the  limit  cycle  as  shown  in  recorder  trace  IXX-A-4.1.   Fne 
system  was  then  excited  with  a  10  volt  step  input  and  was  found  t :■  g  : 
the  same  limit  cycle  as  shown  in  recorder-  trace  XXX-A-4.2.   From  these 
recorder  traces  the  limit  cycle  is  seen  to  be  stable  with  an  angular  fre- 
quency of  0.628  rad/sec  and  an  amplitude  of  5.7  volts. 

The  phase  portrait  of  the  system  is  shown  in  Figure  IIX-A-4.3. 


39 




I  . 


■         I   ":     ri 

1 1 1 ;  'i rn  i  . 

,i  ,i,i  1 1 1 1 
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Recorder  Trace  IXI-A-4.1 

Ideal  Relay  in  Error  Channel,  Loop  Transfer  Function  with  No  Zeros, 
1  Volt  Step  Input  ■  * 


SCALES : 

£  -  0,5  Volts/  Line 


£    -  0.5  Volts/Li/f.e       X  -.   0.5  Volts/1 

Time  -  5  Seconds/Division 
1  Second/Division 
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■NTS  DIVISION  OF  CLEVITE  CORPORATION  CLEVELAND.  OHIO 


PRINTED  IN  US  A 


Recorder  Trace  E-iP-A-4.2 

?nef  ifc?'  ^  Err°r  Cha"nel>  ^°P  mnsfer  Kinetics  it]      Zeros. 
10  Volt  Step  Input. 


SCALES.: 
£  -  0,5  Volts/Line 


£  -  0.5  Volts/Lii        x  -  0.5  v  3 
Time  -  5  Seconds /Division 
1  Second/Division 
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IXX-A-5  Comparison  of  the  Four  Methods  in  Predicting  Limit  Cycles 

the  System  of  Figure  III-A. 


Method 

Angular  F:e,  -  - 

••-■  p  - 

rad/sec 

Volts 

Mitrovic8  s 

0.6315 

4.53 

Root  Locus 

0.6314 

4.!  5 

Describing  Functions 

D.635 

5.8 

Analog  Simulation 

0.628 

5.7 
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III-B  Loop  Transfer  Function  with  One  Zero. 

The  system  to  analyze  is  shown  in  the  block  diagram  of  Figure 
III-B. 


Qe 


Figure  III-B 
Ideal  Relay  in  Error  Channel  Loop 
Transfer  Function  with  One  Zero 
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III-B-1.   Mitrovic's  Method. 

Assigning  N  as  the  instantaneous  gain  of  the  relay,  the  character- 
istic Equation  can  be  written  as 

(iii-b-1.  i; 

Assuming  the  last  two  coefficients  are  variable,  the  characteristic 
equation  becomes, 


S*+  3lSff+3i0s4  4 /240S3+ 19345*48,5  4  60=  0 


(III-B-1 


where 


and 


B      =-4000  N  (III-B-1.  --, 

o 


B     -   1024  +  8000  N  (III-B-1. 

From  Appendix  A,    Mitrovic's  equations    for   B     and   B-    are, 

+  Z\U)f<foM+    LOri^fo^J]        (III-B-1. 5) 
and 

+  3!uJ^^)-f^^^J  •      Cm-B-] 

Substituting  the  values  of  the  functions  from  Appendix  A  for  %   -  0 
gives  the.  parametric  equations  of  the  stability  curve  ass 

80-  /9S4  uS-VOUinU^b  (III-B-1.7) 

and 

8,=     ia^0  06^-   3/oJk)4  (111-6=1.3) 

To    find  the    intersection  between  the  M-poJ  (Equations    . 
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and  III-B-1.4)  and  the  stability  curve,  the  graphical  approach  is  chosen, 

Table  III-B-1.1  gives  the  values  of  B  and  B,  for  the  stability  curve  as 

o      1  J 

plotted  in  Figure  III-B-1.1. 

TABLE  III-B-1.1 


Wn 

0 

Bi 

0.5 

227 

308 

0.8 

1143 

780 

1.0 

1675 

1209 

1.2 

2219 

1722 

1.4 

2708 

2311 

1.6 

3067 

2962 

1.8 

3200 

3695 

2.0 

3040 

4464 

2.2 

2454 

5275 

2.3 

1968 

5692 

2.5 

544 

6540 

To  obtain  the  M-point  locus  solve 
substitute  into  Equation  III-B-1.3. 


ion  III-B-1.4  for  N  and 


B  =  0.5  B.,  -  512 


(III-B-1.9) 


Equation  III-B-1.9  is  seen  to  be  a  straight  line  with  a  slope  of  0.5  an 
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a  EL  intercept  of  1024. 

It  is  seen  from  Figure  IXI-B-1.1  that  the  intersection  occurs 

roughly    AJp.  =  2,24.   To  more  accurately  determine  the  intersection, 

values  of  B  and  B,  are  calculated  for   LO <a   close  to  2.24.   The  value 
o      1  n 

of  B1  is  then  substituted  into  Equation  IXI-B-1.9  and  the  intersection 
occurs  when  the  two  B's  are  equal.   These  additional  values  are  given 
in  Table  XII-B-1.2. 


TABLE  III-B-1.2 


60^i 

B   (Stability 
curve) 

B- 

JL 

B  (M-point 

Locus) 

' 

2.24 

2266 

5439 

2208 

2.249 

2233.5 

5478.8 

2227.4 

2.25 

2228.7 

5483.0 

2229.5 

Thus  from  Table  III-B-1.2,  the  frequency  of  the  limit  cycle  is 


£dn  =2.25  rad/se 


c  p 


Substituting  B  at    LAA   =  2.25  into  Equation  III-B-1.2,  the  "average* 
gain  over  a  cycle  is  determined  as 


K 


2229 


=  0.5573  , 


ave     4000 
from  which  the  approximate  amplitude  is  determined  as 


max    0.5573 


=  8.973  volts, 
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Mitrovic9 s  method  predicts  a  stable  limit  cycle  by  the  same  argu- 
ments applied  to  Figure  III-B-1.1  as  were  given  in  Chapter  X„  Section 
I-A-l. 
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III-B-2   Root  Locus  Method. 

Substituting  s  =  j(jj  into  the  characterie 
and  evaluating  the  higher  powers  of  J   yields,, 

-u>Sj3lu/+3IOu^-J/a40oo3H9?4<o:L 

Requiring  that  the  real  and  imaginary  parts  of  £  B-2.1 

go  to  zero  independently,,  gives 

and 

Dividing  Equation  III-B-2.3  by  Jui  and  sol  ring  this  ft  r  4000  is 
g  ive  s  s 

406&N  =  ClOuf-m-lf.rto4*  Cm-" 

Combining  Equations  III-B-2.3  and  III-B-2.4  gives 

2 

To  solve  for  LO      from  Equation  IIILB-2.5S  a  digital  lomptater  pj 

gram  was  used  and  gave  the  roots, 


2 

-   5.062 

2 
60 

-  -   0.55 

2 
U) 

39.55 

Only  the  first  value  is  of  use  and  gives 

(jj   =2.50  rad/sec. 

2 

Substituting  (j^       into  equation  III-B-2.4  gives 

of  the  relay  over  a  cycle  as 
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from  which  the  approximate  amplitude  of  the  limit  cycle  is  determined  to 
be 

£     =  8.972  volts 

max 

A  graphical  presentation  of  the  solution  is  sketched  in  Figure 
1 

III-B-2.1.  The  root' locus  method  predicts  a  stable, limit  cycle  by 
applying  to  Figure  III-B-2.1  the  same  arguments  given  in  Chapter  I, 
Section  II-A-2. 


■LIMIT  CYCLE 
*W=0.5f74 


+<T 


-.ito 


Figure  III-B-2.1 
Root  Locus  Solution  for  Ideal  Relay  in  Error  Channel, 
Loop  Transfer  Function  with  One  Zero. 
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IXX-B-3  Describing  Function  Method. 

The  describing  function  for  the  ideal  relay  Is  given 
III-A-3.1.   The  linear  loop  transfer  function  for  the  gi     system  in 
frequency  response  form  is 

rim-     .       .  3^/(2jujfi) 


(III-B-  3.1; 


Table  III-B-3.1  gives  the  magnitude  and  phase  of  Eq    I  >n  XXX-B-3.1 
for  the  different  frequencies  which  were  used  tc      th<  linear  loop 
transfer  function  on  the  gain-phase  curve  of  Figure  III-B=3.1.   These 
values  are  obtained  from  a  Bode  diagram  plot  of  Equation  IXX-B-3. 1. 


TABLE    III-B-3.1 

uo 

I  GOto)  1   <U 

jGQuSl    ^e9 

1.97 

7 

-168 

2.1 

6 

-174,5 

2.24 

5 

-180 

2.38 

4 

-183 

2.52 

3 

-190 

! 

Table  III-B-3.2  gives  the  values  of  1/G  in  the  area   : 
obtained  from  Equation  III -A- 3.1  for  various  values  of  err.  » 


TABLE  III-B-3.2 


E 
max 

GD 

1/G, 

D 

1/G^dfe 

u 

8 

0.796 

1.2  57 

1 .  99 

10 

0.634 

1.571 

3.: 

11 

0.579 

1.729 

4.75 

11.3 

0.563 

1.775 

4.98 

12 

0.531 

1 . 8  • 

103 


By  interpolation  of  Figure  IXI-B-3.1.,  the  intersection  and  thus 
the  limit  cycle  occurs  at 

60  s  2,24  tad/sec, 

and 

E    -  11.31  volts  , 
max 

The  describing  function  method  predicts  a  stable  limit  cycle  by 
applying  to  Figure  III-B-3.1  the  same  arguments  given  in  Chapter  Is 
Section  I-A-3. 
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III-B-4  Analog  Computer  Simulation, 

The  analog  simulation  of  the  given  nonlinear  system  is  shown  in 
Figure  IXX-B-4.1.   The  coefficient  pot  settings  with  associated  re- 
sistances and  capacitances  for  real  time  and  real  magnitude  scaling  are 
given  in  Table  III-B-4.1. 

TABLE  III-B-4.1. 


Pot  Setting 

Associated  Elements 

E. 

1 

For  a  Step  Input  of  E  Volts 
R  =  1  Meg,     R  „  =  1  Meg 

ai  "  100 

a2  =  1.0 

R2  =  1  Meg,     Rfl  =  1  Meg 

a   =  (See  discussion) 

R3  =  0.1  Meg,   R  2  =  1  Meg 

a.  =  0,5 
4 

R4  =  1  Megs    Cfl  -  1  uf 

a   =  1.0 

j 

R5  =  1  Meg,     Cfl  -  1  uf 

a,    m   1.0 
6 

R6  =  1  Meg,    R£3  =  1  Meg 

a   =  1.0 

R?  -  1  Meg,     R    -  1  Meg 

a8  -  1.0 

RQ  =  0.1  Meg,   Cro  =  1  uf 

o                     f  I. 

a9  =  1.0 

Rg  =  0.5  Megs  Cf«  ~  l   uf 

«10  =  i.o 

R1Q  =  0.1  Megs  Gf3  -  1  uf 

an  =  0.4 

R.,  =0.1  Meg,  Cro  =  1  uf 

11        OJ   t3 

a12  =  U0 

R12  =  0.1  Meg,  Cf4  =  1  uf 

a13=0.8 

R13  =  0.1  Meg,  Cf4  =  1  uf 

au  =  0.8 

R14  =  1  Meg,    Gf  =  0.1  uf 

a,_  =  0.8 

15 

R15  -  0.5  Meg,  Cf5  =0.1  uf 

a16  -  1.0 

R16  =  1  Meg»   Cf6  "  l   uf 

a17  '  U0 

R17  =  1  Meg,    Rf4  -  1  Meg 

a  ■  0.05  (Initial 
Setting) 

For  5  Volt  relay  voltage 
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The  back-biasing  pots  "a"  and  pot  a_  were  set  in  the  same  way 
described  in  Section  III-A-4.   The  static  characteristic  curve  is  also 
the  same  as  shown  in  Figure  III-A-4. 1. 

The  system  was  first  excited  with  a  5  volt  step  input  signal  and 
was  found  to  build  up  to  the  limit  cycle  shown  in  recorder  trace  III- 
B-4.1.   The  system  was  then  excited  with  a  20  volt  step  signal  and 
found  to  go  into  the  same  limit  cycle  as  shown  in  recorder  trace  III-B- 
4.2.   From  these  two  recorder  traces  the  limit  cycle  is  seen  to  be  stable 
with  angular  frequency  of  2.22  rad/sec  and  an  amplitude  of  11.7  Volts. 

The  phase  portrait  of  the  system  is  shown  in  Figure  III-B-4.2, 
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III-B-5.   Comparison  of  the  Four  Methods  in  Pred  cting  Lit      les  £0.1 
the  System  Shown  in  Figure  III-B. 


1 

Method 

Angu la r  Fre que n zy 

Amp , 

rad/sec 

'V ;  . ';.  S 

Mitrovic" s 

2/4 

8.973 

Root  Locus 

2 .  2 ;  D 

8.972 

Describing  Function 

2,24 

11.31 

Analog  Simulation 

2.22 

11.7 

:,: 


CHAPTER  IV 
RELAY  WITH  DEAD  ZONE  IN  ERROR  CHANNEL 
IV-A.   Loop  Transfer  Function  with  No  Zeros. 
»    The  block  diagram  of  system  chosen  is  shown  in  Figure  IV-A. 
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+  5 
-5 

X 

f7£k-           SO 

Sb   J 
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+r  e 

1    > 

U    J    SCS+05£+tf5f2XS+4X*fc) 

Figure   IV- 
Relay  with  Dead   Zone   in 
Transfer  Function  with 

A, 

Error  Channel  Loop 
No   Zeros 

IV-A-1  Mitrovic's  Method 

The  characteristic  equation  of  the  given  system  is, 

S4,  +  l5.5"5*+77.,5-S4  +  IS&s2  +  124  s2  +  3ZS  +5"0N   -  0. 


(IV-Arl.l) 


Assuming  the  last  two  coefficients  variable,  the  characteristic 
equation  becomes 


5    t  <i75"5'r+77.5V+l5\rs3-tl2  4  s2"  +  3,  s'  +  6o=0 


(IV-A-1.2) 
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where 


and 


B  =  50  N  (IV-A-1.3) 

o 


Bt  =  32  (IV-A-1.4) 


From  Appendix  A,  Mitrovic's  equations  for  B  and  B  are, 


(IV-A-1.5) 


and 

Substitution  of  the  values  of  the  0   functions  from  Appendix  A  for 

^  =  0  gives  the  parametric  equations  of  the  stability  curve  as, 

BQ  =  124  u;0    -   77.5  U)„    +  U)n  (IV-A-1.7) 

and 

B-  =  155  U)n2      -   15.5  U0n   4  9  (IV-A-1.8) 

Table  IV-A-1.1  gives  the  points  of  stability  curve  plotted  in  Figure 

IV-A-1.1.   These  points  are  calculated  from  Equation  IV-A-1.7  and  IV-A-1.8 

for  various  LO^  . 

From  Equations  IV-A-1.3  and  IV-A-1.4,  the  M-point  locus  is  seen  to 

be  a  straight  line  parallel  to  the  B  axis  with  a  maximum  of  3  =  50. 

o  o 

From  Figure  IV-A-1.1,  the  intersection  of  the  M-point  locus  with  the 
stability  curve  is  seen  to  occur  at  roughly  Ct\-   0.46  rad/sec.   To  more 
accurately  determine  the  intersection  more  points  in  the  area  of  the 
intersection  are  calculated  in  Table  IV-A-1.2. 
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TABLE  IV- A- 1.1 


UOn 

Bl 

B 

o 

0.2 

6.18 

4.76 

0.4 

24.0 

17.9 

0.5 

37.8 

26.3 

0.6 

53.8 

34.6 

0.7 

72.2 

42.4 

0.8 

92.9 

48.1 

0.9 

115.4 

50.1 

1.0 

139.5 

47.5 

1.1 

164.8 

38.3 

1.2 

190.9 

25.5 

TABLE  IV-A-1.2 


0>n 

B 

o 

h 

0.45 

21.94 

30.75 

0.46 

22.78 

31.59 

0.461 

22.86 

32.24 

0.47 

23.62 

33.48 

, . 

By  interpolation  of  Table  IV-A-1.2  using  Equation  IV-A-1.4,  the 
intersection  and  thus  the  limit  cycle  occurs  at 


and 


tOr\    ~  0.4607   rad/sec 


B     -  22.84    # 
o 
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From  Equation  IV-A-1.3,  the  "average"  nonlinear  gain  over  a   cycle 
is  determined  as, 

N,  =  0.457 
from  which  the  approximate  amplitude  of  the  limit  cycle  is  determined  to 


be 


£     =  10.95  Volts  , 
w  max 


Mitrovic's  method  predicts  a  stable  limit  cycle  by  the  following 
reasoning.   The  input  signal  must  be  larger  than  the  dead  zone  for  the 
system  to  respond  at  all.   Once  the  input  signal  becomes  larger  than  the 
dead  zone,  the  characteristic  of  the  relay  is  for  practical  purposes 
similar  to  saturation.   Thus  the  same  arguments  given  in  Chapter  I> 
Section  I-A-l  can  be  applied  to  Figure  IV-A-1.1. 
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IV- A- 1.2   Root  Locus  Method. 

Substituting  S    -   JOJ   in  the  characteristic  EquatJ      A- 1 . 1 
evaluating  the  higher  powers  of  J  gives, 

-(0^4  J  \£.5oJ +  17Slo/I-)ISSujz-  tltLO1 +  )3lti  +JTON  -0. 

(XV-A-2.1) 

Requiring  that  the  real  and  imaginary  parts  of  Equation  IV-A-2.1  go 
to  zero  independently,  yields  the  following  two  equations,, 

-60  S  77, 6"^  -WUJZ  +SON   -£>  (IV-A-2.2) 

and 

^iJVvTfc/  -MSSco3  +  J31^>  =  O  .  (IV-A-2.3) 

Dividing  Equation  IV-A-2.3  by  ](j)     and  using  the  quadratic  formula 


2 
to  solve  for  60   yields, 


2. 
to  =  0.210899 


as  the  only  meaningful  root.   Thus  the  frequency  of  the  limit  cycle  is 

iO   =  0.4592  rad/sec. 

2 
Substituting  60   back  into  Equation  IV-A-2.2  and  solving  for  N 

as  the  "average"  gain  of  the  nonlinaaritys  gives 

N    -  0.4543 
ave 

from  which  the  approximate  amplitude  is  found  to  be 

£  =  11.01  volts  . 
m 

A  graphical  presentation  of  the  root  locus  solution  is  shown  in  Figure 
IV-A-2.1. 

The  root  locus  method  predicts  a  stable  limit  cycle  by  the  following 
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arguments.   For  an  input  signal  larger  than  the  dead  zone,  the  relay 
characteristic  is  similar  to  saturation  and  the  arguments  given  in 
Chapter  I,  Section  I-A-2  can  be  applied  to  Figure  IV-A-2.1. 


uneaR  systfm  CLOSED  LOOP  ROOT 


LIMIT  CYCLE 


Figure  IV-A-2.1. 
Root  Locus  Solution  for  Relay  with  Dead  Zone 
in  Error  Channel,  Loop  Transfer  Function  with 
No  Zeros 
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IV-A-3.   Describing  Function  Method. 

The  describing  function  for  a  relay  with  dead  zone  is 


GD  = 


4V 


<-ferX    zs! 


(IV-A-3. 1) 


where 


V  is  the  signal  controlled  by  the  relay 
^^is  the  relay  input  signal 
d-J  is  the  total  dead  zone  minus  to  plus. 
The  linear  loop  transfer  function  for  the  given  system  in  frequency 
response  is 


G  (J60)- 


/,£*££ 


Jw(2Ju)+/J(Jio+i)(a5jw+i))j3,25'j(o+i)(o.iarj£o+0    .,      (iv-a-3.2) 


Table  IV-A-3. 1  gives  the  magnitude  and  phase  of  Equation  IV-A-3.2  to: 
the  various  frequencies  that  are  used  to  plot  the  linear  loop  transfer 
function  on  the  gain-phase  plane  of  Figure  IV-A-3. 1.   These  values  are 
obtained  from  a  Bode  diagram  plot  of  Equation  IV-A-3.2. 

TABLE  IV-A-3. 1 


to 

0.4 

9 

-172 

0.43 

8 

-177.5 

0.455 

7.5 

-180 

0.463 

7 

-182 

0.49 

6 

-186.5 

Table  IV-A-3.2  gives  the  values  of  1/G  which  are  plotted  in 
Figure  IV-A-3. 1.   These  values  are  obtained  from  Equation  IV-A-3. 1, 
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TABLE  IV -A-  3, 


&CVUX/<. 

GD 

G 

10 

0.55*3 

■ 

12 

0.4823 

2.0 

6,33 

14 

0.4247 

2.  354 

.44 

15 

0.4026 

2.484 

• 

Inspection  of  Figure  IV-A-3.1  shows 
limit  cycle  to  occur  at., 

CO   =  0.4 5 . 

and 

cLa*l/  =  14. 1  •  Its 
Ihe  describing  function  method  predicts  a  stab 
the  following  arguments.   For  an  input  signs   greater  I 
zone,  the  relay  characteristic  is  similar  to  satura  -    - 
ments  given  in  Chapter  I>  Section  t-A-3  can  be  applied  I   i  guts    Iv-a-j.  "'. . 
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IV- A -4  Analog  Computer  Simulation 

The  analog  computer  simulation  of  the  gl vein  i.s 

shewn  in  Figure  IV-A-4.1.  The  coefficient  pot  sett Inge 

resistances  and  capacitances  for  real  t:  i  lagnitudle 

are  given  in  Table  IV-A-4.1. 

I ABLE  IV-A-4.1 


Pot    Settings 


i 

a;   "     100 

a ,    =   1.0 


a      -    1.0 


a,    -   1.0 

4 

a_     r-    \.Q 
3 

a..    =0.5 


=    i.O 


a0    =   1.0 

8 

a„  =  1.0 


a.,    =   i.O 


all 


0. 


a12  .0.4 


13 


14 


=   1.0 


-   0.8 


a15  .   1.0 


a   «b  0.05    (Initial 

Setting) 

a1    =   0.05 
(Initial   Setting) 


" 


As; 


fioc:  at     :    _  ?;.i.i 


''OS     an    j       ■  .         S    8  3    Of    £ 


i.        -      I 


R2   =   -    . 


-'.  '»      ^     ;• 


- 


R3    -   D.l    Megs      Rf2    -        .-. 

i    H&& 


12 


R  =  0.  1    Meg, 
4 

R  s  i  Mei 

j 

R6  =  I   Megs 

R?  -.  I   Megs 

^8  -  -    Mteg, 


ho  =  °'S  Megs   Cf3 

R,,    ^0.1    Meg,    C        =    I 

11  :  - 


i< 

.^ 

Jl 

ux 

c .„ 

-■ 

, 

- 

1 

. 

f3 

=2 

■-  i 

I 


£.4 


f 


K 


13 


0.1   Meg,    C 

C  2 


;,.   -   o.i       .. 


-      :! 


15 


-  i 


. 


For  Dead    ■  5  Voli 

.-  i    ■    of   Z e  - 

For    5    volts    Relay   Oucp^c 
Volt 
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^i^ngH 


M 

■ 

- 

j. 
•H 


: 


I 


! 
■ 


Q 


IF 
1*4 


I 


SO 


3 

i 

■  . 

- 

': 

! 


-M 


12:4 


Io  accurately  set  the  back-biasn, 
generatoi  was  connected  to  the  inp  t  of  1 

An  X-Y  recorder  was  then  calibrated  to  2  volts  p 
The  input  and  output  of  the  relay  simulator  were  th«  i  c         i   Lh 
respective  X  and  Y  axis  terminals.   The  pots  "a"  ...  - 

justed  until  the  X-Y  recorder  plot  closely  approximated  the  static 
characteristic  of  the  chosen  relay  with  dead  zone.  T\ 
istic  obtained  by  the  X-Y  recorder  is  shown  in  Figui   IV-A-4.2, 

The  system  was  first  excited  by  a  iO  yolt  step  input  signal  and  th 
system  was  found  to  build  up  into  a  limit  cycle  as  shown  in  r«  :ordi  i 
IV-A-4.1.   The  system  was  then  excited  by  a  20  volt  step  input  signal 
found  to  go  into  the  same  limit  cycle  as  shown  by  recorder  trace  iV-A-4.20 
From  these  recorder  traces  the  limit  cycle  was  found  to  have  an  angular 
frequency  of  0.455  rad/sec  and  an  amplitude  of  14.3  volts. 

The  phase  portrait  of  the  system  is  shown  in  Figure  IV-A-4.3. 
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RECORDER  MARK  H 


iTdet   l 


■      -A-^, 


Relay  with  Zzad   Zoiae  in'-Errc 

with  No  Zeros,  20  Volt  Step 

SCALES »  ■ 

0 

£  -  1  Volt/Lire       £,  -  0.5  Volt/] 


, 


'  s/Lii 


IV-A-5  Comparison  of  the  Four  Methods  :  ..  Limit  - 

a  Nonlinear  System,  of  Figure  IV-A. 


Method 

;.'..,.- 

rad/ £  a 

(  '  •  -  s ) 

Mitrovic" s 

[  ,46(  ■ 

10. 

Root  Locus 

0.4! 

, 

Describing  Function 

0.455 

14.1 

Analog  Simulation 

0.4 

1 4  .  3 

.30 


IV-B  Loop  Transfer  Function  with  One  Zero. 

The  block  diagram  of  the  system  chosen  is  shown  in  Figure  IV-B. 


0^ 


N 


x 

5\ 


-s 


s  e 


■  -fir 


G(s)=   aoo(S»i) 


s(.s+o.s"Xs+ 4x^+8) 


Figure  IV-B 
Relay  with  Dead  Zone  in  Error  Channel  Loop 
Transfer  Function  with  One  Zero 
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IV-B-1  Mitrovic's  Method. 

The  characteristic  equation  of  the  given  system  is 

sV  \2SSZ  +  Z*s\  (IL  +S0ON)  s  +  FOoN  -0. 

(IV-B-1. i) 
Assuming  the  last  two  coefficients  variable,  the  characteristic 

equation  becomess 

S4f  \2SsZ+?>%$1  +B,  5+8o-0  (iv-B-i. 2) 


where 


B  =  500  N  (IV-B-1. 3) 


and 


B  =  16  +  500  N  •  (IV-B-1. 4) 

From  Appendix  A,  Mitrovic's  equations  for  B  and  B  are 

Zo'-i^^4(t^)^nSu^L^)^Uj,%L'i)]  (IV-B-1. 5) 


and 


/3,  -  3?a>»^ffl)  +  iatd-tOii*^£iO^  tofrkOi).  (iv-B-i. 6) 

Substitution  of  the  values  of  the  0   functions  from  Appendix  A  for 
^  =  0  gives  the  parametric  equations  of  the  stability  curve  as, 

B  =  38  Ujfl2  +  0b^  (IV-B-1. 7) 

and 

B  -  12.5  60n  (IV-B-1. 8) 

At  the  intersection  of  the  M-point  locus  and  the  stability  curve, 
equation  IV-B-1. 3  equals  Equation  IV-B-1. 7  and  Equation  IV-B-1. 4  equals 
Equation  IV-B-1. 8,  thus 

500  N  =  38  Ud^2     +  60*4  (IV-3-1.9) 


and 


16  +   500  N  =  12.5  U)A2  (IV-B-1. 10) 


132 


Substituting  Equation  IV-B-1.9  into  Equation  IV-B-1.10  and  combining 
terms  yields 

&o/-2s£X£0*l-/&  =  G    ♦  (iv-B-i.ii) 


2 
Using  the  quadratic  formula  to  solve  Equation  IV-B-1.11  for  LO^ 


yields 


U)A2     =  26.11 


from  which  the  frequency  of  the  limit  cycle  is  determined  to  be 

60^  =  5.11  rad/sec. 

2 

Substituting  (jOn     back  into  Equation  IV-B-1.9  the  "average"  gain  of 

the  nonlinearity  over  a  cycle  is  determined  as 

N  =  0.6208 
from  which  the  approximate  amplitude  of  the  limit  cycle  is  determined 

to  be 

E     =  8.054  volts. 
max 

A  graphical  presentation  of  the  intersection  is  sketched  in  Figure 
IV-B-1.1. 

Mitrovic's  Method  predicts  a  stable  limit  cycle  by  applying  the 
same  arguments  given  in  Section  IV-A-1  to  Figure  IV-B-1.1. 
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A/*  \,0 


2Z7S  12C4 

Figure  IV-B-1.1 

Graphical  Solution  for  Relay  with  Dead  Zone  in  Error 

Channel,  Loop  Transfer  Function  with  One  Zero 


476-  6, 
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IV-B-2   Root  Locus  Method 

Substituting  S  =  Ju)    in  the  characteristic  Equation  IV-B-1.1 
and  evaluating  the  higher  powers  of  J  gives. 

Requiring  that  the  real  and  imaginary  parts  of  Equation  IV-B-2 . 1 

go  to  zero  independently  yields  the  following  equations, 

CO4  -38-60^  -h^O0A/=6  (IV-B-2, 2) 

and 

-Jli^^-t  ^(/^+r^0M)u)  -O  (IV-B-2. 3) 

Eliminating  N  between  equations  IV-B-2. 2  and  IV-B-2. 3  gives, 

Ld^-IFS^  -IL-&  (IV-B-2. 4) 

2 
Using  the  quadratic  formula  to  solve  Equation  IV-B-2. 4  for  Co 


gives. 


602  =  26.11 


from  which  the  frequency  of  the  limit  cycle  is  determined  to  be 

U)   =  5.11  rad/sec. 

2 

Substituting  60   back  into  Equation  IV-B-2. 2  gives,  the  "average" 

gain  of  the  nonlinearity  over  a  cycle  as 

N    =   0.6208 
ave 

from  which  the  approximate  amplitude  of  the  limit  cycle  is  determined 

to  be 

£      =  8.054  volts. 

max 

A  graphical  presentation  of  the  root  locus  solution  is  sketched 
in  Figure  IV-B-2. 1. 

The  root  locus  method  predicts  a  stable  limit  cycle  by  applying  the 
arguments  given  in  Section  IV-A-1  to  Figure  IV-B-2. 1. 
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LINEAR    SYSTEM    CL0$£O    LOOP 


-LIMIT  CYCLE 
S.UO^/ssc 


Figure   IV-B-2.1 

i 

Root   Locus   Solution   for   Relay  with  Dead  Zone   in 
Error  Channel  Loop  Transfer  Function  with  One   Zero 


+<T 
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IV-B-3  Describing  Function  Method. 

The  describing  function  for  a  relay  with  dead  zone  is  given  in 
Equation  IV-A-3.1,  The  linear  loop  transfer  function  for  the  given 
system  in  frequency  response  form  is 


GlM 


3L25"P(J0+|) 


(IV-B-3.1) 


JUJ(e3UHi)(O,25jaj  +  /J(o.|2o-juL>+0 
Table  IV-B-3.1  gives  the  magnitude  and  phase  of  Equation  IV-B-3.1  for 
the  various  frequencies  used  to  plot  the  linear  loop  transfer  function  on 
the  gain-phase  plane  of  Figure  IV-B-3.1.   These  values  are  obtained  from 
A  Bode  diagram  plot  of  Equation  IV-B-3.1. 

TABLE  IV-B-3.1 


CO 

\G(M)\  6b 

/6(xu)  ae9 

4.65 

6 

-176 

4.95 

5 

-178.5 

5.1 

4.6 

-180 

5.3 

4 

-181.5 

5.6 

3 

-184 

5.9 

2 

-187 

le  IV-B-3.2  gives  the  values  of  1/G  which  are  plotted  in  Figure 
IV-B-3.1.   These  values  are  determined  from  Equation  IV-A-3.1. 


TABLE 

IV-B-3.2 

O  W<X;( 

S 

1/GD 

1/G     db 

8 

0.621 

1.61 

4.14 

8.9 

0.592 

1.69 

4.56 

9 

0.588 

1.7 

4.62 

10 

0.551 

1.814 

5.175 
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Inspection  of  Figure  IV-B-3.1  shows  the  intersection  and  thus  the 

limit  cycle  to  occur  at ' 

CO    -   5.1  rad/sec. 

and 

f  -  8.95  Volts. 

v —  max 

The  describing  function  method  predicts  a  stable  limit  cycle  by 
applying  the  arguments  given  in  Section  IV-A-3  to  Figure  IV-B-3.1. 
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IV-B-4  Analog  Computer  Simulation. 

The  analog  computer  simulation  of  the  given  system  is  shown  in 
Figure  IV-B-4. 1.   The  coefficient  pot  settings  with  associated  resistances 
and  capacitances  for  real  time  and  real  magnitude  scaling  are  given  in 
Table  IV-B-4. 1. 

TABLE  IV-B-4. 1. 


Pot  Setting 

Associated  Elements 

E. 

For  an  input  step  of 

E 

i  volts 

J. 

ai  =  Too 

Rx  =  1  Meg,     Rfl  = 

1  Meg 

a2  =  1.0 

R2  =  1  Meg,     Rfl  - 

1  Meg 

a3  =  1.0 

R3  =  0.1  Meg,   R  2  - 

1  Meg 

a,  =  1.0 

R  =0.1  Meg,   Rf3  = 

1  Meg 

ac  -  0.3 

R5  =  0.1  Meg,   Cfl  = 

1  uf 

a6  =  0.4 

Rfi  =  0.1  Meg,   Cfl  = 

1  uf 

a?  =  1.0 

R?  =  0.1  Meg,   Rf^  = 

1  Meg 

a8  =  1.0 

Rg  =  0.1  Meg,   Rf4  = 

1  Meg 

a9  =  1.0 

R9  =  0.1  Meg,   Cf2  = 

1  uf 

a10=0.5 

R1Q  =  1  Meg,    Cf2  = 

1  uf 

ail^0"5 

R   =  0.5  Megs  Cf3  = 

1  uf 

a,,  =  0.3 
12 

R12  =  0.1  Meg,  Cf3  = 

i  uf 

au  =  1.0 

R13  =  1  Meg,    Cf4  = 

1  uf 

a  =  0,05  (Initial 

For  5  volt  dead  zone 

either 

Setting) 

side  of  Zero 

i 
a'  =  0,05  (Initial 

For  5  Volt  relay  output  Voltage 

Setting) 

To  accurately  set  the  back-biasing  pots  "a"  and  "a  ",  the  same  pro- 
cedure as  given  in  Section  IV-A-4  was  followed.   The  static  character- 
istic of  the  relay  is  the  same  as  that  shown  in  Figure  IV-A-4. 2. 
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The  system  was  lust  excited  with  a  6  volt  step  input  signal  and 
was  found  to  build  up  into  the  limit  cycle  shown  in  recorder  tra 
XV-B-4,,1.   The  ssytem  was  then  excited  with  a  12  volt  step  input  signal 
and  the  system  was  found  to  go  into  the  same  limit  cycle  as  shown  in 
recorder  trace  IV-E-4.2.   From  these  two  recorder  traces  the  limit  cycle 
was  determined  to  have  an  angular  frequency  of  5.07  rad/sec  and  an 
amplitude  of  9„0  volts. 

The  phase  portrait  of  the  system  is  shown  in  Figure  IV-A-4.2. 
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IV-B-5  Comparison  of  the  Results  of  the  Four  Method        -  System 

Of    X    r-  .     IV  ~  B  0 


Method 

r  i  ■ 
rad/st 

Amp]  itude 
s) 

Mi t     "  s 

Root  i . 

Describing  Fumct ion 

A: alog  .  . 

5.11 

Soil 

5.1 

5.07 

- 

B  ,  354 

8.! 

CHAPTER  V 
NONLINEARITY  IN  A  FEEDBACK  PATH 
V-A  Saturation  in  Velocity  Feedback  Path. 

The  system  chosen  to  analyze  is  given  in  Figure  V-A. 


<a 


-A 


-A 


y 


N 


+  3 


Y 


-a 


6  5 


■>- 


figure    t-a 
saturation    in   velo^/ty  feedback   path 

V-A-l  Mitrovic's  Method. 

Assigning  N  as  the  instantaneous  gain  of  the  nonlinearity,  the 
characteristic  equation  can  be  written  as, 

s4  +  7s3  +  14s2  +  (  8  +  120N)S  +  30  =  0.  (V-A-l. 1) 

Assuming  the  last  two  coefficients  variable,  the  characteristic 


equation  is  rewritten  as, 


where 


4     3      2 
s  +  7s  +  14s  +  B,S  +  B  =  0 

1     o 


B    =  8  +  120N 


and 


B   =   30  . 
o 

From  Appendix  A,  Mitrovic's  equations  for  B  and  B1  are, 


(V-A-l. 2) 


(V-A-l. 3) 


(V-A-l. 4) 
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Bo  =  -DW0(flS)  +  7<«>..3<kW  +u>*  tit*}"]  (V-A-1.5) 

and 

j3,=   I^(14«)+7W,^3(^)+^3^W).  (V-A-1.6) 


Substituting  the  values  of  the  0  functions  from  Appendix  A  for  <£   -   0 
gives  the  parametric  equations  of  the  stability  curve  as 


and 


Bq   =   14  OOflt  UJki  (V-A-1.7) 


B,    =      7\£A     .  (V-A-1.8) 


From  equations  V-A-1.3  and  V-A-1.4  the  M-point  locus  is  seen  to  be 
a  straight  line  segment  parallel  to  the  B.  axis.   At  the  intersections 
of  the  M-point  locus  with  the  stability  curve,  equation  V-A-1.4  equals 
V-A-1.7  and  equation  V-A-1.3  equals  equation  V-A-1.8.   Thus  at  the  inter- 
sections 

4  1. 

60^4    14^       =   36  (V-A-1.9) 

and 

2.  - 


760^=  &  +  IXON  ♦  (v-A-i.io) 


3 

Solving  equation  V-A-1.9    for      Uj^     gives 


and 


from  which 


CoS  2.6411 


U)fi  =  11.3589, 


and 


U\    ■   1.625   rad/sec 

U)^=   3.370   rad/sec. 
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Since  the  system  is  fourth  order,  the  stability  curve  is  a  para- 

bala  with  its  axis  parallel  to  the  B  axis.   Consequently,  both  values 

o 

of  (X)  f\      must  be  considered  since  the  M-point  locus,  a  straight  line 

parallel  to  the  B.  axis,  can  intersect  the  parabola  in  two  places. 

For    LOr\»  *  1«625  rad/sec,  the  average  gain  of  the  nonlinearity 

over  a  cycle  is  found  from  equation  V-A-1.10  to  be 

N,    =  0.0874, 
lave 

from  which  the  approximate  amplitude  of  signal  X  in  Figure  V-A  is  deter- 
mined to  be 

V .    =  34.325  volts. 
f\  lmax 

For  (jtL=  3.370  rad/sec,  the  average  gain  of  the  nonlinearity  over 

a  cycle  is  determined  to  be 

N_    =  0.5958, 
2ave 

from  which  the  average  amplitude  of  signal  X  is  determined  to  be 

X0    =  5.035  volts, 
/max 

However.,  the  signal  X  is  proportional  to  the  output  velocity  and 
not  the  output  angle.   To  obtain  an  approximate  amplitude  of  the  limit 
cycle  in  terms  of  output  angle,  O^  ,  one  must  assume  the  signal  X  is 
relatively  sinusoidal.   With  this  assumption,  the  amplitude  of  the  output 
angle  can  be  obtained  from 

9_    Am** 
6  "   fc.&OU  (V-A-l.ll) 

From  equation  V-A-l.ll,  the  approximate  amplitudes  of  the  output 

angle  for  both  limit  cycles  are 

A  =  3.52  volts 

^Oimax 

and 

0o2max=0'249  VoltS' 
149 


%  ■ 


The  graphical  presentation  of  the  solution  by  Mitrovic's  method 
is  shown  in  Figure  V-A-l.l. 
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Mitrovic's  method  predicts  that  the  limit  cycle  at  ^.=  1.625  rad/sec 
is  an  unstable  one  and  that  the  limit  cycle  at  (J^y\^~   3.37  rad/sec  is  a 
stable  limit  cycle  by  the  following  arguments.   For  the  stable  limit  cycle 
at  k-)fl2.=  3.37  the  same  arguments  used  in  Chapter  I,  Section  I-A-l  can 
be  applied  to  Figure  V-A-l.l.   For  the  unstable  limit  cycle  assume   Qji 
is  such  that  the  first  peak  of  signal  X  causes  the  instantaneous  M-point 
to  be  at  ML.   The  system  is  then  stable  and  the  oscillations  decrease.   As 
the  oscillations  decrease  N  increases  moving  the  M-point  to  the  right. 
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Eventually,  the  "average"  M-point  will  lie  at  the  intersection  of  the 
stable  limit  cycle  at  6\\i_-   3.37.   Now  assume  that  ^  is  large  enough 
to  cause  the  first  peak  of  signal  X  to  be  large  enough  to  place  the 
instantaneous  M-point  at  M  .   The  system  is  unstable  and  the  oscillations 
increase.   As  the  oscillations  increase,  N  decreases  causing  the  M-point 
to  move  even  farther  to  the  left  into  the  unstable  region.   Consequently 
the  oscillations  increase  indefinitely. 
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V-A-2   Root  Locus  Method. 

Substituting  5  -- J  lO   into  equation  V-A-l.l  and  evaluating  the 
higher  powers  of  ^  gives, 

604-J7603-  /4C02-  +j(8-H20N)60+30=  6.  (V-A-2.1) 

By  requiring  that  the  real  and  imaginary  parts  of  equation  V-A-2.1 
go  to  zero  independently,  then 

U>4  -  14  U)2-  irliO  =0  (v-a-2.2) 

and 

-0760^  +  ^fl+JOOAlW  -°  •  (V-A-2. 3) 

Dividing  equation  V-A-2.2  by  Jul)  the  above  two  equations  are 

exactly  the  same  as  equations  V-A-1.9  and  V-A-1.10  developed  in  Mitrovic's 

method.   The  solutions  of  these  simultaneous  equations  are  again  given 

£0  |  =  1.625  rad/sec 

(JO   2.=  3.370  rad/sec 

N.    =  0.0874 
lave 

N0    =  0.5958 
2ave 


X,    =  34.325  volts 
lmax 


and 


X0    =  5.035  volts. 

zmax 

Also  applying  equation  V-A-l.ll  gives 

Q0  .  =  3.52  volts 

w  lmax 

and 

gL  .,    =   0-249  volts. 
O  2max 

The  graphical  solution  of  the  root  locus  solution  is  given  in  Figure 
V-A-2.1.  In  the  figure,  the  "open  loop  poles"  are  the  poles  of  the  tran- 
sfer function  obtained  when  the  position  feedback  loop  is  closed  on  the 
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plant,  first  manipulating  the  1.5  gain  block  out  of  the  error  channel. 


l/near  System  closed  locp  root 


STABLE  LIMIT  CYCLE 
UNSTABLE  L/MIT  CYCLE 


ROOT  LOCUS    SOLI/TION     FOR 
SATURAT/ON    XN    VELOCITY 
FEEDBACK     PATH 


+<r 


The  root  locus  method  predicts  the  limit  cycle  at  60|  =  1.625  rad/sec 
will  be  an  unstable  limit  cycle  and  the  limit  cycle  at  k)^=   3.37  rad/sec 
will  be  a  stable  one.   The  arguments  of  Chapter  I,  Section  I-A-2  can  be 
applied  to  the   ^0^=  3.37  rad/sec  crossover  to  determine  a  stable  limit  . 
cycle  at  that  point.   At  the    60,  =  1.625  rad/sec  crossover,  assume  the 
input  signal,  Qj^,    ,  is  large  enough  to  cause  the  first  peak  of  signal 
X  to  be  large  enough  to  place  the  instantaneous  root  at  r,  .   The  system 
is  stable  and  the  oscillations  decrease.   As  the  oscillations  decrease, 
N  increases  causing  the  instantaneous  closed  loop  root  to  move  up  the 
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locus  tcward  the  linear  closed  loop  root. 

Eventually,,  the  average  root  location  over  a  cycle  will  lie  at  the 

stable  limit  cycle  crossover.   If  $/L  *-s  iarge  enough  to  cause  the  first 

peak  of  signal  X  to  be  large  enough  to  place  the  instantaneous  root  at 

r  ,  then  the  system  is  unstable  and  the  oscillations  increase.   As  the 
c 

oscillations  increase,,  N  decreases,  driving  the  closed  loop  root  farther 
into  the  right-half  plane  toward  the  open  loop  root.   Consequentlys  the 

oscillations  will  increase  indefinitely. 
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V-A-3  Describing  Function  Method. 

The  describing  function  for  saturation  is  given  in  Equation  I-A-3.1 
of  Chapter  I. 

To  obtain  an  equivalent  loop  transfer  function,  the  gain  of  the 
error  channel  was  first  made  unity  by  moving  the  gain  block  to  the  input 
channel  and  to  the  unity  feedback  path  as  shown  in  Figure  V-A-3. 1. 

Reduction  of  the  outer  position  feedback  thus  yields 


20 


<=>  0)    -  5«+  7s24  Msz-r  Ss+    SO     • 


(V-A-3. 1) 
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FIGURE  Y-A-3,./ 

A  digital  computer  program  was  used  to  find  the  roots  of  the  de- 
nominator of  equation  V-A-3. 1,  which  are, 
r.  =  0.2168  -  jl.3854 
r2  =  0.2168  +  j  1.3854 

r,  =  -  3.7168  -  j  1.2011 

3  J 

r,  =  -  3.7168  +  j  1.2011. 

4  J 

Thus  the  equivalent  open  loop  transfer  function  becomes, 


ZOS 


GM(s)  -    (s-0.2168  ±j  j.38^  )(s  f  3 tpca  ±Ji.ao/i) 


(V-A-3. 2) 
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With  S  =  Jul)  ,  the  magnitude  and  angle  of  equation  V-A-3.2  were 
calculated  using  a  digital  computer  program  with  fjd  varying  between  0.1 
rad/sec  and  100  rad/sec.   The  results  of  this  program  with  angles  in  the 
neighborhood  of  -180°  are  given  in  Table  V-A-3.1.   These  results  are  plot 
ted  on  the  gain-phase  of  Figure  V-A-3.2. 


TABLE  V-A-3.1 


CL> 

l&'HQutil  Jh 

/6'HM    ^ 

1.60 

21.65 

-182.57 

1.61 

21.46 

-181.51 

1.62 

21.27 

-180.50 

1.63 

21.08 

-179.54 

1.64 

20.89 

-178.64 

1.65 

20.71 

-177.79 

3.30 

4.903 

-179.02 

3.35 

4.605 

-179.72 

3.36 

4.556 

-179.86 

3.37 

4.498 

-179.995 

3.38 

4.441 

-180.134 

3.39 

4.384 

-180.27 

3.40 

4.327 

-180.41 

Values  of  the  1/G  locus  are  given  in  the  two  areas  of  interest  in 
Table  V-A-3.2.   These  values  are  calculated  from  Equation  I-A-3.1. 
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TABLE  V-A-3.2 


R 

Go 

1/6  o   J>b 

0.060 

0.07635 

22.34 

0.066 

0.08397 

21.52 

0.068 

0.08651 

21.26 

0.070 

0.08905 

21.01 

0.072 

0.09159 

20.53 

0.074 

0.09413 

20.53 

0.46 

0.5643 

4.969 

0.47 

0.5756 

4.798 

0.48 

0.5868 

4.630 

0.484 

0.5935 

4.532 

0.488 

0.5957 

4.499 

0.490 

0.5979 

4.467 

0.50 

0.6090 

4.308 

By  interpolation  of  Figure  V-A-3.2  and  the  tables,  two  limit  cycles 
are  seen  to  exist  at 

CO,    m   1.625  rad/sec 
/?,  =  0.0686 


and 


^irtef  43,7  volts 

UO^  ■   3-37  ^ad/sec 
Rx  =  0.488 

^imcf  6'15  volts      » 
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To  determine  the  approximate  amplitude  of  the  output  angle,  sinusoidal 

variations  of  the  signal  X  must  be  assumed.   Thus  equation  V-A-l.ll  can  be 

used  to  give,,  for  the  limit  cycle  at  U)\    =  1.625  rad/sec 

G0  ■,    ■  4.48  volts 
°  Imax 

and  for  the  limit  cycle  at  U)n~   3.37  rad/sec 

Q*   .    =  0.304  volts. 
wO  2max 

The  describing  function  method  predicts  a  stable  limit  cycle  at 
COx  ~   3.37  rad/sec  and  an  unstable  limit  cycle  at  60(  =  1.625  rad/sec. 
At   tOo-  3.37  the  limit  cycle  is  determined  to  be  stable  by  applying  the 
same  reasoning  given  in  Chapter  I,  Section  I-A-3  to  Figure  V-A-3.2  using 
the  signal  X.   At  d0|  =  1.625  the  limit  cycle  is  determined  to  be  unstable 
by  the  following  reasoning.   If  the  input,  @^     ,  is  large  enough  to 
cause  the  first  peak  of  signal  X  to  be  large  R  will  be  small  and  the  in- 
stantaneous critical  point  will  lie  at  say  R  in  Figure  V-A-3.2.   The 
system  is  unstable  due  to  a  negativephase  margin  and  the  oscillations  in- 
crease.  As  the  oscillations  increase,  R  becomes  even  smaller  driving  the 
critical  point  farther  up  the  -180°  axis  and  increasing  the  negative  phase 
margin.   Thus  the  system  is  always  unstable  and  the  oscillations  increase 
indefinitely.   If   £J\.   is  small  enough  such  that  the  first  peak  of  signal 
X  is  small  enoughs  the  instantaneous  critical  point  will  lie  at  say  R  of 
Figure  V-A-3,2.   The  system  is  stable  due  to  positive  phase  margin  and  the 
oscillations  decrease.   As  the  oscillations  decrease  R  becomes  smaller 
driving  the  critical  point  down  the  -180°  axis  and  increasing  the  positive 
phase  margin.   Eventually  the  positive  phase  margin  begins  to  decrease  due  to 
the  shape  of  the   G  H   O  COj      curve  and  the  system  reverts  to  the 
stable  limit  cycle  at   60o  =  3.37  rad/sec. 
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V-A-4  Analog  Computor  Simulation. 

The  analog  computor  simulation  of  the  given  system  is  shown  in 
Figure  V-A-4. 1.   The  coefficient  pot  settings  and  the  associated  resist- 
ances and  capacitances  for  real  time  and  magnitude  scaling  are  given  in 
Table  V-A-4. 1. 

TABLE  V-A-4. 1 


Pot  Setting 


Associated  Elements 


al  = 


0.75  Ei 
100 


a2  =  0.75 


1.0 


a.  =  1.0 
4 


ac  =  1.0 

5 


1.0 


a7  =  1.0 
a8  =  1.0 
a9  =  1.0 
a1()  =  0.4 

aH  =  1'° 
a12  =  °'6 

■13-i.o 

a^  =  1.0 

a  =  0.03 
Initial  Setting 


For  input  of  Ei  Volts  and  Ke  «  1.5 


R  =0.5  Meg, 
R2  =  0.5  Meg, 
R3  =  1  Meg, 
R4  =  1  Meg, 
R5  =  1  Meg, 
R6  =  1  Meg, 
R  =  0.5  Meg, 
Rg  =0.5  Meg, 
Rg  =  0.1  Meg, 
R1Q  =  0.1  Meg, 
Ru   =  1  Meg, 
R12  =0.1  Meg, 
R13  =  1  Meg, 
R14  =  1  Meg, 


Voltage 


Rfl  =  1  Meg 
Rfl  =  1  Meg 
Rf2  =  1  Meg 
Rf2  =  1  Meg 
Cfl  =  1  uf 
Cfl  =  1  uf 
Cf2  =  1  uf 
C  2  =  1  uf 
Cf3  =  1  uf 
Cf3  =  1  uf 
Rf3  =  1  Meg 
Rf4  =  1  Meg 


R£5  = 

1 

Meg 

Cf4  = 

1 

uf 

uration 

To  accurately  set  the  back-bias  pots  "a",  a  sine  wave  generator 
was  connected  to  the  input  of  the  nonlinearity  simulation  and  the  output 
monitored.   Pots  "a"  were  then  adjusted  until  the  output  was  limited  to 
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3  volts.   The  static  characteristic  curve  is  the  same  as  that  in  Figure 
I-A.4v2j,  except  the  saturation  voltage  is  3  volts  instead  of  5  volts. 

A  step  input  signal  of  1  volt  was  first  impressed  on  the  system  and 
the  system  was  found  to  build  up  into  the  limit  cycle  shown  in  recorder 
trace  V-A-4.1.   A  step  input  signal  of  5.8  volts  was  then  impressed  and 
the  system  was  found  to  go  into  the  same  limit  cycle  as  shown  in  recorder 
trace  V-A-4.2.   As  expanded  portion  of  the  stable  limit  cycle  is  shown  in 
recorder  trace  V-A-4.3. 

A  visual  presentation  of  the  unstable  limit  cycle  was  difficult  to 
obtain.   For  a  step  input  of  6.21  volts  the  system  was  only  slightly 
damped  before  reverting  to  the  stable  limit  cycle  as  shown  in  recorder 
trace  V-A-4.4.   For  step  input  of  6.22  volts  the  system  was  slightly  un- 
stable with  the  oscillations  eventually  increasing  indefinitely  as  shown  in 
recorder  trace  V-A-4.5. 

From  the  recorder  traces,  the  stable  limit  cycle  was  determined  to  bes 

OD o  =  3.35  rad  sec 

V  ,    =6.2  volts 
/x  lmax 

/--n  -    =  0.30  volts. 
Co  lmax 

The  unstable  limit  cycle  was  determined  by  averaging  recorder  traces 
V-A-4.4  and  V-A-4.5  with  measurements  taken  on  the  first  few  cycles  only^ 

COf     =  1.618  rad/sec 
Xi**x=  42«3  v°lts 
Qo\  mx}=   4.29  volts. 
The  phase  portrait  of  the  system  is  shown  in  Figure  V-A-4.2. 
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B  II US  II      LVSTR  IJINI  YZ  MTS  division  of  clevite  corporation 


Recorder  Trace  V-A-f . 1 
Saturation  in  Velocity  Feedback  Path,  1  Volt  Step  Input. 
SCALES;  £-0.1  Volts/Line    £  -   0.1  Volts/Line    X  -  0.5  Volts/Line 
Y  -  0.2  Volts/Line     Time  -  1  Second/Division 
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-  CLEVITE  CORPORATION  CLEVELAND.  OHIO  PRINTED  IN  USA 


i  ' 


Recorder  Trace  V-A-4,2 
Saturation  in  Velocity  Feedback  Path,  5.8  Volt  Step  Input. 
SCALES:£- 0.5  Volts/Line   £.  -0.5  Volts/Line     X  -  5  Volts/Line 
Y  -  0.2  Volts/Line     Time  -  1  Second/Division 
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Recorder  Trace  V-A-4.3 
Saturation  in  Velocity  Feedback  Path,  Limit  Cycle. 

SCALES:  £  -0.02  Volts/Line    £  -  0.1  Volts/Line   '  X  -  0.5  Volts/Line 
Y  -  0.2  Volts/Line   Time  -  0.2  Seconds/Division 
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Recorder  Trace  V-A-4.4 
Saturation  in  Velocity  Feedback#  Path,  6.21  Volts  Step  Input 

£  ;  °'n  Tl'li^      fi  "  °°5  Volts/Lin£      X  -  5 
Volts/Line     Time  -  1  Second/Division 
— JL66 


Recorder  Trace  V-A-4.5 
Saturation  in   Velocity  Feedback  Path,  6.22  Volts  Step  Input, 
SCALES-  £"- -0.5  Volts/Line      £-  0.5  Volt s/Line     X  -  5  Volts/l 

Y  -  0.2  Volts/Line        Time  -  1  Second/Division 
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V-A-5   Comparison  of  the  Four  Methods  in  Predicting  Limit  Cycles  of  the 
:en  -  -  Figun  ii 


Method 



Angular  Freiq  jiei     ■- 

Amplitude    '  ■>■:  Its) 

rad/sec 

X  max 

n 

E  max 

Stable.. 

.  .  ITnfltahle 

...„.  Stable 

Unstable 

Stable 

Unstable 

Mit rov:  ."s 

Root   Locus 

Describing 
Functic  c 

Analog 

Simulation 

3,37 
3o37 
3.37 

3.35 

525 
1,625 
1. 

1.618 

5.03 

5 .  03 

,15 

6.2 

34.33 
34.33 
43.7 

42.3 

0.249 
0.249 

0.304 

0.30 

3.52 

3.52 
4.48 

4.29 

V-B  Saturation  in  Acceleration  Feedback  Path. 

The  system  chosen  to  analyze  is  shown  in  Figure  V-B. 


7o 


Figure  V-B 
Saturation  in  Acceleration  Feedback  Path 


V-B-l.   Mitrovic's  Method. 

Assigning  N  as  the  gain  of  the  non-linearity  the  characterise 


equation  is: 

4   .,  3  2 

s  +  7s  +  (14  +  10N)s  +  8s  +  20  =  0  4 

1      2 
Assuming  the  s  and  s   coefficients  to  be  variable  the 

equation  becomes 


(V-B-l. 1) 
characteristic 


where 


and 


4,3      2 
s  +  7s*  +  B2s  +  B.s  +  20  =  0 


B  =  8 


I 


B9  =  14  +  10  N   „ 


(V-B-l. 2) 


(V-B-l. 3) 


(V-B-l. 4) 
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From  Appendix  A>,  MItrovic's  equations  for  B-  and  B«  are 


and 


Substitution      :  I  ions  frou   apt-  ndix   A  fos       ^     =0  gives 

the  parametric  equatic  •     ■     ibi  arve  as 


1 

and 


*  ^OOki  (V-B-1.7) 


1  =  "ST!    +  <On    .  Cr-B-1.8) 

The  shape  of  the  stability     -  on  the  B-  versus  B«,  plane  is  not 
easily  recognized.   H        it  the  J     secti  n  of  the  stability  curve 
and  the  Mhpoint  locusa         V°B=1«,  mist    squal  equation  V=B=1.3  or 

7<Hnx   -  8 

from  which  the  frequency   £  I     nit  cycle  is  determined  to  be 

60^  *  1.069  rad/sec. 

Also  at  the  intersection  equa  ion  V-B~1.8  must  equal  equation  V=B-1.4 
or 

14  +  ION  =  »v  -J-  00*  n    m 

Substitution  oi    CO    yields  the        |  ain  of  the  nonlinearity  over 
;•■  .    - 

N    =  0.4614 
from  which  the  approx:  signal  X  is  determined  to  be 


X    =  10.82  volts 
max 

To  obtain  an  approximate  amplitude  of  the  angular  oscillations,  one 
must  assume  X  to  be  sinusoidal.   Then  the  amplitude  of  the  error  signal 
may  be  determined  by 


&max 


or 


0,S    U)f 


(V-B-1.9) 


$,         =  18.94  volts. 
The  graphical  solution  by  Mitrovic's  method  is  sketched  in  Figure 
V-B-l.l. 
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Mitrovic's  method  predicts  an  unstable  limit  cycle  by  the  following 
reasoning,,   If  the  input  signal  is  large  enough  such  that  the  first 
peak  of  signal  X  is  large  enough 9  the  instantaneous  nonlinear  gain  will  be 
small  enough  to  place  the  M-point  at  M  in  Figure  V-B-l.l  the  system  is 

3 

unstable  and   the  oscillations  increase.  As  the  oscillations  increase  N  de- 
creases moving  the  M-point  further  into  the  unstable  region.   Thus  the 
average  M-point  can  never  be  in  the  stable  region  and  the  oscillations 
increase  indefinitely.   If  the  input  signal  is  small  enough  such  that  the 
first  peak  of  signal  X  is  small  enough,  N  will  be  relatively  large  placing 
the  instantaneous  M-point  at  M,  in  figure  V-B-l.l.   The  system  is  stable 
and  the  oscillations  decrease.   As  the  oscillations  decrease  N  increases 
approaching  one  and  moving  the  average  M-point  farther  into  the  stable 
region.   Thus  the  system  is  always  stable  and  all  oscillations  eventually 
die  out. 
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V-B-2.   Root  Locus  Method. 

Substitution  of  S  -  J(JO   into  the  characteristic  equation  and 
evaluation  of  the  higher  powers  of  J   yields 

C04-^76u3-04  +  ION)Lo2  +  ^  8co-r  20  -O.  (v-b-1.1) 

Requiring   that   both  the   real  and   imaginary  parts   of  equation  V-B-1.1 
go  to   zero   independently 

Oj4  _   (\A  -H0A/)cjo2  4  20    =-0  ^V-B-1.2) 

and 

~M  UJ3  4  SlU  =  0,  (V-B-1.3) 

By  dividing  equation  V-B-1.3  by    JLO  ,    the    frequency  of  the   limit   cycle 
is   obtained  directly  as 

uozu   (fyy    =    I.0L9      ntd/SeCi 

Substituting  (j£    into  equation  V-B-1.2,  the  average  gain  of  the  non- 
linearity  is  obtained  as 

N    =  0.4614 
ave 

from  which  the  approximate  amplitude  of  signal  X  is  obtained  as 

X    =10.82  volts, 
max 

To  get  the  approximate  amplitude  of  the  error  signal,  X  is  assumed 
sinusoidal  and  equation  V-B-1.9  is  used, 

£    =  18.94  volts. 

max 

The  root  locus  solution  is  sketched  in  Figure  V-B-2. 1.   Before 
sketching  the  figure,  the  unity  feedback  loop  was  first  reduced  and  the 
poles  of  the  reduced  transfer  function  were  used  as  the  poles  of  the 
equivalent  loop  transfer  function. 

174 


(-3.587  +J0.K5) 


-<r 


^w 


(-3.^87  -JOflW) 


L/N£flfl  StfTtM  (CLOSED  LOOP  ffCOT 


+JU> 


(0.  OSGrt + J 1. 2  OJ) 
'LIMIT  CVO£ 

tO-I.O^I-ft^ec 


FIGURE"  F-B-2.J 

ROOT  LOCUS  SQLUTIOhl   OF 

SATURATION  XW   ACCELERATION 

FEEDBACK  PATH 


+  <r 


*  (ooBflhtyaiO 


-iiJ 


The  root  locus  predicts  an  unstable  limit  cycle  by  applying  the  same  , 
reasoning  given  in  Section  V-A-2  to  Figure  V-B-2.1. 
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V-B-3.   Describing  Function  Method. 

The  describing  function  for  saturation  is  given  in  equation  I-A-3.1 
of  Chapter  I. 

To  obtain  an  equivalent  loop  transfer  function  the  unity  feedback 
loop  was  first  reduced  which  gives, 


&'(s)  s  s«^s»+?4°ssas-no        .  <V-B"3-1) 


A  digital  computor  program  was  used  to  find  the  roots  of  the  de- 
nominator of  G°(s). 

r.  =  0.0869  +  j  1.201 
1  J 

r2  =  0.0869  -  j  1.201 
r3  =  -3.5869  +  j  0.9652 
r4  =  -  3.5689  -  j  0.9652 
Thus  the  equivalent  open  loop  transfer  function  becomes 

r'uf  \    -  '^^ 

03  H(5J  "  (5-0.08^iJ|<20l)(Sr-3.5'8^±J0.^5aJ   , 

(V-B-3.2) 
With   5=J^   i  the  magnitude  and  angle  of  equation  V-B-3.2  were 
calculated  by  a  digital  computor  program  with  OO  varying  between  0.1 
rad/sec  and  100  rad/sec.   The  results  of  this  program  wit;  angles  in  the 
neighborhood  of  -180°  are  given  in  Table  V-B-3. 1.   These  results  are  also 
plotted  on  the  gain-phase  plane  of  Figure  V-B-3. 1  as  the  linear  loop  trans- 
fer function. 
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TABLE  V-B-3.1 


» 

60 

fG'tfooal  <u 

/G'AKM  ^o 

1,05 

5.599 

-182.94 

1.06 

6.153 

-181.48 

1.07 

6.721 

-179.83 

1.08 

7.301 

-177.97 

3.2. 


The  values  of  1/GD  in  the  area  of  interest  are  given  in  Table  V-B- 


TABLE  V-B- 3. 2 


R 

GD 

1/Gn  db 

0.30 

0.3762 

8.493 

0.35 

0.4364 

7.203 

0,37 

0.4601 

6.743 

0.372 

0.4625 

6.698 

0.374 

0.4648 

6.654 

0.38 

0.4719 

6.523 

0.40 

0.4954 

6.101 

By  interpolation  of  Figure  V-B-3.1  and  the  tables,  the  limit  cycle 


xs  at 


and 


from  which 


6l>  =  1.609  rad/sec 


R  =  0.374 


X    =  13.37  volts, 
max 
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To  obtain  an  approximate  amplitude  of  the  error  signal,  X  is 

assumed  sinusoidal  and  equation  V-B-1.9  used, 

C    =  23.4  volts, 
max 

The  describing  function  method  predicts  an  unstable  limit  cycle  by 
applying  the  same  reasoning  given  in  Section  V-A-3  to  Figure  V-B-3.1. 
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V-B-4.   Analog  Computer  Simulation. 

A  good  simulation  of  the  acceleration  signal  was  difficult  to  obtain. 
The  simulation  that  was  finally  decided  upon  was  to  use  the  velocity  signal 
and  approximately  differentiate  it  by  simulating  s/(s  +  100).   By  using 
this  transfer  function  an  extra  gain  of  100  had  to  be  simulated  in  the 
feedback  path  in  order  to  maintain  the  D.C.  gain  constant  of  the  path  un- 
altered.  Thus  the  simulated  gain  in  the  feedback  path  is  100  K  or  50 
instead  of  0.5. 

The  analog  computor  simulation  is  shown  in  Figure  V-B-4. 1.   The  co- 
efficient pot  settings  and  associated  resistances  and  capacitances  for 
real  time  and  magnitude  scaling  are  given  in  Table  V-B-4. 1. 

The  diode  back-biasing  pots  were  set  in  the  same  way  described  in 
Chapter  I.   The  static  characteristic  curve  is  the  same  as  Figure  I-A-4.2. 

An  input  signal  of  17  volts  was  first  applied  and  the  system  was 
found  to  be  stable  but  only  slightly  damped  as  shown  in  recorder  trace 
V-B-4. 1.   An  input  signal  of  23  volts  was  then  applied  and  the  system  was 
found  to  be  unstable  as  shown  in  recorder  trace  V-B-4. 2.   The  unstable 
limit  cycle  was  found  to  occur  for  an  input  step  of  19.8  volts  and  is 
shown  in  recorder  trace  V-B-4. 3. 

From  recorder  trace  V-B-4. 3,  the  limit  cycle  is  determined  to  have  an 

amplitude  of  X    =  10.6  volts  and  a  frequency  of   LQ    =   1.03  rad/sec. 
max 

The  phase  portrait  of  the  system  is  shown  in  Figure  V-B-4. 2. 

The  simulation  of  the  acceleration  signal  is  admittedly  poor,  especial- 
ly since  the  amplitude  of  X  differs  considerably  from  that  predicted  by 
the  describing  function. 
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Table  V-B-4.1 

Pot  Setting 

Associated  Elements 

Ei 
ai  =  100 

For  a  step  input 
R:  =  1  Meg, 

of  Ei  Volts 
Rfl  -  1  Meg 

a2  =  1.0 

R2  =  1  Meg, 

Rfl  =  1  Meg 

a3  =  1.0 

R3  =  1  Meg, 

Rf2  =  1  Meg 

a4  =  1.0 

R4  =  1  Meg, 

Rf2  =  1  Meg 

a5  =  1.0 

R5  =  1  Meg, 

Cfl  =  1  uf 

a.    ■  1.0 
6 

R6  =  1  Meg, 

Cfl  =  1  uf 

a  -  1.0 

R_  =  0.5  Meg, 

Cf2  =  1  uf 

ag  =  1.0 

Rg  =0.5  Meg, 

Cf2  =  1  uf 

a9  =  1.0 

R9  =  0.1  Meg, 

Cr„  =  1  uf 
f3 

aio  =  °'4 

R1Q  =  0.1  Meg, 

Cf3  =  1  uf 

an  =  1.0 

Rn  -  0.1  Meg, 

C,.  =  0.1  uf 
f4 

a12  =  U0 

R12  =  0.1  Megs 

Cc,    =   0.1  uf 
f4 

a^  =  1.0 

R13  =  0.1  Meg, 

Rf3  =  1  Meg 

a^  -  1.0 

R14  =  0.1  Meg, 

Rf3  =  1  Meg 

•15  =  0.5 

R15  =  0.1  Meg, 

Rf4  =  1  Meg 

a.,  -  1.0 
16 

R16  =  1  Meg, 

Rf5  =  1  Meg 

a„  =  1.0 

R1?  =  1  Meg, 

Cf5  =  1  uf 

alg  =  1.0 

R18  =  1  Meg, 

Rf6  =  1  Meg 

a  =  0.05  Initial  Setting 

For  a  saturation 

of  5  volts 

Although  the  nonlinearity  input  X  is  not  strictly  sinusoidal,  it  is 
not  too  far  from  a  sinusoid.   The  large  discrepancy  can  not  be  accounted 

for  solely  by  wave  shape.   However,  if  a  root  locus  solution  similar  to 

2 
that  of  Section  V-B-2  is  done  on  the  system  assuming  50s  /(s  +  100)  in 

the  inner  feedback  path  the  following  results  are  obtained: 
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CO    =  1.0725 


X    -  13.6 

max 


When  the  results  are  compared  with  those  of  the  root  locus  solution, 
Section  V-B-2,  it  is  seen  that  the  extra  pole  in  the  feedback  path  causes 
a  slight  increase  in  frequency  and  a  significant  increase  in  the  approxi- 
mate amplitude.  Consequently,  whatever  it  was  in  the  actual  analog  simu- 
lation that  caused  the  frequency  to  be  less  (1.03  rad/sec)  than  the  pre- 
dicted value  (1.069  rad/sec)  also  caused  the  amplitude  to  be  considerably 
less  than  predicted  by  the  describing  function. 
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SCALES 


£.    2   volts/line 
y  -  0.5  volts/line 


S      -  2  volts/line 
X  -  1  volt/line 


Time  -  1  second/division 


RECORDER  MARK  II         chart  no.  ra 2921  32         brush    instruments 


Recorder  Trace  V-B-4.1 
Saturation  in  Acceleration  Feedback  Path,  17  volt  Step  Input 
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RECORDER  MARK  II         chart  no.  RA2921  32  brush    instruments 


SCALES 


£  -  2  volts/line 

Y  -  0.5  volts/line 


£  -  2  volts/line 
X  -  1  volt/line 


Time  -  1  second/division 


RECORDER  MARK  II         chart  no  RA2921  32  brus 


Recorder  Trace  V-B-4,2 
Saturation  in  Acceleration  Feedback  Path,  23  volt  Step  Input 


185 


SCALES 


£ 

Y 


2  volts/line 
0o5  -volts/line 


<£  -  2  volts/line 
X  -  1  volt/ine 


Time  -   1  second/division 


Recorder  Trace  V-B-4.3 
Limit  Cycle,  Saturation  in  Acceleration  Feedback  Paths  19,8  Volts  Step  Input 
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V-B-5.   Comparison  of  the  Four  Methods  in  Predicting  Limit  Cycles  of  the 

System. of  Figure  V-B. 


Method 

Frequency 

Amplitude   (Volts) 

rad/sec 

X-Max 

£,  Max 

Mitrovic's 

1.069 

10.82 

18.94 

Root  Locus 

1.069 

10.82 

18.94 

Describing  Function 

1.069 

13.37 

23.4 

Analog  Simulation 

1.03 

10.6 

20.0 
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PART  II 
A  STUDY  OF  LIMIT  CYCLES  OF  A  SYSTEM  WITH 
TWO  GAIN  VARIABLE  NONLINEARITIES  USING  MITROVIC'S  METHOD 

Introduction 
The  author  had  no  experience  in  the  problem  of  a  system  with  two 
non-linearities  and  consequently  had  no  idea  of  what  performance  to  expect. 
Thus  a  system  was  arbitrarily  chosen  and  an  analog  computer  study  of  the 
system  was  started. 

The  system  chosen  is  shown  in  Figure  VI.   This  particular  configura- 
tion was  chosen  primarily  because  each  of  the  nonlinear  gains  appears  alone 
in  two  of  the  coefficients  of  the  characteristic  equation,  which  is, 


4     3      2 
s  +  7s  +  14s  +  (8  +  KK  N0)s  +  KK  N.  -  0 

t  2  el 
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Figure  VI 
System  with  Two  Saturation  Nonlinear it ies 
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Assuming  the  last  two  coefficients  variable,  the  characteristic  equa- 


tion becomes 


where 


and 


§4  +  7  s3  +   14S2  +  b,8  +  B  =0  (Vl-b) 

1     o 


B  =  KK  N.  (VI-c) 

o     el 


B-  -  8  +  KK  N2  (VI- d) 

The  equations  of  the  stability  curve  are  the  same  as  those  in  Chapter 
Vs  Section  V-A-l  and  are  repeated  below 

B  =  14  OJrx   -  U)n  (Vl-e) 


and 


=  7CJUn2  (VI- f) 


This  particular  arrangement  easily  lends  itself  to  presentation  on  the 

B  versus  B„  plane  in  that 
o        1 

lo  The  stability  curve  is  independent  of  the  gain  K  since  the 

4   3       2 
stability  curve  only  depends  on  the  coefficients  of  the  s  ,  s  ,  and  s   terms, 

2„   Once  K  is  chosen  a  second  rectangular  grid  of  constant  K  and 

t 

K  can  be  drawn, 

e 

3.   Once  K  and  K  are  chosen  then  any  instantaneous  or  average 
t      e 

M-point  is  easily  plotted  using  the  K  -  K  grid.   For  example^  if  K  =6 

and  K  =4  and  by  some  means  an  M-point  is  determined  by  N„  -  0.5  and  N„  = 
t  12 

0.259  then  this  M-point  is  easily  plotted  at  the  intersection  of  the  K  -  3 

line  and  the  K  -  1  line  of  the  secondary  grid. 

t 

In  the  studys  three  terms  are  used  to  describe  certain  parts  of  the 
system  shown  in  Figure  VI  and  need  to  be  defined.   These  are: 

Plant  -  that  portion  of  the  system  represented  by  the  linear  trans- 


fer function        G  {^\  —    H 


SCS+JXS  +  ZXS  +  4) 
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ar  System  -  the  equivalent  linear  closed  loop  transfer 

function  with  both  feedback  loops  reduced  and  no  nonlinearities  present, 
or  the  linear  system  described  by  the  characteristic  equation  Vl-a  with 

Nj  and  N2  equal  to  1.0, 

Linear  M-point  -  that  point  on  the  B  versus  B,  plane  which  is 

o        x 

determined  by  equations  VI-c  and  Vl-d  with  N  and  N2  equal  to  1.0. 
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Chapter  VI 

ANALOG  COMPUTER  STUDY  AND  STATEMENT  OF  RESULTS 
Besides  the  poles  of  the  plant,  there  are  six  parameters  which  directly 
or  indirectly  affect  the  two  assumed  variable  coefficients  of  the  character- 
istic equation.   These  are  Q^  s  K  ,  E  .,  K,  K  ,  and  E  9.   As  a  starting  point 
the  plant  gaim9  K9  was  made  20  and  both  saturation  voltages  were  made  equal 
to  3  volts .   For  these  conditions  the  analog  simulation  is  shown  in  Figure 
VI- 1.  The  settings  of  the  coefficient  pots  and  their  associated  resistances 
and  capacitances  are  given  in  Table  VI- 1. 

The  step  input  signals  E„  or  ©^  ,  was  obtained  by  a  pot  reduction  of  a 
100  volt  d.c.  source  through  an  isolating  amplifier.  To  change  K  various 
combinations  of  a  8  R  ,  a-s  and  R„  were  used9  howevers  both  channels  had  to 
be  the  same.   To  change  K  various  combinations  of  a   and  R   were  used. 

Because  the  diode  resistance  is  not  zero  during  conduction  and  pot 
loading  effects  of  pots  "a"  and  "a"',  the  back  bias  pots  were  readjusted 
slightly  as  the  amplitude  of  the  limit  cycle  changed.   The  readjustments 
were  made  by  eye  such  that  on  the  average  the  output  signals  of  both  satura- 
tion simulators  were  3  volts.   A  typical  output  signal  of  a  saturation  simu- 
lator for  a  large  input  signal  is  sketched  in  Figure  VX-2. 

First  the  stability  curve  was  drawn  on  the  B  versus  B,  plane  and  then 

o         1 

a  second  grid  of  constant  K  and  K  lines  was  drawn  as  shown,,  in  Figure  VI-3. 

From  this  figure «,  it  is  apparent  that  one  should  get  a  good  indication  of  the 

limit  cycles  possible  by  taking  many  linear  starting  points  located  by  K  and 

K  and  observing  the  resulting  limit  cycles. 

The  initial  study  took  exactly  this  form,  that  is,  a  K  was  chosen  and 

several  responses  to  a  step  input  were  observed  for  various  values  of  K  on 

t 
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CABLE  VI- 1 

Pot  Setting 

Associated  Elements 

a,  =  K  for  K  4  1.0 
1    e      e 

Rx  =  1  Meg5     Rfl  =  1  Meg 

a,  -  fl  for  K  >  1.0 
1   10       e  ' 

Rx  =0.1  Meg,   Rfl  =  1  Meg 

a„  =  K  for  K   4  1.0 

2    e      e  ^ 

R2  -  1  Megs    Rfl  -  1  Meg 

a„  =  ~  for  K  S   1.0 
2   10      e  S 

R2  -  0.1  Meg,   Rfl  »  1  Meg 

a  ■  1.0 

R3  =  1  Meg,    Rf2  -  1  Meg 

a.  ■  1.0 
4 

R4  -  1  Megs    Rf3  -  1  Meg 

a  ■  1.0 

R5  -  1  Megs    Rf3  =  1  Meg 

a,  ■  1.0 
6 

R6  =  1  Meg,    Cfl  -  1  uf 

a_  -  1.0 

7 

R?  -  1  Meg,    Cfl  =  1  uf 

ag  =  1.0 

Rg  =0.5  Meg,   Cf2  =  1  uf 

«,  -  1-0 

Rg  -  0.5  Meg,   Cf2  -  1  uf 

aK>  '  1-° 

R1Q  =  0.1  Meg,  Cf3  =  1  uf 

"11  =  °'4 

R11  -  0.1  Megs  Cf3  -  1  uf 

au  .  l.o 

R12  -  1  Meg,   Rf4  =  1  Meg 

a13  "  Kt  for  Kt  <■  U0 

R13  *  1  Meg    Rf5  =  1  Meg 

313  "  I  f0r  Kt  >   1-° 

R13  =0.1  Meg,  Rf5  -  i  Meg 

a.,  m   1.0 
14 

R .,    =  1  Meg,    Rf6  =  1  Meg 

a,c  =  1.0 

13 

R15  =  1  Meg,    Cf4  -  1  uf 

a16  .  1.0 

R16  -  1  Meg,   Rf?  =  1  Meg 

a  =  0.03  Initial  Setting 

For  a  3  volt  saturation  level 

a"  =  0.03  Initial  Setting 

For  a  3  volt  saturation  level 
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a  horizontal  traverse  across  the  plane ,   In  the  study,  K  was  varied  from  K 

r  J  J      e  e 

^  0,5  to  K  =  6.0.   K  was  varied  from  I  ■  0.5  to  K  =10.0.  About  every 
e         t  t  t  7 

other  value  of  K  9  the  input  step9  E.  y  was  varied  over  several  values  to  see 
if  the  resulting  limit  cycle  was  dependent  on  the  magnitude  of  the  step  in- 
put signal.   The  only  limitation  placed  on  E.  was  overloading  of  any  one  of 
the  computer  amplifiers;  consequent ly^  for  certain  linear  starting  points 

Es    was  as  large  as  200  volfs.   The  results  of  these  computer  runs  are 
imax 

tabulated  in  Table  I  of  Appendix  B» 

From  this  first  portion  of  the  analog  computer  study  several  things 
were  evident ,   These  were: 

1.   No  limit  cycles  were  obtained  from  a  step  input  signal  when 
the  linear  system  M-point  was  inside  the  stable  region  on  the  B  versus  !„ 
plane  no  matter  what  the  input  signal  magnitude,, 

2a  Limit  cycles  were  only  obtained  for  step  input  signals  when  the 

linear  system  M=point  was  outside  the  stable  region  and  these  limit  cycles 

were  independent  of  input  signal  magnitude  and  only  dependent  on  the  linear 

system  M-point  position,, 

3o  The  unstable  region  of  the  B  versus  B„  plane  can  be  divided  in- 

°  o        1 

to  three  regions  such  that  only  one  of  the  nonlinear  elements  saturate  or  both 
elements  saturate  for  the  linear  system  M- point  lying  in  a  given  region. 

Several  runs  were  made  for  various  linear  M-points  inside  the  stable 
region  and  for  various  step  input  signal  magnitudes.   From  the  results  of 
Section  V-A  of  Chapter  v\,  one  might  expect  for  linear  M-points  close  to  the 
left  portion  of  the  stability  curve  that  the  system  could  be  driven  unstable 
for  large  step  input  signals.   However^  all  linear  M-points  chosen  showed 
some  damping  even  for  the  linear  point  extremely  close  to  the  stability 
curve,  A  typical  response  for  a  linear  M-point  very  close  to  the  left 


portion  of  the  stability  curve  is  shown  in  Recorder  Trace  VI -1  for  a 
small  step  signal  input  and  Recorder  Trace  VT-2  for  a  large  st  sp  signal 

input.   Typical  responses  for  a  linear  M- point  very  close  to  the  right 
side  of  the  stability  curve  are  shown  in  Recorder  Trace  VI-4  for  a  large 
step  signal  input.  Typical  responses  ::"  i    i  linear  M-point  in  the  well 
damped  portion  of  the  stability  region  are  shewn  in  Recorder  Trace  VI- 5 
for  a  small  step  signal  input  and  Recorder  Trace  Vl-6  for  a  large  step  signal 
input. 

Although  it  has  been  stated  that  no  limit  cycles  were  observed  for  a 
linear  M=point  In  the  stable  region9  the  system  could  be  driven  into  a  limit 
cycle  for  repeated  disturbances  in  the  input  signal  after  the  initial  step 
input.   This  type  of  response  was  only  noted  for  linear  M-points  with  a  very 
small  damping  ratio  and  close  to  the  left  portion  near  the  maximum  of  the 
stability  curve.,  A  typical  response  of  this  type  is  shown  in  Recorder  Trace 
VI-7  for  extremely  little  damping.   In  this  trace,  a  limit  cycle  was  achie- 
ved with  only  one  disturbance  after  the  initial  step.   Recorder  Trace  VI-8 
shows  the  same  type  of  response  but  with  more  damping.,   In  this  trace  three 
disturbances  were  required  to  achieve  a  limit  cycle. 

To  obtain  these  responses  an  initial  step  input  signal  was  applied. 
Then  with  the  system  in  a  very  oscillatory  slightly  damped  condition  the 
disturbance  was  applied.   To  apply  the  disturbance.,  the  input  pet  was  given 
a  small  quick  twist  to  approximate  a  step  change  in  input  signal.   It  . 
found  that  the  disturbances  had  to  be  (1)  generally  small  in  amplitude  such 
that  the  error  channel  did  not  appreciably  saturate*,  (2)  applied  at  a  defin- 
ite time  during  a  cycle  and  in  a  specific  direction.   It  was  found  that  the 
maximum  effect  was  achieved  when  a  positive  disturbance  was  applied  when  the 
error  was  passing  through  zero  from  minus  to  plus.   Alsos  the  maximum  eff- 


was  observed  for  a  negative  disturbance  as  the  error  was  passing  through 
zero  from  plus  to  minus.   This  type  of  system  performance  was  only  found 
late  in  the  experimental  phase  of  the  study  and  was  not  investigated 
thoroughly. 

Recorder  Traces  VI-9  and  VI-10  show  that  the  resulting  limit  cycle  for 
a  linear  M-point  in  the  unstable  region  is  independent  of  the  magnitude  of 
the  step  input  signal.  These  are  typical  responses ,  and  at  all  points  where 
the  effect  of  the  step  input  signal  was  investigated,  the  limit  cycle  was  the 
same  for  different  input  magnitudes.   Comparison  of  Recorder  Trace  VI  =  9  or 
VI-10  with  Recorder  Trace  VT-11  illustrates  the  effect  of  the  relative  posi- 
tion of  the  linear  M- point  and  the  stability  curve  on  the  resulting  limit 
cycle.   Recorder  Trace  VI-11  shows  a  limit  cycle  of  higher  frequency  and  much 
lower  amplitudes  than  the  one  in  Trace  VI-10. 

The  most  striking  result  of  the  initial  computer  study  was  that  the 
unstable  region  of  the  B  versus  B.  plane  can  be  divided  into  three  sub- 
regions.   These  subregions  are  shown  on  the  coefficient  plane  of  Figure 
VI-3.   For  the  linear  M-point  in  Region  I,  only  the  error  channel  saturates 
and  the  limit  cycle  can  be  predicted  from  the  intersection  of  the  K  co- 
ordinate of  the  linear  M-point  and  the  stability  curve.   For  the  linear 
M-point  in  Region  II  both  channels  saturate,  but  the  resulting  limit  cycle 
frequency  always  indicated  an  intersection  of  an  average  M-point  locus  with 
the  stability  curve  at  a  point  to  the  left  of  the  maximum  point  of  the 
stability  curve.   For  the  linear  M-point  in  Region  III.,  only  the  tacho- 
meter channel  saturates  and  the  limit  cycle  can  be  predicted  from  the  inter- 
section of  the  K  coordinate  of  the  linear  M-point  and  stability  curve. 
e 

The  dividing  line  which  proceeds  upward  and  to  the  left  from  the  sta- 
bility curve  will  be  referred  to  as  dividing  line  No.  1.   The  dividing  line 
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which  proceeds  almost  horizontally  to  the  right  from  the  top  of  the 

stability  curve  will  be  referred  to  as  dividing  line  No.  2. 

The  approximate  position  of  dividing  line  No.  1  was  obtained  from  the 

results  tabulated  in  Appendix  B.   To  fix  more  accurately  the  position  of 

the  dividing  line,  more  runs  were  taken  on  the  analog  computer  simulation 

by  choosing  a  K  and  varying  K  .   A  point  on  the  dividing  line  was  then 
e  t 

determined  when  the  amplitude  of  signal  X~  in  the  resulting  limit  cycle  was 
equal  to  3  volts,  the  saturation  voltage.   The  results  are  tabulated  in 
Table  2  of  Appendix  B.   Typical  changes  in  the  limit  cycle  as  K  was  in- 
creased across  No.  1  dividing  line  are  shown  in  Recorder  Traces  VI-12  and 
VI-13.   The  main  feature  to  note  from  these  two  traces  is  the  slight  in- 
crease in  the  amplitude  of  signal  X«  from  just  under  3  volts  to  just  over 
3  volts. 

The  approximate  position  of  dividing  line  No.  2  was  obtained  from  the 
results  of  Table  I  in  Appendix  B.   To  fix  more  accurately  the  position  of 

the  dividing  line,  more  computer  runs  were  made  by  fixing  K  and  varying 

t 

K  .   The  results  are  tabulated  in  Table  3  of  Appendix  B.   Typical  changes 
in  the  response  as  K  was  increased  above  the  dividing  line  are  shown  in 
Recorder  Traces  VI-14  and  VI-15.   Recorder  Trace  VI-14  is  a  couple  of 
oscillations  of  the  limit  cycle  resulting  from  a  linear  M-point  at  K  =  8.0 

and  K  =2.38.   However,  sufficient  time  was  allowed  to  ensure  that  the 

e  ' 

system  did  not  exhibit  the  same  response  as  shown  in  Recorder  Trace  VI-15. 

Recorder  Trace  VI-15  is  a  series  of  three  sections  taken  out  of  the  response 

for  a  linear  M=point  at  K  =8.0  and  K  =2.39  to  show  the  salient  features 

t  e 

of  the  entire  response.   The  time  intervals  between  the  sections  are  in- 
dicated on  the  trace.   Section  one  shows  the  build  up  to  an  apparently 
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stable  limit  cycle  with  only  the  tach  channel  saturating.   However,  closer 
inspection  reveals  that  the  amplitudes  of  signals  X.  and  X~  are  increasing. 
Section  2  starts  at  approximately  the  same  amplitudes  as  occur  at  the  end 
of  Section  1,  but  quickly  builds  up  to  the  final  stable  limit  cycle  as  shown 
in  Section  3  at  an  expanded  time  scale.   The  salient  features  of  these  two 
recorder  traces  are  as  follows: 

1.  The  sudden  build  up  from  an  apparent  limit  cycle  to  a  final 
stable  limit  cycle  for  a  linear  M-point  just  above  dividing  line  No.  2. 

2.  The  large  difference  in  signal  magnitudes  between  the  limit 
cycles  of  the  two  traces  and  the  difference  in  their  frequencies. 

3.  The  change  in  frequency  of  the  oscillations  during  the  build 
up  to  the  final  limit  cycle  in  Recorder  Trace  VI-15.   This  point  is  not 
very  noticeable  on  Trace  VI-15  but  it  is  noticeable  on  traces  with  an  ex- 
panded time  scale  using  dividers  to  measure  the  period. 

4.  The  frequency  and  amplitude  of  Trace  VI- 14  is  very  nearly  the 
frequency  and  amplitude  of  the  apparent  limit  cycle  of  Trace  VI-15. 

Since  the  foregoing  results  allow  one  to  predict  the  limit  cycles  for 
a  large  number  of  linear  system  M-points  for  3  volt  saturations,  the  next 
logical  step  was  to  see  if  the  dividing  lines  are  independent  of  saturation 
level.   Thus  both  error  channel  and  tach  channel  saturation  voltages  were 
adjusted  to  10  volts.   (pots  "a"  and  "a"'  equal  to  0.1).   The  same  analog 
computer  procedures  used  to  determine  the  dividing  lines  for  E   =  E  ~  =  3 
volts  were  again  used.   The  resulting  limit  cycles  for  the  linear  starting 
points  are  tabulated  in  Tables  4  and  5  of  Appendix  B.   A  comparison  of  the 
dividing  lines  for  the  3  volt  and  10  volt  saturation  levels  is  given  below 
in  Tables  VI-2  and  VI-3. 
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TABLE   VI -2 


No.  1  Dividing  Line 

10  volt  Saturation  Level 

3  volt  Saturation  Level 

K 
e 

Kt 

K 

e 

K 
t 

2.1 
2.2 

2.4 
2.5 
3.0 
4.0 
5.0 
6.0 

0.875 
0.855 

0.84 

0.83 

0.795 

0.785 

0.775 

0.77 

2.1 

2.3 

2.5 
3.0 
4.0 
5.0 
6.0 

0.92 
0.89 

0.85 
0.83 
0.81 
0.80 
0.785 

TABLE  VI -3 


No.  2  Dividing  Line 

10  volt  Saturation  Level 

3  Volt  Saturation  Level 

K 
e 

K 
t 

K 
e 

K 
t 

2.445 

2.415 

2.395 
2.385 
2.335 

2.5 

4.0 

6.0 

8.0 

15.0 

2.425 
2.415 
2.405 
2.395 
2.385 
2.335 

3.0 
3.0 
5.0 
6.0 
8.0 
15.0 
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Tables  VI-2  and  VI-3  indicate  that  both  dividing  lines  are  essential- 
ly independent  of  saturation  level,  although  dividing  line  No.  1  is  more 
doubtful  than  dividing  line  No.  2.   In  order  to  definitely  show  this 
point,  more  saturation  voltages  will  have  to  be  used  and  more  sophisticated 
measurement  procedures  should  be  used. 

Practically  all  of  the  data  given  in  Appendix  B  is  graphically  dis- 
played in  Figure  VI-4  through  VI-9.   Figure  VI-4  is  a  plot  of  the  limit 
cycle  amplitude  of  signal  X~  as  K  is  varied  with  K  fixed.   Figure  VI-5 
is  a  plot  of  the  limit  cycle  amplitude  of  signal  X1  for  like  K  variations. 
Figure  VI-6  is  a  plot  of  the  limit  cycle  amplitude  of  signal  X«  as  K  is 
varied  with  K  fixed.   Figure  VI-7  is  a  plot  of  the  amplitude  of  signal 
X  for  like  variations  of  K  .   Figure  VI-8  is  a  plot  of  the  limit  cycle 
frequency  as  K  is  varied  with  K  fixed.   Figure  VI-9  is  a  plot  of  the 
limit  cycle  frequency  as  K  is  varied  with  K  fixed. 

These  six  figures  essentially  show  the  results  as  already  stated.   By 
correlation  of  these  figures  with  the  location  of  the  linear  M-point  in 
Figure  VI-3,  one  can  see  the  regions  in  which  only  one  element  saturates. 
For  example,  the  upper  right  portion  of  Figure  VI-8  shows  Region  III  of 
Figure  VI-3  in  which  one  would  expect  the  frequency  to  remain  constant  as 
K  is  increased.   Also,  the  transition  across  dividing  line  No.  2  is  very  ap' 
parent  in  all  the  figures  because  of  the  large  discontinuity  in  the  curve. 
Some  curves  for  relatively  small  K  and  K  abruptly  stop  and  then  start  up 
again  in  a  different  portion  of  the  graph.   This  is  due  to  the  variable 
being  inside  the  region  of  stable  linear  M-points  in  Figure  VI-3,  thus  no 
limit  cycle.   An  important  result  to  note  which  has  not  yet  been  pointed  out 
is  that  for  K  greater  than  3  and  K  greater  than  3,  the  frequency  of  the 
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resulting  limit  cycle  is  essentially  independent  of  variations  in  either 

K  or  K  . 
e     t 

Figure  VI-10  shows  various  instantaneous  and  "average"  quantities  plotted 

against  the  number  of  elapsed  oscillation  peaks  for  a  linear  M-point  just 

above  No.  2  dividing  line.   A  better  abscissa  name  would  be  "number  of 

elapsed  half  cycles"  although  both  are  related  to  the  elapsed  time  by  the 

characteristics  of  the  system.   The  quantities  plotted  in  the  figure  are 

defined  below: 

B.min  =   The  minimum  instantaneous  value  during  a  half  cycle 

of  the  s   coefficient  calculated  by  equation  Vl-d 

using  the  minimum  instantaneous  gain  N«  .  .   N„  . 

2min    2min 

occurs  when  the  signal  X«  is  maximum  during  the  half 


1 


lave 


omtn 


oave 


cycle. 

The  "average"  value  of  the  s1  coefficient  over  a  half 

cycle  as  calculated  from  equation  Vl-d  using  the  "average" 

gain  of  the  nonlinearity  N_    .  The  describing  function 
°  J      2ave 

was  used  to  relate  N0    to  the  maximum  value  of  X« 

2ave  2 

during  the  half  cycle. 

The  minimum  instantaneous  value  of  the  s°  coefficient 
calculated  by  equation  VT-c  using  the  minimum  instan- 
taneous gain  N,  .  .   N,  .   occurs  when  the  signal  X, 
°     lmin    lmin  °     1 

is  maximum  during  the  half  cycle. 

The  "average"  value  of  the  s°  coefficient  calculated 

from  equation  VI- d  using  the  "average"  gain  of  the  non- 

linearity  N,    .   The  describing  function  was  used  to 
lave 

relate  N,    to  the  maximum  value  of  X,  during  the  half 
lave  1 


cycle. 
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CO     ■   *  The  average  frequency  of  the  oscillation  over  a  half 

oVc 

cycle  obtained  by  measuring  the  half-period  of  the 
oscillations. 
Justification  for  using  the  describing  function  in  the  manner  indicated 
above  will  be  indicated  in  Chapter  IX. 

If  the  "average'1  M-point  is  assumed  to  indicate  the  stability  of  the 
system  over  the  half  cycle,  then  it  is  interesting  to  note  that  the  "aver- 
age" M-point,  as  defined  by  B,    and  B    ,  moves  through  a  portion  of 

lave      oave 

the  stable  region  on  a  constant  B  line.   It  is  also  noted  that  under  limit 

o 

cycle  conditions  at  the  extreme  right  side  of  Figure  VI-10,  the  values  of 

B,  and  B  calculated  using  U)    in  Mitrovic's  stability  curve  equations 
1      o  °   ave 

are  approximately  equal  to  the  end  values  of  B,    and  B     curves  which 

^  lave  oave 

were  calculated  using  the  describing   function  and  the  amplitude  of  the   limit 
cycle.     A  comparison  is  made  below. 

From  Stability  Curve  From  Figure  VI-10 


Bx  =  30.9 

B     -  43.2 
o 


B.  =  31.4 

lave 

B  =  42.8 

oave 


+3!Mfe  — 


-3V«lf$ 


FIGURE     TI-2L 
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15  RUSH      INSTRUMENTS 


DIVISION  OF  CLEVITE  CORPORATION  CLEVELAND.  OHIO 


m 


Recorder  Trace  VI- 1 
Linear  M-Point  Close  to  Left  Side  of  Stability  Curve, 


Ke  =  2.4,  K.  -  1.7,  1.2  Volts  Input  Step 

t 

SCALES : 


Xt  =  0.2  Volts/line 


Y  -  0.2  volts/line 


X2  -  0.2  volts/line 


Y2  -  0.2  volts/line   Time  -  5  Seconds/Division 
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Recorder  Trace  VI-2 

Linear  M-Point  Close  to  Left  Side  of  Stability  Curve,  K  =  2.4 

e  , 

K  =  1.7,  62.5  Volts  Input  Step 

SCALES : 
X  -  10  Volts/Line 

-0.05  Volts/Line 

Y?  -  0.2  Volts/Line 
-0.1  Volts/Line 
-0.05  Volts/Line 

Time  -  5  Seconds/Division 
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Y  -  0.2  Volts/Line 
-0.05  Volts/Line 


X2  -  0.5  Volts/Line' 

-0.1  Volts/Line 
-0.05  Volts/Line 


Recorder  Trace  VI- 3 
Linear  M-point  Close  to  Right  Side  of  Stability  Curve, 

K  =  1.0.  K  =  3.9,  2.6  Volts  Input  Step. 
e        t 


SCALES : 

X  -0.2  Volts/Line 
-0.01 


Y.  -0.2  Volts/Line    Y2-0.2  Volts/ 
-0.01  Volt  s/Line  Line 


X2-0.2  Volts/Line 


Time  -  5  Seconds/Division 
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Recorder  Trace  VI-4 

Linear  M-point  Close  to  Right  Side  of  Stability  Curve,  K  =  1.0, 
K  =  3.9,  95  Volts  Input  Step. 


SCALES 

X.  -  5Volts/Line 
-0.05  Volts/Line 

X2  -  0.3  Volts/Line 


Y.  -  0.2  Volts/Line 


Y2-0.2  Volts/Line 
-0.05  Volts  Line 


Time  -  5  Seconds/Division 
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Recorder  Trace  VI-5 

Linear  M-point  in  Well  Damped  Region,  Kg  =  0.5,  K  =  2.0,  5.4  Volts 
Input  Step. 


SCALES : 

X  -0.1  Volts/Line 


Y.-0.2  Volts/Line 


X2-0.1  Volts/Line 


Y„-0.2  Volts/Line     Time  -  5  seconds/Division 
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Recorder  Trace  VI-6 
Linear  M-point  in  Well  Damped  Region,  K  =  0<,5S  K  =  2.0 
60  Volt  Input  Step. 
SCALES : 
Xx-  1  Volt/Line     Y  -0.2  Volt/Line   X2-0.2  Volt/Line 

Y.-0.2  Volt/Line    Time  -  5  Seconds/Division 
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Recorder  Trace  VI- 7 

Repeated  Sisturbances,  M-Point  at  K  =  2.3,  K  =  1.5 

SCALES : 
X.  -0.2  Volts/Line     Y.-0.2  Volts/Line     X2-0.5  Volts/Line 

Y2-0.2  Volts/Line   Time  -  5  Seconds/Division 
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Recorder  Trace  VI-8 


Repeated  Disturbances,  Linear  M-Point  at  K  =  2.3S  K  =  2.0 

6  t 

SCALES : 

X  -0.2  Volts/Line    Y^O.2  Volts/Line   X2~0.5  Volts/Line 

Y?-0.2  Volts/Line    Time  -  5  Seconds/Division 
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RECORDER   MARK  II         chart  no.  RA2921  32 


Recorder  Trace  VI-9 
Unstable  Linear  M-point  at  K  =  39  K  =  4,  One  Volt  Step  Input 

SCALES : 


X.-0.5  Volts/Line 
Y2-0.2  Volts/Line 


Y.-0.2  Volts/Line     X2~l  Volt/Line 

Time  5  Seconds/Division 
-0.2  Second/Division 
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Recorder  Trace  VI- 10   

Unstable  Linear  M-Point  at  K  =  3,  K  =  4,  10  Volt  Step  Input. 


SCALES : 

X  =  2Volts/Line 
1   -0.5  Volts/Line 


Y.-0.2  Volts/Line 


X„-l  Volt/Line 


Y2-0.2  Volts/Ltne 


Time  -  5  Seconds/Division 
-0.2  Seconds/Division 
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Recorder  Trace  VI-11 

Unstable  Linear  M-Point  Close  to  Stability  Curve  at  K  =  1.5,  K  =  4.5 

e        t 


X  -  1  Volt/Line 
-0.05  Volts/Line 


Y  -0.2  Volts/Line 
-0.05  Volts/Line 


X2-0.5  Volts/Line 


Y2-0.2  Volts/Line 


Time  -  5  Seconds/Division- 
-0.2  Volts/Division 
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SCALES 


Recorder  Trace  VI-12 

Limit  Cycle  for  Unstable  Linear  M-Point  to  Left  of  No.  1 

Dividing  Line,  K  =  4,  K=  0.8 
e       t 

Xx    -0.5  Volts/Line   Yj-0.2  Volts/Line    X2~0.2  Volts/Line 

Y2-0.2  Volts/Line    Time  -0.2  Seconds/Division 
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Recorder  Trace  Vl-13 

Limit  Cycle  for  Unstable  Linear  M-Point  to  the  Right  of  No.  1 
Dividing  Line,  K  =  4S  K  ■  0.82. 

SCALES : 

X,  -0.5  Volts/Line     Y,-0.2  Volts/Line    Xo-0.2  Volts/Line 
1  \  I 

Y2-0.2  Volts/Line    Titne-0.2  Seconds/Division 
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DIVISION  OF  CLEVITE  CORP 
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Recorder  Trace  VI- 14 

Limit  Cycle  for  Unstable  Linear  M-point  Below  No.  2 

Dividing  Line  at  K  =  2.38,  K  =8. 

e         t 

SCALES : 

X1-0.1  Volts/Line   Y^O.2  Volts/Line    X^l  Volt/Line 
Y2-0.2  Volts/Line    Time-0,2  Volts/Line 
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Recorder  Trace  VI-15 

Unstable  Linear  M-Point  Above  No.  2  Dividing  Line  at  K  =  2  40  K  =8 

e       s   t 

SCALES : 

X^O.2  Volts/Line    Y^O.2  Volts/Line     X2~2  Volts/Line   Y2~0,2  Volts/Line 

Time  -  5  Seconds/Division 

-0.2  Second/Division  -  — — 
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CHAPTER  VII 
EXPLANATION  OF  Th.     ITS  AND  SOME  ME  I  OF 

PREDICT!  -•  IS 

An  attempt   is  ma  li  hiapter  to  explain  the  syst--      -  sp 

•  VI  by  a  ■■■ .         .  btained  in 

•  -    M=poi  '  ion  on  plan« 

. 

i  -■  sxplan  \    ..       f   of  the  results   -  t  -Ipful 

Ight  in  ens  can  t   obtained  by  first  looking  at  the  system  as 

'       --    ion  and  then  looking  at  the  system  as  only 

having  1  saturation.   First 9  assuming  only  err-    annel  satura 

ipplying   •  !  arguments  given  in  Chapter  19  Section  I-A-l  to  Fig- 

••  :>  t   VI-  nined  fox       <C  4„5  the  error  saturation  produces  a 

stable  li:'.       ,   For  K  y   4,5  the  linear  system  can  never  be  made  stable 

fey  varying  K  ;  consequently,,  there  is  no   possibility  of  a  limit  cycle  wi 

erros  channel  saturation  only  and  1  y    4050  This  fact  is  easily  shown  by 

t 


root  locus  analysis.   The  root  locus  of  the  inner  loop  is  sketched 
in  Figure  VII=10   For  a  F.  large  enough  to  place  the  inner  loop  closed  loop 
t    -     he  right  half  planea  no  amount  of  K  can  make  the  root  locus  swing 


■  .  -    me,  I    i'.rting  with  the  closed  loop  roots  of  the  ii 

loop  as  ■        p    a   loop  roots  of  the  outer  loop9  the  assymptotes  remain  ati60' 

I  is  moved  to  the  righto   A  rough 

the  angle  ®  from  the  upper  pole  shows  i  hat  this   g    §  a  fi  • 

quads  yi/j     axis  corssover  of  the  inner  loop  root  lo  s  Ls 

determined  to  be  F  =  4„5»   Thus  the  outer  loop  root  locus  can  only  precede 

t  '    " 

as   indicated   in  Eigui  K  5*    4„5  the   system  is  always    .         sfele. 


Outer   toop    Root  locus 


OPBh)    LOOP    POLE      OF  OUTER    LOOP 


XNNER  LOOP    f?OOT    LOCUS 


. < ^ 

-4 

INNER  LOOP  Asy^PTOTF 


OUTER   LOOP   ASYMPTOTE 


Figure  VII- 1 
Root  Locus  of  the  Tachometer  Feedback  Path 


thHn  the  Position  Feedback  Path. 
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Assuming  only  tach  channel  saturation  and  applying  the  arguments 
of  Chapter  Vs  Section  V-A-l  to  Figure  VI-3,  it  is  determined  that  for  a  K  ^ 
2.45; 

1.  The  side  of  the  stability  curve  to  the  right  of  the  maximum 
(  U)<\  ^  2.645)  yields  stable  limit  cycles. 

2.  The  side  of  the  stability  curve  to  the  left  of  the  maximum 

(  U)fl  ^    2.645)  yields  unstable  limit  cycles,  at  least  for 
0.75  <   Ke<  2.45. 

3.  Which  limit  cycle  governs  the  response  is  entirely  dependent 
upon  the  input  signal  magnitude. 

For  K   >   2.45  Figure  VI-3  indicates  the  linear  system  is  always  unstable 
no  matter  what  the  K  ,  thus  no  limit  cycle  is  possible.   For  K  K.     0.75, 
the  absolute  minimum  gain  of  the  nonlinearity,  i.e.,  N~  =  0,  the  M-point 
is  still  within  the  stable  region,  and  thus  an  unstable  limit  cycle  is 
impossible. 

From  the  above  discussion,  one  can  safely  say  that  with  both  satura- 
tions present  a  stable  limit  cycle  will  always  result  no  matter  where  the 
linear  M-point  lies,  even  though  individually  each  saturation  has  certain 
areas  in  which  the  single  nonlinearity  system  is  absolutely  unstable.   This 

statement  is  true  because  for  large  K  and  K  the  tach  channel  saturation 

°   e      t 

moves  the  "average"  M-point  to  the  left  such  that  the  error  channel  satura- 
tion can  move  it  downward  to  lie  on  the  stability  curve,  or  vice  versa. 

A  result  of  the  Chapters  of  Part  I  will  also  be  helpful  in  being  able 
to  predict  the  results  of  Chapter  VI.   The  magnitude  of  the  describing  func- 
tion at  the  intersection  at  the  intersection  with  the  linear  loop  transfer 
function  curve  on  the  gain-phase  plane  is  essentially  equal  to  the  magni- 
tude of  the  "average"  gain  of  the  nonlinearity  as  obtained  by  Mitrovic's 
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and  root  locus  methods.   Table  VII-1  gives  a  comparison  of  N    and  G^ 

°  ave      D 

under  the  limit  cycle  condition  for  all  the  systems  analyzed  in  Part  I. 


TABLE  VII-1 


Nonlinearity 

No.  of 
Zeros 

Root  Locus 

N 
ave 

Mitrovic1 s 

N 
ave 

Describing 

Function(G^) 
D7 

Saturation 

Dead  Zone 

Ideal  Relay 

Relay  with 
Dead  Zone 

Tach  FB 
Saturation 

Acceleration 
FB  Saturation 

0 

1 

0 

1 

0 

1 

0 

1 

0 
0 

0.60 
0.8051 

0.5951 
0.4089 

1.1053 
0.5573 

0.4568 
0.6208 

0.08737 
0.5958 

0.4614 

0.60 
0.8052 

0.590 
0.4072 

1.1032 
0.5573 

0.4543 
0.6208 

0.08739 
0.5959 

0.4614 

0.604 
0.7953 

0.5989 
0.3987 

1.0976 
0.5629 

0.4225 
0.5757 

0.08802 
0.5957 

0.4648 

The  number  of  significant  figures  listed  for  the  describing  function 
in  Table  VII-1  is  misleading  since  G  was  calculated  using  at  best  three 
figure  accuracy  from  a  graphical  plot  of  an  intersection.   Even  with  this 
uncertainty  in  G  s  the  comparisons  only  differ  greatly  in  the  last  two  or 
three,  significant  figures.   With  more  accurate  calculation  of  the  inter- 
sections of  the  three  methods,  the  numbers  should  be  exactly  equal.   A 
non-rigorous  proof  of  this  equality  follows: 

The  describing  function  method  applies  the  Nyquist  criterion  and  is 
basically  a  frequency  response  technique  involving  a  variable  system  criti- 
cal point.   For  this  proof  assume  that  the  symbol,  G  ,  and  the  words,  des- 
cribing function,  are  undefined  in  terms  of  the  amplitude  of  the  input 
signal  to  nonlinearity.   Then  for  gain  variable  nonlinearities,  the  Nyquist 
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criterion  becomes 

gd  GUu)}  =  -i  (vii-i) 

This  general  equation  can  then  be  applied  to  any  one  of  the  systems  anal- 
yzed in  Part  I. 

For  convenience  apply  the  root  locus  criterion  to  the  system  of  Section 
I-As  Chapter  I  which  yields 


10  N  =  -1 


s(s  +  l)(s  +  2) 

But  in  equation  VXI-2 


(VII-2) 


. 10 =  G(s) 

s(s  +  l)(s  +  2) 

which  is  the  linear  loop  transfer  function.   Thus 

NG{s)=  -  1  (VII- 3) 

where  G(s)  can  now  be  any  linear  loop  transfer  function. 

In  the  root  locus  method  of  predicting  limit  cycles  for  gain  variable 
nonlinearitieSj  one  is  only  interested  in  the  limit  of  stability  of  the 
linear  loop  transfer  function.   Another  way  of  saying  the  same  thing  is  that 
one  is  interested  in  finding  the  point  on  the  linear  system  root  locus  at 
which  self  sustained  oscillations  are  possible.   Of  course,  this  occurs  when 
the  linear  root  locus  crosses  the  JtO  axis  or  when  s  =  ^oJ)  .   Sub- 
stituting ^U)  for  s  in  equation  VII-4, 

NG(  ^U)  )   =  -1.  (VII-4) 

It  is  a  well  known  fact  that  both  the  frequency  response  technique 
and  the  root  locus  accurately  predict  the  absolute  stability  of  a  system. 
Consequently;,  equations  VII-1  and  VI-4  must  be  equivalent  and  therefore  N  and 
G  must  be  equal  since  the  G(^LO  )  is  the  same  in  both  equations  for  any  given 
system. 
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The  fact  that  N    and  G„   are  equal  under  limit  cycle  conditions  is 
ave      D 

really  only  useful  when  G  as  a  number  is  related  to  the  amplitude  of  the 
input  signal  to  the  nonlinearity.   Of  course  this  relationship  is  strictly 
valid  only  when  the  input  signal  is  sinusoidal.   Thus  when  the  use  of  the 
describing  function  can  be  justified  in  a  given  system,  root  locus  and 
Mitrovic's  methods  will  also  accurately  predict  the  limit  cycle  amplitude. 
The  limit  cycle  can  then  accurately  and  quickly  be  determined  by  substitut- 
ing s  «  3U3  in  the  characteristic  equation,  calculating  the  amplitude 
from  the  describing  function  relationship. 

PROBABLE  QUALITATIVE  EXPLANATIONS 
1.   Stable  responses  to  step  inputs  for  which  the  linear  system 
M-point  lies  within  the  stability  region. 

a.   For  small  signal  inputs  in  which  the  error  channel  does 
not  saturate,  Recorder  Traces  VI-1,  VI-3,  and  VI-5  are  representative. 
Two  of  these  traces  show  the  first  peak  of  signal  X,,  is  large  enough  to 
saturate  the  tach  channel.  However,  after  the  first  velocity  peak,  the  X- 
and  X~  signals  are  well  below  the  saturation  levels  of  3  volts  and  the 
oscillations  die  out  according  to  the  position  of  the  linear  M=point.   In- 
sight into  an  explanation  is  obtained  if  the  instantaneous  minimum  M-points 
are  calculated  from  X«    of  the  first  peak  in  Traces  VI-3  and  VI-5  and  then 

plotted  on  the  B  and  versus  B.  plane.   These  instantaneous  M-points  are  shown 
o  1 

as  M.  and  M_  respectively  in  Figure  VII-2.  The  starting  point  of  the  arrow 
is  the  linear  M-point  and  the  head  of  the  arrow  is  the  instantaneous  mini- 
mum M-point.   One  trace  which  was  not  included  in  the  results  of  Chapter  VI 
showed  that  the  first  peak  of  X~  saturated  the  tach  channel  sufficiently  to 
drive  the  instantaneous  minimum  M-point  into  the  unstable  region  as  shown 

232 


by  M  in  Figure  VII-2„   However,  the  tachometer  channel  did  come  out  of 

saturation  and  the  system  was  stable.  Thus  for  linear  M-points  close  to 

the  left  portion  of  the  stability  curve  the  tach  channel  may  saturate  but 

the  magnitude  of  signal  X~  is  not  large  enough  to  cause  the  instantaneous 

M-point  to  be  in  the  unstable  region  for  a  sufficient  portion  of  the  first 

oscillation  to  drive  the  average  M-point  to  the  stability  curve.   This  was 

thought  to  be  due  to  the  relatively  small  values  of  K  .   As  K  is  increased 

in  an  attempt  to  raise  the  magnitude  of  the  first  peak  of  signal  X«»  two 

effects  combine  to  keep  the  response  stable.   First  as  K  increases,  the 

linear  M-point  moves  to  the  right  and  the  linear  system  damping  is  increased, 

at  least  the  damping  is  increased  until  the  M-point  moves  approximately  across 

the  axis  of  the  parabolic  stability  curve.   The  increased  damping  tends  to 

nullify  the  expected  increase  in  X~    .   Secondly  and  most  important,  the 

relative  distance  the  average  M-point  must  move  to  reach  the  left  side  of 

the  stability  curve  increases.  The  left  portion  of  the  stability  curve  is  the 

governing  stability  limit  for  the  linear  M-point  in  the  stable  region  because 

only  by  tach  saturation  can  the  system  be  driven  unstable.   Thus,  the  increase 

in  X„    is  not  sufficient  to  offset  the  increase  in  this  relative  distance 
2max 

such  that  the  instantaneous  M-point  can  lie  in  the  unstable  region  for  a 
sufficient  portion  of  the  first  oscillation  to  drive  the  average  M-point  to 
the  stability  curve. 

As  the  linear  M-point  moves  to  the  right  past  the  approximate  axis  of  the 
stability  curve  the  damping  begins  to  decrease  but  the  relative  distance  be- 
tween linear  M-point  and  the  left  side  of  the  stability  curve  is  quite  large 
and  obviously  overrides  the  increase  in  peak  X~  and  keeps  the  system  stable. 
This  relative  distance  is  in  some  manner  related  to  the  degree  of  saturation 
necessary  to  cause  limit  cycle  operation  since  this  is  the  graphical  distance 
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the  average  M-point  must  move  to  reach  the  unstable  limit  cycle  governing 
the  tach  channel. 

b.   For  large  input  steps  in  which  the  error  channel  in- 
itially saturates,  the  same  explanations  given  above  also  apply  after  the 
error  channel  comes  out  of  saturation.   Recorder  Traces  VI-2,  VI-4,  and 
VI-6  are  applicable.   It  is  constructive,  however,  to  analyze  the  response 
prior  to  the  error  channel  becoming  linear.   The  effects  to  be  discussed 
are  most  apparent  in  Recorder  Trace  VI-4.   With  the  error  channel  initially 
saturated;,  the  initial  driving  voltage  of  the  plant  is  3  volts.   Since   X~ 
does  not  change  instantaneously  (it  is  proportional  to  velocity)  the  driv= 
ing  voltage  remains  at  3  volts  for  a  short  instant  of  time.   The  trace  in- 
dicates that  X^,  increases  very  rapidly  and  that  Y«  essentially  jumps  to  the 
saturation  value  of  3  volts.  Also,  the  output  angle  and  X.  have  not  changed 
during  this  short  time  interval.   Thus  because  of  the  signs  of  the  feedback 
voltages,  the  driving  signal  of  the  plant  goes  to  zero  after  a  very  short 
time.   This  leads  to  the  initial  plant  driving  signal  shape  shown  in  Figure 
VII-4a  as  a  first  approximation.   Second  approximations  might  be  as  shown 
in  Figure  VII-4b  and  c.   The  solution  for  the  output  angle  when  the  first 
signal  approximation  is  applied  to  the  plant  is  of  the  form, 

0o  =  Clt  +  Cie"t+c5e"2t-+Qer4t+C)S^-t,)  +C6  e"(tr"*'^ 
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Figure  VII-2 


M-point  Loci  for  Stable  Linear  Systems  for  Small  Step  Input 

Magnitude 
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Figure  VII- 3 
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M-point  Loci  for  Stable  Linear  Systems  for  Large  Step  Input 

Magnitude 


Differentiation  gives  the  output  velocity  as 

•     „      i     -t     ^    *     -It       -   i    -4t   „    ~  I  -ft-t,) 


<► 


Since  X?  =  K^-  ©<p    ,  the  initial  portion  of  signal  X?  is  easily  explained 
by  the  sum  of  the  expontential  terms  given  above.   Time  t.  is  the  time  at 
which  the  tach  channel  saturates.  At  time  t2s  indicated  on  Recorder  Trace 
VI-4,  the  response  of  the  plant  to  the  initial  input  has  ended  and  the 
system  is  coasting  because  the  plant  driving  signal  is  zero.   At  time  t  , 
the  velocity  of  the  coasting  system  has  decreased  such  that  X«  is  less  than 
3  volts  and  the  tach  channel  becomes  linear.  As  X«  slowly  decreases  below 
3  volts9  a  slowly  increasing  positive  voltage  is  applied  to  the  plant  be- 
cause the  error  channel  is  still  transferring  a  positive  3  volt  signal. 
This  is  indicated  on  the  Recorder  Trace  by  the  fact  that  X~  levels  out  at 
about  2.8  volts  instead  of  remaining  in  an  expontential  decay.  The  relative- 
ly sudden  decrease  in  X0  after  time  t.  and  before  X„  levels  out  can  not  be 
j  2  3  2 

explained,  except  by  noting  that  the  saturation  simulations  do  not  strictly 
limit  their  outputs  to  three  volts.   This  is  due  to  a  finite  diode  resist- 
ance and  pot  loading  effects  when  the  diode  conducts.   Signal  Y.  is  initial- 
ly 3.6  volts  and  signal  Y~  is  3.1  volts  after  saturation.   Thus  the  plant  of 
the  analog  cimulation  is  actually  being  driven  by  about  0.5  volts  or  less 
during  the  period  in  which  the  author  has  assumed  zero  driving  voltage.   This 
fact  may  account  for  the  relatively  sudden  decrease  in  X~  after  time  t_. 
Eventually  the  output  drives  to  within  3  volts  of  the  input  and  the  error 
channel  comes  out  of  saturation  at  time  t,  .   From  this  time  on  the  system 
responds  according  to  the  linear  M- point. 
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Figure  VII-4 
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Possible  Waveshapes  of  Initial  Plant  Driving  Signals  for  Stable  , 
Linear  M-point  with  Large  Input  Steps  Applied.  > 

The  same  response  in  Recorder  Trace  VI-4  can  also  be  qualitatively  ex- 
plained strictly  by  the  system  M-point  motion.   With  the  error  channel  in- 
itially saturated  the  M-point  immediately  jumps  downward  along  the  linear 
M-point  K  coordinate.   For  Recorder  Trace  VI-4,  this  initial  M-point  jump 

moves  the  M-point  to  K  =  0.0316  or  B  =  0.632  by  equation  VI-c.   This  M- 

e  o 

point  motion  increases  the  system  damping  a  small  amount  over  the  damping 
of  the  linear  system.   Within  a  very  short  time  interval,  the  tach  channel 
saturates  and  rapidly  moves  the  M-point  to  the  left,  essentially  along  a 
constant  K  line  since  the  error  has  not  appreciably  changed  to  cause  an 
upward  component  in  the  M-point  motion.   The  rapid  movement  of  the  M-point 
to  the  left  greatly  increases  the  system  damping,  possibly  close  to  criti- 
cal damping.   With  the  M-point  moving  into  a  highly  damped  region  within  a 
short  time,  very  few  oscillations  would  be  expected  and  would  be  damped  out 
very  quickly  by  say  time  t2  in  Recorder  Trace  VI-4.   The  peak  velocity  is 
more  difficult  to  explain  but  is  certainly  dependent  upon  how  soon  the  tach 
channel  saturates  since  this  is  the  controlling  factor  of  the  time  required 
for  M-point  motion  into  the  highly  damped  region.   After  time  t«,  it  is 
difficult  to  explain  the  response  strictly  by  M-point  motion  until  time  t. 
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is  reached  when  the  error  channel  comes  out  of  saturation.   At  this 
time  the  system  is  completely  linear  and  the  oscillations  occur  according 
to  the  *o  and  L&r)    of  the  linear  M-point.   These  oscillations  will  be  small 
since  the  disturbance  as  the  error  channel  comes  out  of  saturation  is 
small.   The  estimated  M-point  motion  corresponding  to  Recorder  Trace  VI-4 
is  sketched  in  Figure  VII-3  along  with  those  of  Recorder  Traces  VI-2  and 

VI-6.   The  corresponding  times  from  recorder  trace  VI-4  are  indicated  along 

.-•rv' 
the  path.   The  lines  of; constant  U    were  not  calculated  but  their  shape  was 

sketched  from  those  calculated  for  a  different  fourth  order  system.   For  the 

given  fourth  order  system,  the  constant  3  lines  will  not  differ  greatly 

from  those  sketched.   The  objective  of  the  constant  ^   lines  is  to  show 

that  the  instantaneous  damping  can  increase  or  decrease  either  a  small 

amount  or  a  large  amount  for  both  horizontal  and  vertical  motions  of  the 

T 
system  M-point.   he  entire  figure  is  not  to  scale  but  the  numbers  listed 

are  correct. 

c.   The  fact  that  the  system  could  be  driven  into  a  limit 

cycle  by  repeated  disturbances  and  linear  M-points  close  to  the  left  portion 

of  the  stability  curve  cannot  be  explained  satisfactorily  by  M-point  motion. 

However 9  from  Recorder  Trace  VI-8,  it  was  noted  that  between  disturbances 

and  after  the  first  saturated  peak  of  X~,  the  signals  X-  and  X~  were 

slightly  damped  and  below  3  volts  but  at  a  generally  higher  amplitude  than 

before  the  disturbance.   With  little  linear  system  damping  the  magnitude  of 

X~  oscillations  decreases  very  little  between  disturbances.  After  several 

disturbances  the  X^  oscillations  are  just  below  the  saturation  level.   The 

next  disturbance  causes  the  first  peak  of  X~  to  saturate  N„  to  the  degree 

required  to  move  the  average  M-point  from  the  linear  point  to  the  stability 

curve  and  the  oscillations  build  up  since  tach  saturation  gives  an  unstable 
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limit  cycle.   The  final  limit  cycle  is  stable  because  the  error  channel 
eventually  saturates  as  the  oscillations  build  up.   It  was  stated  in 
Chapter  VI  that  the  magnitude  of  the  disturbance  had  to  be  small  such  that 
the  error  channel  does  not  saturate  appreciably,  otherwise  the  magnitude  of 
X?  oscillations  after  the  disturbance  was  less  than  the  magnitude  of  the 
oscillations  before  the  disturbance.   This  result  is  easily  explained  by 
the  M-point  motion.   When  the  error  channel  saturates  the  M-point  moves 
rapidly  downward,,  greatly  increasing  the  damping.   The  increased  damping 
damps  out  the  oscillations  before  the  output  has  time  to  come  into  corres- 
pondence with  the  new  input  demand.   Consequent lys  the  disturbance  causing 
the  oscillations  after  the  error  channel  comes  out  of  saturation  is  much 
less  than  the  disturbance  which  caused  the  error  channel  to  saturate.   The 
response  is  very  similar  to  the  response  to  a  large  step  input  applied  at 
time  zero  as  described  previously. 

2„   Limit  cycles  were  obtained  only  for  linear  M-points  in  the 
unstable  region. 

a.  The  result  that  the  limit  cycles  were  independent  of  the 
magnitude  of  the  input  signal  was  suspected  because  of  the  arguments  given 
in  the  introduction  to  Part  I.   For  saturation  in  the  error  channel  limit 
cycles  are  only  possible  for  an  unstable  linear  system.   Thus  the  existance 
of  the  limit  cycle  is  entirely  dependent  on  the  linear  system  and  not  the 
step  input  signal  because  the  linear  poles  and  zeros  determine  the  possibit 
ity  of  instability  and  the  linear  gain  determines  whether  or  not  the  linear 
system  is  stable  or  unstable.   For  saturation  in  the  tach  channel,  an  un- 
stable limit  cycle  is  possible  with  a  stable  linear  system.   However,  the 
existance  of  this  limit  cycle  is  only  dependent  on  the  linear  system.   The 
magnitude  of  the  step  input  signal  determines  only  whether  or  not  the  non- 
linear system  response  is  unstable  or  stable.   Thus  if  the  existance  of  a 
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limit  cycle  for  each  of  the  nonlinearities  individually  is  indepdent  of 
the  magnitude  of  the  input  step  function,  then  the  existance  of  a  limit 
cycle  with  both  nonlinearities  in  the  same  system  is  independent  of  the  step 
input . 

b.  The  result  that  the  limit  cycle  amplitudes  of  X,  and  X2 
are  dependent  on  the  relative  distance  between  the  linear  M-point  and  the 
stability  curve  is  explained  by  degrees  of  saturation.   Degree  of  satura- 
tion means  the  relative  magnitudes  of  N,    and  N0    which  are  required 

°  lave      2ave 

to  place  the  average  M-point  on  the  stability  curve.   For  systems  with 

linear  M-points  close  to  the  stability  curve,  the  average  M-point  has  a 

small  distance  to  travel  either  horizontally  and/or  vertically  to  reach  the 

stability  curve.   For  sinusoidal  input  waveshapes  of  signals  X.  and  X?,  it 

has  been  shown  that  N    can  be  related  to  the  amplitude  of  the  input  signal 

ave 

through  the  describing  function.   For  small  average  M-point  motions  N, 

lave 

and  N„    will  be  just  less  than  1.0  from  equations  VI-c  and  Vl-d  because 
2ave  n 

B     and  B,    will  be  very  close  to  the  linear  M-point  values.   By  the 
oave      lave 

describing  function  relation  for  values  close  to  one,  the  input  amplitude  is 
just  slightly  greater  than  the  saturation  voltage.   Thus  the  amplitudes  of 
X1  and  X«  will  be  slightly  greater  than  3  volts.   For  linear  M-points  re- 
latively distant  from  the  stability  curve,  the  values  of  B     and  B. 

oave      lave 

will  be  greatly  different  from  the  linear  M-point  coordinates.   Thus  N, 

lave 

and  N~    are  small  fractions  and  through  the  describing  function  relation- 
ship the  amplitudes  of  X,  and  X«  are  correspondingly  large. 

c.  The  result  that  for  K  /^  3  and  K   ^3  the  frequency  of 
the  limit  cycle  is  essentially  independent  of  K  and  K  variations  can  be 
explained  only  in  terms  of  other  results  which  are  apparent  in  Figures  VI-4 
through  VI-7.   From  Figure  VI-5  it  is  noted  that  for  constant  K  above  the 
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stability  curve  and  for  K  ^2,  the  amplitude  of  X  is  independent  of  K 

c  1  c 

variations  and  only  dependent  on  K  variations.   Thus  as  K  increases  above 
K  =  2,  the  vertical  distance  the  average  M-point  must  move  to  reach  the 
stability  curve  is  relatively  constant  for  K  .constant.  This  fact  can  be  ex- 
plained as  follows.   No  matter  what  the  K  or  K  ,  the  initial  driving 
voltage  of  the  plant  remains  at  3  volts,  thus  the  initial  changes  in  ©0 
and  ©0  are  unaffected  by  K  and  K  .   Also   ©o  *s  subject  to  an  extra 
time  lag  which   ©o    is  not.   Thus  for  K  ">•   K  >  3     »  it  can 
safely  be  said  that  the  tach  channel  saturates  before  the  error  channel  and 
the  resultant  average  M-point  motion  is  to  the  left  instead  of  downward. 
Consequently,  it  is  difficult  to  concieve  of  a  limit  cycle  occurring  on  the 
right  side  of  the  stability  curve.   This  fact  is  shown  in  experimental  re- 
sults since  all  the  limit  cycles  obtained  from  linear  M-points  in  Region  II 
exhibited  frequencies  on  the  left  side  of  the  stability  curve.   Thus  the 
tach  channel  first  saturates  to  move  the  average  M-point  to  the  left  past 
the  maximum  of  the  stability  curve  along  a  constant  K  line.   Once  the 
average  M-point  is  in  this  position  the  vertical  distance  it  must  move  to 
reach  the  stability  curve  by  error  channel  saturation  is  relatively  con- 
stant as  long  as  the  ^coordinate  of  the  linear  M-point  is  kept  constant. 
As  K  increases,  Figure  VI-7  shows  that  for  K  >   3,  the  limit  cycle  ampli- 

K  6 

tude  increases  linearly  with  K  and  that  K  variations  have  little  effect 

e  t 

when  K   >   2.   The  linear  variations  of  X,    with  K  is  reasonable  since 
t  lmax       e 

it  has  been  shown  that  the  relative  vertical  distance  the  average  M-point 

must  move  is  independent  of  K  for  a  given  K  .   Thus  the  B     coordinate  of 

t  e  oave 

the  average  M-point  of  the  limit  cycle  is  relatively  independent  of  the 

variations  in  the  linear  M-point  for  K  and  K  greater  than  3.   Since  the 

B     coordinate  is  always  the  same,  it  can  intersect  the  left  portion  of 
oave  '  v 
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the  stability  curve  at  only  one  point  and  therefore  the  frequency  of  any 
limit  cycle  resulting  from  a  linear  M-point  located  by  K   ^  3  and  K  *}    3 

is  unaffected  by  variations  of  K  and  K  .  These  same  conclusions  could  be 

J  e      t 

obtained  by  using  X„    and  Figures  VI-4  and  VI-6.   However,  it  must  be 

2max 

pointed  out  that  the  limit  cycle  frequency  and  the  signal  amplitudes 
are  not  necessarily  independent  in  Region  II  of  Figure  VI-3,  thus  the  above 
explanation  is  somewhat  trivial.   In  Regions  I  and  III,  the  frequency  is 
independent  of  amplitudes  and  the  above  result  indicates  that  it  is  also 
essentially  true  in  Region  II  for  K  y    3  and  K   ^»  3. 

3.   Division  of  the  Unstable  Region  into  Subregions. 

The  signals  X.  and  X~  are  not  independent  of  one  another. 
That  is,  except  for  constants  X~  is  the  time  derivative  of  X..   Consequent- 
ly it  is  entirely  reasonable  for  the  system  in  limit  cycle  operation  that 
the  limit  cycle  amplitude  of  one  of  the  signals  is  not  large  enough  to 
allow  the  limit  cycle  amplitude  of  the  other  signal  to  be  great  enough  to 
saturate  the  other  channel.   For  example,  a  limit  cycle  with  small  velocity 
excursions  will  not  drive  the  output  angle  far  enough  before  the  velocity 
changes  sign  to  allow  signal  X.  to  be  greater  than  the  saturation  level  of 
the  error  channel.   This  is  especially  apparent  in  Recorder  Trace  VI- 14 
where  the  limit  cycle  amplitude  of  X~  is  14.3  volts  whereas  the  amplitude  of 
X-  is  only  1.61  volts. 

a.   The  explanation  of  dividing  line  No.  1  is  this.   Record- 
er Traces  VI-12  and  VI-13  are  applicable.   Recorder  Trace  VI-12  shows  the 
limit  cycle  resulting  from  a  linear  M-point  just  to  the  left  of  dividing 

line  No.  1  at  K  =4  and  K  =  0.8.   The  amplitude  of  X,  is  8.1  volts  and 
e         t  1 

the  amplitude  of  X„  is  2.95  volts.   Recorder  Trace  VI-13  shows  the  limit 
cycle  resulting  from  a  linear  M-point  just  to  the  right  of  dividing  line 
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No.  1  at  K  =  4  and  K  =  0.82.   The  amplitude  of  X,  is  unchanged  at  8.1 
e         t  1 

volts  within  the  readable  accuracy  of  the  trace  whereas  the  amplitude  of 

X?  has  increased  to  3.02  volts.   In  this  case  the  amplitude  of  the  velocity 

and  output  angle  probably  did  not  change  enough  to  be  readable  had  these 

signals  been  recorded  because  Qa^  =  2.95/0.8  =  3.69  volts  and  Q>|2  =  3.02/ 

0.82  =  3.69  volts.   However,  the  increase  in  K  was  enough  to  cause  X„ 

'  t         °  2max 

to  be  greater  than  3  volts. 

b.   Dividing  line  No.  2  is  due  to  an  entirely  different  pheno- 
menon since  from  Figure  VI-7  for  K  ^  2.4  the  amplitude  of  X-  is  well  below 
saturation  at  about  1.5  or  1.6  volts.   Also,  since  the  maximum  of  the  stabil- 
ity curve  is  at  K  =  2.45  and  the  amplitude  of  X.,  is  about  half  the  satura- 
tion voltage  at  K  =  2.38,  one  would  expect  something  quite  drastic  to  occur 

as  K  is  increased  above  2.45  because  the  error  channel  must  saturate  for 
e 

larger  values  of  K  .   The  phenomenon  which  occurs  can  be  explained  as  fol- 
e 

lows.   For  a  linear  M-point  just  below  No.  2  dividing  line,  Recorder  Trace 

VI-14,  only  the  tach  channel  saturates  and  the  limit  cycle  average  M-point 

lies  on  the  right  side  of  the  stability  curve.  A  quick  calculation  using 

the  amplitude  X~  and  equation  Vl-d  shows  the  minimum  instantaneous  M-point 

of  the  limit  cycle  is  at  B,  .   =41.5  and  K  =2.38  where  as  an  unstable 

lmm  e 

tach  channel  limit  cycle  exists  at  B,  =41.0  for  K  =2.38.   Thus  the  mini- 

1  e 

mum  instantaneous  M-point  never  crosses  into  the  left  unstable  region  and  a 

dynamic  equilibrium  condition  is  established  about  an  average  M-point  on 

the  right  side  of  the  stability  curve  at  the  K  =2.38  line.   As  K  is  in- 

e  e 

creased  to  place  the  linear  M-point  above  No.  2  dividing  line,  two  effects 

add  together  to  cause  the  response  and  limit  cycle  of  Recorder  Trace  VI-15. 

These  two  effects  are  sketched  in  Figure  VII-5  along  with  B,  ,   of  Recorder 

lmin 

Trace  VI-14.   The  first  effect  is  that  the  increase  in  K  increases  the  re- 

e 

lative  distance  the  average  M-point  must  travel  to  reach  the  stability  curve 
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from   Aa  to  /±a\     As  indicated  previously,  this  causes  the  limit  cycle 
amplitude  of  X«  to  increase  slightly.   Secondly,  the  increase  in  K  causes 
the  graphical  distance  between  a  stable  tach  channel  limit  cycle  and  an 
unstable  tach  channel  limit  cycle  to  decrease  from   A  b  to    Ab'. 
This  change  in    £i  b  is  easily  calculated  from  Equation  Vl-e.   For  K  = 


2.4,  Equation  Vl-e  gives 


co*4-/4o)ni--4  8=o 


from  which  %  , 

LOn       =   <+ 

1 
The  smaller  of  the  values  is  the  desired  one.   Substituting  for  60*  in 


Equation  VI- f  gives 


BL  =  42. 


For  K  =2.38,  Equation  Vl-e  becomes 

eoll4-l4o»n1-47.t=o 

or  COw'2*=  &'  ^'^ 

Substituting  the  smaller  value  in  Equation  VI- f  gives 

Bx  =  40.705 
Since  a  parabola  is  symetric  about  its  axis  the  change  in   A  b  is  then 

Change  in   Ab  =  2(42  -  40.705)  =  2.59 
It  can  be  seen  from  the  figure  that  the  instantaneous  minimum  M-point  for 
K  =2.4  would  lie  in  the  unstable  region  even  if  the  amplitude  of  X^  did 
not  change  by  the  first  effect.   Consequently,  the  small  increase  in  K 
causes  the  instantaneous  minimum  M-point  to  move  past  the  left  portion  of 
the  stability  curve  into  an  unstable  region.   This  in  turn  causes  the  dynamic 
equilibrium  that  would  be  attained  rJuiut  at  point  1  in  Figure  VII-5  to  be 
upset  and  the  signals  increase  in  amplitude.   As  the  signal  X  increases 
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the  average  M-polnt  moves  through  the  stable  region  along  a  constant  Kg 
line.  Eventually  the  error  channel  saturates  allowing  the  system  to 
reach  a  new  limit  cycle  described  by  an  average  M-point  on  the  stability 
curve  somewhere  to  the  left  and  below  point  2  in  Figure  VII-5. 
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Figure  VII-5  &\   ^  ^t 

Sketch  of  Phenomenon  Occurring  at  No.  2  Dividing  Line 
Another  indication  that  the  above  explanation  correctly  describes  the 
phenomenon  is  taken  from  Figure  VI-4.   It  was  noted  that  the  constant  K 
lines  are  linear  with  varying  K  for  M-points  below  No.  2  dividing  line. 
The  line  for  K  =2.4  was  extrapolated  to  K  =  8  which  should  indicate 
X_    if  the  limit  cycle  occurs  on  the  right  side  of  the  stability  curve. 
If  this  X?    is  used  to  obtain  a  B.  .   from  equation  VT-d,  one  gets 

(   3     ) 


B,  .   =8+160 
lmin 


=  35.8 


2  max 


If  this  B-  .   is  now  plotted  on  the  coefficient  plane,  it  falls  just  to 
the  left  of  the  left  side  of  the  stability  curve  which  has  a  B.  coordin- 
ate of  42  for  K  ■  2.4.  Although  this  B^  .   is  further  to  the  left  than 
would  be  expected  it  could  be  considerably  in  error  because  the  extrapolated 
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point  is  far  from  the  last  experimental  point  and  the  lower  portion  of 

the  K  =  2.4  line  in  Figure  VI-4  is  based  on  only  four  experimental  points. 

The  above  explanation  either  greatly  conflicts  with  certain  concepts 
already  used  or  the  author  has  failed  to  see  the  points  that  would  resolve 
the  conflicts.   For  instance,  the  concept  that  the  average  M-point  over  a 
cycle  or  half  cycle  indicates  the  stability  of  the  system  over  the  period 
or  half  period  which  the  average  M-point  is  valid.  The  average  M-point  has 
to  move  through  the  stable  region  to  reach  the  left  side  of  the  stability 
curve  since  the  linear  value  of  K  is  the  upper  vertical  limit  any  M-point 
can  attain.   Average  M-point  motion  through  the  stable  region  indicates 
that  the  signals  should  decrease  slightly  in  amplitude  before  increasing  to 
their  final  limit  cycle  values.   But  all  of  the  Recorder  Traces  exhibiting 
the  same  response  as  Trace  VI-15  showed  increasing  amplitudes  until  the 
final  limit  cycle  amplitude  was  reached.   None  of  the  traces  showed  even  a 
slight  decrease  from  one  peak  to  the  next.   One  explanation  that  might 
resolve  the  conflict  is  that  the  fictitous  average  M-point  jumps  from  the 
right  side  to  the  left  side  of  the  stability  curve.   This  is  not  very  reason- 
able either  since  the  average  M-point  must  in  some  way  be  related  to  the  in- 
stantaneous amplitude  of  signal  X~  and  this  does  not  change  appreciably  from 

peak  to  peak.   Also  the  B,    curve  of  Figure  VI- 10  indicates  a  smooth  trans- 
r  lave 

ition  across  the  stable  region,  although  the  use  of  describing  function  when 
the  system  is  not  in  a  limit  cycle  is  questionable. 

c.   Dividing  lines  independent  of  the  magnitude  of  the  satura- 
tion voltages. 

In  Region  I  just  to  the  left  of  dividing  line  No.  1, 
only  the  error  channel  saturates.   For  linear  M-points  in  this  region  the 
frequency  of  the  limit  cycle  is  completely  determined  by  the  stability  curve 
which  is  independent  of  the  saturation  voltage.   For  linear  M-points  in 
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Region  I,  the  M-point  locus  is  always  the  vertical  line  which  is  the  K 
coordinate  of  the  linear  M-point  which  again  is  independent  of  the  satura- 
tion voltages.   Thus  the  limit  cycle  is  completely  determined  by  the  inter- 
section of  a  curve  which  is  dependent  on  60^and  a  straight  line  which  is 
only  dependent  on  the  linear  M-point.   From  the  intersection,  the  fre- 
quency of  the  limit  cycle  and  a  quantity  called  N..    are  determined.   Only 

when  N,    is  related  to  X,    does  the  saturation  voltage  play  a  part, 
lave  lmax 

Thus  the  relative  distance  the  average  M-point  must  move  to  reach  the 
stability  curve  is  unaltered  by  E   as  long  as  the  linear  M-point  remains 

S  i- 

unchanged  and  in  Region  I.   As  shown  in  the  statement  of  results  and  in 

the  explanation  of  No.  1  dividing  line,  the  dividing  line  occurs  when  X~ 

equals  E  «  and  K  is  increased  from  low  values.   Whatever  the  exact  re- 
n      s2      t 

lationship  between  X,    and  Xn    in  the  limit  cycle  is,  it  should  remain 

lmax      2max 

unaltered  as  long  as  E  .  =  E  „.   Consequently,  No.  1  dividing  line  should 
be  independent  of  the  saturation  voltages  as  long  as  they  remain  equal. 

The  same  arguments  can  be  applied  in  explaining  why  dividing  line  No. 
2  is  apparently  independent  of  the  saturation  voltages.   In  explaining  the 
phenomenon  that  occurs  at  dividing  line  No.  2  it  was  stated  that  in  the  pro- 
cess of  building  up  toward  a  stable  tach  channel  limit  cycle  on  the  right 
side  of  the  stability  curve,  the  amplitude  of  X~  became  large  enough  with 
increasing  K  to  drive  the  instantaneous  minimum  M-point  past  the  left  por- 
tion of  the  stability  curve.   The  pertinent  relative  distances  in  this  ex- 
planation are  independent  of  E  „,.   that  is  distances   A  a  +  &  b  and  Aa 

S/J 

in  Figure  VII-5.   N-    as  a  number  is  completely  determined  by  the  dis- 
tance  ^a.   N0  .   is  defined  as  E  0/X0    and  a  specific  N„  .   is  defined 
2mm  s2   2max        r        2min 

by  an  infinite  number  of  combinations  of  E  -  and  X„    .   It  is  noted  that 

s2      2max 

N?    is  defined  exactly  the  same  as  the  quantity  R  in  the  describing  func- 
tion of  Equation  I-A-3.1  in  Chapter  I.   Thus  once  N~    is  determined  for 
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a  limit  cycle  the  N~  .   or  R  can  be  uniquely  determined  by  the  describing 

2mm 

function.   The  dividing  line  is  then  determined  when  N„  .   as  found  in  the 

°  2min 

above  manner  is  equal  to  the  N„  .   required  to  place  the  instantaneous 

2min 

minimum  M-point  on  the  left  side  of  the  stability  curve.   Since  both  N0  .  '  s 
r  J  2min 

are  independent  of  E  «  and  only  on  the  characteristics  of  the  linear  sys- 
tem. No.  2  dividing  line  must  be  independent  of  E  .. 

PROPOSED  METHODS  OF  PREDICTING  THE  RESULTS  OF  CHAPTER  VI 
1.   Linear  M-point  Inside  Stable  Region. 

a.  Small  step  inputs  which  do  not  initially  saturate  the  error 
channel. 

In  order  to  predict  the  response  of  the  nonlinear  system 
under  these  conditions,  the  initial  peak  of   ^^   for  the  linear  system 
must  be  found.  This  can  always  be  found  by  factoring  the  total  characteris- 
tic equation  with  N.and  N~  equal  to  1.0  and  then  solving  the  linear  differ- 
ential equation  for  a  step  of  E  .  or  less.   Certain  linear  system  predic- 
tion aids  may  be  available  such  that  the  actual  solution  is  unnecessary. 
Then  calculate  X„    from  » 

2max  x       =  k  en 

«max    t   o  max 

If  X„    is  greater  than  E  0  then  calculate  B,  .   from  N„  .  .   If  B,  . 
2max    °  s2  lmin       2min       lmin 

plots  to  the  left  of  the  left  side  of  the  stability  or  on  it,  a  limit  cycle 

could  result  and  the  linear  system  should  be  redesigned.   If  B,  .   plots  to 

lmin 

the  right  of  the  left  side  of  the  stability  curve  then  the  response  of  the 
nonlinear  system  is  essentially  that  of  the  linear  system. 

b.  A  large  step  input  such  that  the  error  channel  initially  satur- 
ates. 

The  first  velocity  peak  again  must  be  found.   This  could  easily 
be  done  by  using  one  of  the  shapes  of  Figure  VII-4  as  applied  to  the  plant 
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if  the  time  t  could  be  determined.   However,  more  study  of  the  problem  is 
required  to  determine  reasonable  values  of  t. .  Since  time  t.  is  indefinite 
an  alternative  is  to  solve  the  differential  equation  by  a  two  part  piece- 
wise  linear  solution.   In  this  method  assume  a  step  input  of  magnitude 

E  ,  is  applied  to  the  plant  with  the  linear  tach  channel  feedback  included, 
si  r 

The  differential  equation  obtained  applies  up  to  time  t  Qf  v  hen  X-  =  E  «. 

This  differential  equation  must  be  solved  for   &o      and  its  derivatives 

up  to  the  point  of  the  previous  criteria.   The  values  of  &q    and  its 

derivatives  at  this  time  are  then  used  as  the  initial  conditions  in  the 

second  differential  equation  which  is  obtained  from  the  plant  alone.   The 

plant  driving  voltage  is  zero  when  the  second  differential  equation  applies. 

With  the  first  velocity  peak  known  B,  .   can  be  calculated  and  plotted  on 

1mm 

the  coefficient  plane.   If  B,  .   plots  to  the  left  of  the  left  side  of  the 

lmin 

stability  curve  a  limit  cycle  may  be  possible.   If  B.  .   plots  in  the 
stable  region  the  response  will  be  similar  to  those  of  Recorder  Traces  VI-2, 
VI-4,  and  VI-6. 

A  second  approximate  method  to  find  the  first  velocity  peak  is  now 
proposed.   The  author  has  not  had  time  to  think  through  the  method  but  it 
does  show  possibilities.   The  method  basically  defines  a  new  equivalent 
instantaneous  transfer  function  which  should  be  applicable,  approximately 
at  least,  through  the  first  velocity  peak.   First  the  lines  of  constant  ^ 
with  (X)y\     as  a  Parameter  must  be  calculated  from  Mitrovic's  equations  of 
Appendix  A  and  then  plotted  on  the  coefficient  plane.   Knowing  the  step  in- 
put magnitude,  the  initial  degree  of  error  channel  saturation  can  be  cal- 
culated from 

I.     ESl 


li    Ei 
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Knowing  N. .  the  Initial  M-point  which  governs  the  system  can  be  calculated 
by  equation  VI-c,  the  other  coordinate  is  the  K  of  the  linear  M-point 
since  N„  =  1.0  initially.   Plotting  this  M-point,  the  &  and  (X>^   which 
initially  govern  the  response  can  be  found  from  the  constant  (/.       lines. 
This  ^  and   (0^  are  for  the  dominant  closed  loop  roots  only.   To  find 
the  remaining  closed  loop  roots  the  instantaneous  characteristic  equation, 

s  +  7s3  +  14s2  +  (8  +  KK  )s  +  KK  N.  .  =0 

t       e  li 


and  the  quadratic   factor, 


s2  +   2      ^Ujfl         s+      UD^ 


must  be  formed  and  the  first  divided  by  the  second.   Thus  all  of  the  roots 
of  an  equivalent  instantaneous  closed  loop  transfer  function  are  known. 
The  gain  of  the  equivalent  closed  loop  transfer  function  is  just  KK  N.. 
This  equivalent  transfer  function  is  shown  in  Figure  VII-6.   Two  assumptions 
must  now  be  made  in  order  to  use  the  initial  equivalent  transfer  function. 
These  are : 

1.  The  nonlinear  error  channel  gain  N   cannot  change  apprecia- 
bly until  the  first  velocity  peak  is  reached. 

2.  The  nonlinear  tach  channel  gain  N~  cannot  change  appreciably. 
The  first  assumption  is  shown  to  be  approximately  true  in  the  recorder 
traces.   The  second  assumption  is  not  true  but  must  be  made  in  order  to 
proceed.   The  primary  question  now  becomes  what  value  of  E.1  to  use  with 
the  equivalent  transfer  function  of  Figure  VII-5.   The  author  has  no  answer. 

A  guess  might  be 

E.1  =  E.  — K  E  . 
l     x         e  si 

since  the  initial  error  channel  saturation  has  already  been  taken  care  of  by 

formulating  the  transfer  function. 
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Figure  VII-6 

Initial  Instantaneous  Equivalent  Closed  Loop  Transfer  Function 
for  Stable  Linear  M-point  with  Large  Input  Steps  Applied. 

The  reason  for  including  the  above  method  was  to  illustrate  an  idea 
which  may  or  may  not  be  applicable  to  the  two  nonlinearity  problem.   Only 
further  study  and  development  will  show  if  the  idea  has  merit.   Certainly 
the  instantaneous  equivalent  closed  loop  transfer  function  is  easily  ob- 
tained from  Mitrovic's  method. 

c.   The  author  was  unable  to  predict  when  the  system  could 
be  driven  into  a  limit  cycle  by  repeated  disturbances.   However,  it  was 
noted  that  for  any  reasonable  damping  the  system  should  remain  stable.   One 

trace  was  taken  for  K  =  2  and  K  =1.5,  which  is  still  quite  close  to  the 

e  t      '  ^ 

left  side  of  the  stability  curve.  The  ^   at  this  point  was  roughly 
^   =0.1.   The  disturbances  applied  were  considered  quite  drastic.   They 
consisted  of  6  to  10  disturbances  alternating  in  sign  on  succeeding  half 
cycles.  The  tach  channel  saturated  to  quite  a  high  degree  but  the  system 
did  remain  stable. 

2.   Prediction  of  Dividing  Lines. 
a.   No.  1  Dividing  Line. 

The  criterion  governing  No.  1  dividing  line  stems  from  the 
relationship  that  must  be  true  between  signals  X  and  X~  under  steady  state 
limit  cycle  operation.   This  relationship  is 
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The  only  quantity  available  from  Mitrovic's  or  the  root  locus  methods 
which  is  related  to  the  input  signal  to  the  nonlinearity  is  N   .   Without 
knowing  the  waveshapes  of  X1  and  X9,  the  above  equation  is  practically  use- 
less for  prediction.   If  sinusoidal  waveshapes  can  be  assumed,  the  above 
equation  can  be  stated  in  a  useful  form.   That  is,  on  an  amplitude  basis 

ift- 

X.   =  ~  UJ  X. 
2m    Ke    lm 

For  sinusoids,  X1   and  X9  occur  a  quarter  cycle  apart  but  this  is  im- 
material in  the  dividing  line  prediction  because  of  the  manner  in  which 
the  dividing  lines  are  defined.   From  the  amplitude  relationship,  the 
criterion  for  No.  1  dividing  line  is  easily  formed  as 

—■  60  X.   =  E  9 
Ke    lm    s2 

Consequently,  it  is  apparent  that  an  accurate  relationship  between  X.   and 

N,     is  needed.   If  the  relationship  used  in  Part  I,  N,    =  E  ,/X,  ,  is 
lave  lave    si   lm* 

used,  the  amplitudes  obtained  differ  greatly  from  the  observed  amplitudes. 

However,  using  this  relationship  in  the  criterion,  the  dividing  line  is 

predicted  to  be  a  vertical  line  starting  from  £q  =  2.0  on  the  stability 

curve,  although  the  use  of  iO     in  the  criterion  then  becomes  questionable. 

In  a  search  for  a  more  accurate  relationship,  the  equivalency  of  N    and 

ave 

G  in  Part  I  was  noted.  Also  the  use  of  the  describing  function  is  indicat- 
ed by  the  fact  that  sinusoidal  signals  had  to  be  assumed  to  put  the  criter- 
ion into  a  useable  form.  Thus  the  criterion  for  No.  1  dividing  line  can  be 

written  as 

Kt      -1 

*  ^  GDI  <Xlm>  "  Es2 

where  G„, (X.  )  =  N.    and  the  symbol  G^,   (X-  )  is  the  inverse  function  of 
Dl   lm     lave  J  Dl     lm 
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of  G  , (X,  )  or  just  X,  (G  ).      To  find  the  point  on  the  stability  curve 
Dl   lnr     J     lmv  Df 

where  the  dividing  line  starts  the  same  criterion  is  used  but  written  in  a 
different  way,  that  is 

Kt  CO  Esl 

K      "  s2 

e 

E   replaces  X.   by  the  following  reasoning.   For  linear  M-points  very 

close  to  the  stability  curve,  the  degree  of  saturation  required  to  move 

the  average  M-point  to  the  stability  curve  is  very  small.   Thus  X.     is 

essentially  equal  to  E  ..   It  is  then  just  a  matter  of  trial  and  error  to 

find  the  point  at  which  K  ,  K  ,  and  iO     satisfy  the  criterion.   It  should 

e   t 

be  noted  that  for  this  particular  point  on  the  dividing  line  waveshapes  are 
unimportant  if  LO      is  considered  the  average  limit  cycle  for  frequency 
obtained  from  the  average  period. 

In  applying  the  criterion  to  calculate  the  dividing  line,  the  follow- 
ing procedure  was  used.   A  linear  M-point  was  chosen  and  error  channel  satura- 
tion only  was  assumed.   From  the  intersection  of  the  K  coordinate  of  the 
J  t 

linear  M-point  and  the  stability  curve  the  (JQ   and  N,    of  the  resulting 
r  lave  & 

limit  cycle  were  read.  A  table  of  X.   and  the  resulting  G  was  calculated 

by  a  digital  computer  program.   Entering  the  table  with  N.    ,  the  limit 

cycle  amplitude  X,   was  obtained.   Using   60  and  X,   along  with  the  linear 
1m  lm 

M-point  coordinates  in  the  criterion,  a  value  of  X_  was  obtained.   If  this 

zm 

value  of  X-  was  less  than  E  „,  K  was  increased  and  a  new  X„  calculated 
2m  s2'   t  2m 

using  the  same  procedure.   This  procedure  was  continued  until  E  0  was 

s2 

bracketed  by  X„  .   The  procedure  rapidly  converges  to  give  the  K   for  the 
2m  t 

chosen  K  which  results  in  X„   =  E  „.   Table  VI-1  gives  a  comparison  of  the 
e  2m    s2  or 

calculated  and  experimental  linear  M-points  which  define  dividing  line  No.  1. 
The  results  of  Table  VI-1  were  considered  quite  good  since  the  analog  computer 
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TABLE  VI- 1 


K 

Calculated 

Experimental 

6 

Kt 

Kt 

2,1 

0.883 

0.92 

2.3 

0.855 

0.89 

2,5 

0.838 

0.85 

3,0 

0.818 

0.82 

4.0 

0.79 

0.81 

5.0 

0.783 

0.80 

6,0 

0.778 

0.785 

results  are  subject  to  3-5%  error. 

The  criterion  can  be  put  into  a  form  in  which  once  choosing  a  linear 
M-point  K  ,  the  linear  M-point  K  of  the  dividing  line  can  be  found  direct- 
ly but  the  algebra  is  quite  involved  due  to  getting  X   (G  )  from  equation 
I-A-3.1  of  Part  I.   The  inverse  sine  and  square  root  terms  of  Equation  I-A- 
3.1  can  be  expanded  in  an  infinite  series  form.   If  the  two  infinite  series 

are  added,  reversion  of  the  series  sum  can  be  accomplished  to  give  X.   as 

lm 

a  function  of  G  .   At  the  intersection  of  the  linear  M-point  K  coordinate 

and  the  stability  curve,  N~  assumed  1.0,  Equation  VI-c  equals  Equation  Vl-e 

and  Equation  Vl-d  equals  Vl-f.   From  these  equations  expressions  for 

and  N   s  which  equals  G^.  can  be  obtained  in  terms  of  K  and  K  .   Substi- 
aves        n      D*  e      t 

tution  for  the  various  quantities  results  in  a  statement  of  the  criterion 

in  terms  of  K  and  K  of  the  linear  M-point  only.   The  expression  is  very 

complicated  and  includes  a  power  series  whose  argument  contains  both  K  and 

e 


V 
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b.   Prediction  of  No.  2  Dividing  Line. 

From  Figure  VI-5  and  the  explanation  of  the  phenomenon 

occurring  at  the  dividing  line,  the  criterion  governing  the  dividing  line 
can  be  stated  as  follows: 


"For  the  linear  M-point  on  the  dividing  line,  B.  .  ,  as 

calculated  for  a  stable  tach  channel  limit  cycle  on  the 
right  side  of  the  stability  curve  must  equal  the  B1  co- 
ordinate of  the  stability  curve  on  the  left  side  where 

the  linear  M-point  K  coordinate  intersects." 
r      e 

Calculation  of  the  dividing  line  is  again  a  trial  and  error  procedure.   One 

such  procedure  using  B,  coordinates  is  as  follows.   Choose  a  K  and  a  K  . 

1  t        e 

Reading  the  B.  coordinate  of  the  K  line  intersections  from  the  graph  was 
found  to  be  too  inaccurate  to  give  adequate  results.   Thus  the  B  co- 
ordinates of  the  intersections  of  the  K  line  with  the  stability  curve  had 

e 

Z 
to  be  calculated  from  Equations  Vl-e  and  Vl-f.   From  Equation  Vl-e,  0d*\ 


is  found  by  solving. 


a.     4 

14  tO*     -  k>nT   =  B  =  KK 


o     e 


•2, 
Substitution  of  the  two  values  of  60^       into  Equation  Vl-f  gives  the  de- 
sired B.  coordinates.   The  following  symbols  will  be  used: 

B,   =  The  coordinate  of  the  right  intersection  at  which  a 
lr 

stable  tach  channel  limit  cycle  exists, 

B   =  The  coordinate  of  the  left  intersection  at  which  an 

unstable  tach  channel  limit  cycle  exists. 

Using  B,  ,  N0    can  be  calculated  from  Equation  Vl-d.   From  N 

lr   2ave  n  ave 

=  G_  and  using  the  pre-calculated  table,  the  value  of  N0  .  can  be  deter- 
D  Zmm 

mined.   It  is  noted  that  the  definition  of  L  .  ,  N,  .   =  E  „/X0    ,  is 

Zmin   2mm    s2   2max 

exactly  the  definition  of  the  quantity  R  used  in  the  describing  function  in 

Part  I,  Equation  I-A-3.1.   Having  N_  .  ,  B,  .   can  be  calculated  using 

Zmin   1mm 
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Equation  Vl-d.   A  comparison  of  B,  .   and  B,T  can  be  made.   If  B,  .   is 
n  r  lmin      1L  lmin 

greater  than  B.  *  the  linear  M-point  was  below  the  dividing  line  and  K 

JL  Ju  C 

must  be  increased.   If  B,  .   is  less  than  B,,  ,  the  linear  M-point  was 

lmin  1L 

above  the  dividing  line  and  K  must  be  decreased.   Table  VII-2  gives  a 

e 

comparison  of  the  experimental  dividing  line  and  the  calculated  dividing 
line  using  the  procedure  outlined  above. 

TABLE  VII-2 


K 

t 

Calculated  K 
e 

Experimental  K 
e 

3 

2.431 

2.425 

4 

2.429 

2.415 

5 

2.430 

2.405 

6 

2.429 

2.395 

8 

2.427 

2.385 

10 

2.425 

2.335 

The  last  significant  figure  shown  in  the  calculated  K  was  interpolated 

between  the  differences  existing  in  the  comparison  of  B,  .   with  B,„ 

1mm       1L 

for  the  calculations  at  K  =  2.42  and  K  =  2.43.  The  results  are  surpris- 

e  e  r 

ingly  good  considering  the  calculated  K  is  subject  to  slide  rule  accuracy 
and  considering  the  analog  results  are  inherently  subject  to  the  3-5%  ac- 
curacy.  The  worst  percentage  difference  between  calculated  K  and  experi- 
mental K  is  less  than  4%. 

e 
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CHAPTER  VIII 
SOME  PREDICTIONS  OF  RESULTS  FOR  DIFFERENT  SYSTEMS 
WITH  TWO  GAIN  VARIABLE  NONLINEARITIES 

An  attempt  is  made  in  this  chapter  to  extend  certain  of  the  ideas 
developed  in  study  of  the  system  of  Figure  VI  with  equal  saturation  voltages. 
However,  it  must  be  pointed  out  that  the  extension  of  the  ideas  is  mostly 
speculation  with  no  specific  experimental  results  to  back  up  the  ideas. 

An  attempt  will  first  be  made  to  indicate  the  expected  results  when 
different  parameters  of  the  system  of  Figure  VI  are  changed.   In  Chapter 
VII  it  has  already  been  shown  that  both  dividing  lines  are  independent  of 
the  saturation  voltages  as  long  as  they  remain  equal.   It  is  then  reasonable 
to  expect  the  same  types  of  responses  for  linear  M-points  in  the  stable  re- 
gion as  the  equal  saturation  voltages  are  changed.   The  argument  applied  is 
that  the  relative  distances  between  a  given  linear  M-point  and  the  stability 
curve  remain  unaffected  by  the  saturation  voltages,  thus  the  degree  of  satura- 
tion required  is  unchanged.   Only  the  amplitudes  of  the  signals  are  changed 
to  cause  the  required  degrees  of  saturation. 

For  the  same  linear  system  with  E  „  V  E  ,  certain  results  seem  evi- 

s2  <    si 

dent.   When  the  linear  M-point  lies  in  the  stable  region  and  for  small  in- 
put steps  in  which  the  error  channel  does  not  initially  saturate,  it  should 
become  increasingly  difficult  to  find  a  linear  M-point  at  which  the  tach 

channel  saturates  as  the  difference  between  E  0  and  E  ,  is  increased.   A 

s2      si 

point  should  be  reached  at  which  E  ~  exceeds  E  .  by  such  a  large  amount 

that  the  response  remains  linear  throughout.   By  the  same  reasoning  that 

the  tach  channel  saturates  less  and  less  as  the  difference  between  E  „  and 

s2 

E  ,  is  increased,  it  should  be  more  difficult  to  drive  the  system  into  a 
si  '  J 
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limit  cycle  by  repeated  disturbances.   For  large  signal  inputs  such  that 

the  error  channel  initially  saturates,  it  should  be  found  that  the  entire 

response  approaches  error  channel  saturation  only  as  the  difference  between 

E  „  and  E  ,  is  increased. 
s2      si 

For  E  _  ^  E  . 9  the  dividing  lines  in  the  unstable  regions  should  be 

S  £-  5  X 

predictable  by  the  criterion  as  stated  in  Chapter  VII.   For  dividing  line 
No.  1,  the  criterion  can  be  rewritten  as 

Kt  \XJ  =  Es2 

K  X. 

e  lm 

Assuming  only  error  channel  saturation,  X.   is  independent  of  E  »  since 
the  vertical  distance  from  linear  M-point  to  the  stability  curve  is  indepen- 
dent of  E  .•   Thus  for  a  given  K  ,  one  would  expect  the  dividing  line  to 
shift  to  the  right  where  K  and  OJ    are  larger  as  E  .  is  increased  and 

E  ,  is  kept  constant.   A  minimum  of  calculations  show  that  for  E  ,  =  3  and 
si  si 

E  j   =  5,  the  dividing  line  starts  at  (j)n  =  2.42  on  the  stability  curve  by 
applying  the  criteria  in  the  form 

Kt  U)        E  .  =  E 

— ac-    si    s2 

K 
e 

Another  point  for  the  same  saturation  voltages  was  calculated  for  K  =  5 

and  was  found  to  be  about  K  =1.4.   The  criterion  in  the  form  useful  at  the 

t 

stability  curve,  indicates  that  the  difference  of  E  „  and  E   is  unimportant 

to  the  placement  of  the  dividing  line.   It  does,  however,  indicate  that  the 

ratio  of  E  _  to  E  ,  is  the  quantity  which  determines  the  placement  of  divid- 
sz     si 

line  No.  1.   It  has  already  been  shown  in  Chapter  VII  that  for  equal  satura- 
tion voltages,  a  ratio  of  1.0,  the  dividing  line  is  independent  of  their 
magnitude. 

This  same  reasoning  also  applies  to  the  case  where  E  ~  ^  E  .  only 

the  dividing  line  would  shift  to  the  left  of  the  position  for  E  ,  =  E  «. 
°  si    s2 
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Thus  for  any  given  linear  plant  and  any  saturation  voltages,  dividing 
line  No.  1  should  be  predictable.  The  only  restriction  is  that  the  plant 
and  the  feedback  paths  have  the  same  configuration  as  Figure  VI.   Postulat- 
ing even  farther,  if  the  linear  plant  can  be  normalized,  which  has  already 
been  done  for  a  third  order  plant,   then  a  master  set  of  curves  should  be 
possible  defining  dividing  line  No.  1  in  terms  of  the  ratio  of  E  .  to  E  ., 
For  the  system  of  Figure  VI,  described  by  the  coefficient  plane  of  Figure 

VI-3,  it  is  unclear  what  happens  when  the  ratio  of  E  „  to  E  ,  becomes  large 

r  s2     si  ° 

enough  such  that  the  starting  point  of  dividing  line  No.  1  moves  to  the  right 

of  the  stability  curve  maximum.   At  the  maximum  for  E  .  =  3,  the  criteria 

gives  E  «  =  6.64  or  a  ratio  of  2.213.   For  this  ratio  and  K  =  5,  another 
s2  e 

point  on  the  dividing  line  is  at  K  =  1.72.   The  criterion  applied  at  the 

stability  curve  keeps  giving  reasonable  values  of  E  „  for  points  to  the 

right  of  the  maximum.   For  example,  for  E  .  =  3  and  LOf^  -   3,  the  required 

E  •  is  11  volts. 
s2 

It  was  shown  in  Chapter  VII  that  dividing  line  No.  2  was  independent 

of  the  saturation  voltages.   Although  it  was  not  stated  explicitly  the 

criterion  for  the  dividing  line  is  only  dependent  on  N„    and  N„  .   as  de- 

2ave      2mm 

grees  of  saturation.  Since  the  dividing  line  is  independent  of  the  error 
channel  and  it  has  already  been  shown  to  be  independent  of  E  »  in  Chapter 
VII,  it  is  therefore  the  same  for  any  combination  of  saturation  voltages. 

Since  No.  2  dividing  line  depends  only  on  the  degree  of  saturation  in 
the  tach  channel,,  the  error  channel  may  saturate  as  E  9  becomes  greater  than 

1 

Thaler,  G.  J.  and  Brown,  R.  G.,  Analysis  and  Design  of  Feedback  Control 

Systems,  McGraw-Hills  2nd  Ed.,  1960:  380-383. 
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E  ..   It  was  implied  for  equal  saturation  that  No.  2  dividing  line 

S  1 

separated  linear  M-points  with  tach  channel  saturation  from  linear  M-points 
with  both  channels  saturating.   A  better  description  of  the  dividing  line 
No.  2  is  that  it  separates  linear  M-points  for  which  limit  cycles  occur  on 
the  right  side  of  the  stability  curve  with  the  tach  saturation  producing 
the  stable  limit  cycle  from  linear  M-points  for  which  limit  cycles  occur 
on  the  left  side  of  the  stability  curve  with  the  error  channel  producing 
the  stable  limit  cycle.   Thus  for  linear  M-points  below  No.  2  dividing  line 
it  is  entirely  possible  to  get  both  channels  saturating  for  E  „  \     E  1 
although  the  error  channel  saturation  is  incidental  to  the  stable  limit 
cycle.   Consequently,  a  third  dividing  line  must  be  postulated.   The  des- 
cription of  the  third  dividing  line  would  be  that  it  separates  linear  M- 
points  for  which  only  the  tach  channel  saturates  from  linear  M-points  for 
which  both  channels  saturate.   For  a  linear  M-point  in  the  region  between 
No.  2  dividing  line  and  No.  3  dividing  line  the  limit  cycle  will  occur  on 
the  right  side  of  the  stability  curve  with  the  stable  limit  cycle  produced 
solely  by  the  tach  channel.   The  criterion  for  No.  3  dividing  line  has  to 
be  the  same  as  that  for  No.  1  dividing  line  only  applied  in  reverse.   In 
general  terms  the  limit  cycle  relationship  between  X  and  X«  is 

Assuming  sinusoidal  waveshapes  of  X. (t)  and  X9(t),  then 


K«2 


A  n  rf\   —  A» 


Then  forming  the  dividing  line  criterion  yields, 

At  the  stability  curve  the  criterion  would  simply  be 

-&£-  P   -F 
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Since  dividing  line  No.  1  and  dividing  line  No.  3  are  obtained  by  differ- 
ent ways  of  applying  the  same  criterion,  it  is  reasonable  to  say  they  are 

the  same  dividing  line  at  least  for  an  E  0/E  ,  ratio  large  enough  for  the 

°  s2   si 

right  hand  portion  to  be  below  No.  2  dividing  line.   With  this  assumed  true} 

for  a  few  points  on  the  right  side  stability  curve  the  E  2  voltage  required 

to  lust  make  X,    =  E  ,  =  3  were  calculated  by  the  criterion.   Then  a  few 
J  Imax    si 

points  on  the  dividing  line  at  three  of  these  voltages  were  calculated  at 

K  =  8.   A  very  rough  sketch  of  what  the  complete  set  of  curves  might  look 

like  is  given  in  Figure  VIII-1.   The  points  actually  calculated  are  circled. 

The  dividing  lines  are  labeled  in  terms  of  the  ratio  E  0/E  , ,  but  were  cal- 

s2   si 

culated  assuming  E  ,  =  3  and  E  „  necessary  to  satisfy  the  criterion  at 
°   si  s2 

the  stability  curve.  The  dividing  lines  with  no  calculated  points  other 
than  at  the  stability  curve  are  estimated  from  the  shape  of  the  ones  with 
calculated  points.   The  shape  of  the  dividing  lines  from  the  left  side  of 
the  stability  curve  were  estimated  by  the  shape  for  E  „/E   =1.0  although 

the  value  of  E  „,/E  ,  was  calculated  by  the  criterion  at  the  stability  curve. 

s2   si  J 

If  the  above  postulation  of  dividing  lines  proves  correct,  then  for  a  divid- 
ing line  such  as  E  ~/E  1  =8.77  the  only  meaningful  portion  of  dividing  line 
No.  2  is  the  portion  to  the  right  of  the  8.77  dividing  line.   The  shape  of 
the  dividing  lines  close  to  stability  curve  has  not  been  determined.   Fig- 
ure VIII-2  shows  three  possibilities  which  may  occur.   If  these  postulated 
dividing  lines  can  be  confirmed  by  analog  computer  results,  then  it  should 
be  interesting  to  see  if  the  dividing  lines  for  E  _/E  £~.   2.233  stop  at 
the  stability  curve.   For  E  ~fz    .  ^  2.33  the  dividing  lines  at  worst  have 
a  discontinuity  at  the  stability  curve.   Thus  it  may  be  reasonable  that  the 
other  dividing  lines  may  enter  the  stable  region.   If  this  is  true  then  they 
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Figure  VIII-2 

Possible  Shapes  of  the  Dividing  Lines  at  the  Stability 

Curve  for  Large  E  0/E  ,  Ratios. 
s2  si 
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should  have  more  than  a  casual  association  with  possible  limit  cycle  opera- 
tion of  systems  whose  M-point  lies  in  the  stable  region.   In  the  small 
amount  of  investigation  done  in  driving  a  stable  linear  system  into  a  limit 
cycle  with  repeated  disturbances,  all  the  points  investigated  were  close 
to  the  upper  part  of  the  stability  curve  and  not  much  lower  than  K  =  2.0 
or  ti){\   =  2.0.   This  may  indicate  that  the  dividing  line  for  E  JY.        =1.0 
may  circle  through  the  stable  region  and  leave  at  dividing  line  No.  2.  The 
author  just  did  not  have  the  time  to  see  if  the  dividing  line  criterion 
could  give  answers  for  M-points  in  the  stable  region.   Certainly  from  the 
concepts  of  degrees  of  saturation  and  relative  distances  it  should  be  easier 
to  drive  the  instantaneous  M-point  past  the  unstable  tach  channel  limit  cycle 
for  E  9/E  1   ^  1.0  on  the  first  velocity  peak. 

When  the  plant  gain  K  is  changed,  not  very  much  should  change  as  far  as 
Figure  VI-3  is  concerned.   The  stability  curve  does  not  change  since  it  is 
independent  of  K,   The  scales  on  the  K  -K  grid  will  vary  inversely  with  K. 
That  is,  if  K  is  doubled,  the  scales  will  be  halved.   The  all  important  re- 
lative distances  between  linear  M-point  and  the  stability  curve  will  not 
change  as  K  increases,  they  will  only  have  different  K  and  K  coordinates. 

Since  the  relative  distances  determine  N   ,  then  No.  1  and  No.  2  dividing 

ave'  & 

lines  are  unchanged  as  far  as  the  B  -B.  coordinate  system  is  concerned.   If 
the  above  extension  of  the  dividing  lines  proves  to  be  valid,  then  they  will 
be  unchanged. 

The  author  can  concieve  of  no  reason  why  most  of  the  concepts  discussed 
in  connection  with  the  plant  of  Figure  VI  can  not  be  applied  with  different 
plants  in  the  block  except  for  second  order  systems.   Mitrovic's  method 
does  not  lend  itself  easily  to  the  analysis  of  second  order  systems.   For  a 
third  order  plant  the  stability  curve  is  a  straight  line  and  thus  dividing 
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line  No.  2  could  not  exist.  Also,  with  a  third  order  system,  the  filter- 
ing may  be  insufficient  such  that  the  signal  X~  is  not  sinusoidal  in  which 
case  the  criterion  for  dividing  line  No.  1  would  give  poor  predictions  since 
it  makes  use  of  the  describing  function. 

With  different  types  of  gain  variable  nonlinearities  in  the  blocks 
the  concept  of  dividing  lines  and  a  knowledge  of  the  type  of  limit  cycles 
produced  by  each  nonlinearity  individually  should  give  a  good  indication  of 
the  actual  responses.   For  example,  assume  two  elements  with  dead  zone  in 
the  nonlinear  blocks  of  Figure  VI.   In  this  case  the  entire  stability  curve 
represents  unstable  limit  cycles  for  the  error  channel  dead  zone.   For  the 
dead  zone  in  the  tach  channel,  the  right  side  of  the  stability  curve  represents 
unstable  limit  cycles  and  the  left  side  represents  stable  limit  cycles.   It 
cannot  be  definitely  concluded  that  a  system  with  a  linear  M-point  in  the 
unstable  region  will  always  go  into  a  stable  limit  cycle.   For  instance, 
consider  a  linear  M-point  well  above  the  stability  curve.   If  the  average 
value  of  signal  X.  over  a  half  cycle  is  initially  large,  then  N-    can 
never  be  small  enough  to  move  the  average  M-point  downward  below  the  maxi- 
mum of  the  stability  curve  where  the  tach  channel  N„    can  act  to  move  the 

2ave 

average  M-point  to  the  left  side  of  the  stability  curve.  At  this  same  linear 
M-point s  it  should  be  possible  to  obtain  a  stable  response  for  small  input 
signals,  that  is,  N.    and  N~    never  get  large  enough  to  move  the  in- 
stantaneous M-point  outside  the  stable  region.   Thus  it  can  be  seen  that 
the  analysis  is  much  more  complicated  than  for  two  saturations  primarily 
because  the  initial  M-point  starts  at  the  point  (0,0)  of  the  K  -K  grid. 
However,  with  more  study,  one  should  be  able  to  relate  certain  M-point 
motions  to  the  position  of  the  linear  M-point  and  thereby  be  able  to 
describe  certain  dividing  lines.   If  dividing  lines  can  be  found,  then 
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their  position  should  be  dependent  only  on  a  ratio  of  the  two  amounts  of 

dead  Eone  because  N    and  N    would  only  be  dependent  on  the  appropriate 

ave      max 

relative  distances  between  linear  M-point  and  stability  curve. 

In  all  of  the  previous  discussion,  a  reliance  on  the  describing  func- 
tion is  strongly  implied.   It  should  be  point  out  at  this  point  that  in 
general  the  describing  function  and  the  requirement  of  sinusoidal  wave- 
shapes are  not  necessary.  The  only  necessary  quantities  are  the  relative 
distances  involved  in  the  coefficient  plane  since  these  distances  can  be 
related  to  various  phenomena  occurring  in  the  experimental  results.   The 
relative  distances  can  be  described  and  measured  entirely  by  the  somewhat 

nebulous  quantities  N    and  N  .   or  N   ,   The  author  has  only  been  able 

ave      mm     max 

to  give  N    and  N  .   meaningful  interpretations  for  prediction  purposes 

°      ave       min        or  r  r   r 

through  the  use  of  the  describing  function.   It  should  be  stressed  again 

that  N    is  just  a  symbol  which  is  used  for  lack  of  a  better  symbol.   Even 

ave 

for  a  sinusoidal  signal,  a  straight  integral  time  average  of  the  instantan- 
eous gain,  N  ,  over  a  cycle  does  not  corrolate  with  the  observed  maximum 
signals.   However,  for  sinusoidal  signals  the  describing  function  does  corro- 
late N    with  maximum  observed  signal.   Consequently,  the  relation  between 
the  amplitude  of  the  input  signal  and  N    is  much  more  complicated  than  a 
time  average  of  the  instantaneous  gain  and  for  non  sinusoidal  waveshapes 
this  relation  must  certainly  be  dependent  on  the  waveshape.   If  a  relation- 
ship between  N    and  the  input  signal  can  be  found  which  is  independent 
of  waveshape,  then  all  of  the  prediction  methods  should  be  valid. 
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CHAPTER  IX 
MISCELLANEOUS  IDEAS  AND  CONCEPTS 

In  studying  the  two  nonlinearity  problems,  the  author  had  a  few 
ideas  which  were  not  particularly  concerned  with  the  problem  or  which  did 
not  give  pertinent  results.   These  ideas  are  considered  to  have  possibil- 
ities but  were  not  developed. 

The  first  idea  to  be  discussed  concerns  a  general  method  for  cal- 
culating the  limit  cycles  in  any  system  with  any  number  of  nonlinearities 
if  any  limit  cycles  exist  in  the  system.   If  no  limit  cycle  existss  then 
the  method  should  give  no  non-trival  solutions.   The  only  restrictions  on 
the  method  are  these: 

1.  Each  nonlinearity  must  be  such  that  it  can  be  described  by 
a  describing  function. 

2.  The  linear  system  must  provide  enough  filtering  between  the 
nonlinearities  such  that  the  describing  function  description  can  be  justi- 
fied. 

The  basis  of  the  method  is  application  of  the  stability  criterion  to 
the  total  nonlinear  system.   In  Chapter  VII  it  was  indicated  that  the  crit- 
erion for  the  root  locus  method  degenerates  into  the  Niquist  criterion  at 
the  limit  of  stability  or  when  s  =  J  (jO  .   In  applying  the  root  locus 
criterion  to  the  total  nonlinear  system,  each  nonlinearity  must  somehow  be 
included  as  a  magnitude  and  angle.   At  the  stability  limit  or  under  limit 
cycle  conditions,  the  describing  function  provides  the  description  of  the 
nonlinearity  as  a  magnitude  and  angle.   In  other  words,  the  describing  func- 
tion can  be  thought  of  as  linearizing  the  nonlinearity  on  the  average  over 
a  cycle.   Thus  the  root  locus  criterion  can  be  formed  for  the  total  system 
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including  the  nonlinear  effects  by  the  magnitude  and  angle  of  their 
describing  function.   Since  the  root  locus  criterion  is  just  a  method 
of  rearranging  the  total  characteristic  equation,  one  may  as  well  start 
with  the  characteristic  equation.   It  then  becomes  convenient  to  express 
the  describing  functions  as  a  real  part  and  an  imaginary  part. 

For  illustration  purposes  consider  the  system  configuration  of  Figure 
VI.   Also  let  the  plant  be  fourth  order,  although  any  order  plant  can  be 
used,  and  let  the  two  nonlinear  elements  be  any  general  nonlinear  elements, 
That  is,  the  nonlinear  elements  can  be  gain  variable  only,  hysteretic,  or 
frequency  sensitive.   The  describing  function  of  the  nonlinearities  can  be 
written  as 

%(\m,W      -   Nlr(Xlm,U)>  +  NU  (Xlm>^  > 
and 

N2(X2ni>tO)      =  N2r(X2m>UJ)   ♦     N2.(X2m>W) 

The   total   characteristic  equation   is   then 

a.s      +  a.s      +  a.s      +    (a,    +   KKN0)s   +   KK  N.    =0 
4  3  2  v    1  t   2'  el 

substitution  of  the  above  describing  functions  for  N.  and  N~ 
At  the  stability  limit  for  the  total  system,  s  =  Jlu  ,  thus 

+  KKeQ*ir(Xnfl,O))4-J,Na(Xim,(0)31  =0 

Thus  two  independent  equations  can  be  obtained  from  the  above  equation 

each  containing  the  three  variables,  £0»  3L  .  and  X0  . 

lm       2m 
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Clearly  a  third  equation  must  be  obtained  and  this  equation  must  come 
from  the  physical  configuration  of  the  system  and  must  relate  X1   to  X2  . 
For  other  than  gain  variable  nonlinearities,  the  phase  difference  must  be 
maintained  between  X-  and  X~.   Thus  for  the  configuration  of  Figure  VI, 
the  relationship  is 

or  for  limit  cycle  operation, 

Substitution  of  this  relationship  into  the  characteristic  equation  will 
give  one  equation  in  two  unknowns  which  can  be  broken  up  into  two  indepen- 
dent equations  in  the  same  two  unknowns.   Of  course  the  nature  of  the  des- 
cribing functions  will  make  the  algebra  very  messy  especially  manipulating 
the  terms  containing  J   into  real  and  imaginary  terms.   Conceivably  this 
can  be  done  and  then  the  single  equation  can  be  broken  up  into  the  two  indep- 
endent equations,  the  solutions  of  which  will  not  be  simple. 

Thus,  in  theory  anyway,  the  method  is  not  restricted  to  the  number  of 
nonlinearities  as  long  as  a  sufficient  number  of  relationships  between  the 
individual  nonlinearity  input  signals  can  be  determined  from  the  physical 
configuration  of  the  system.   The  only  limitation  is  the  justification  of 
the  use  of  the  describing  function  for  each  nonlinearity,  which  in  itself 
will  eliminate  a  vast  number  of  multiple  nonlinearity  system  configurations. 

A  second  idea  came  about  in  trying  to  predict  the  limit  cycle  for  the 
configuration  of  Figure  VI  for  a  linear  M-point  in  Region  II  where  both 
channels  saturate.   The  limit  cycle  must  surely  be  described  by  the  inter- 
section of  an  average  M-point  locus  and  the  stability  curve.   The  problem 
becomes  one  of  predicting  the  average  M-point  locus  from  the  linear  M-point 
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to  the  stability  curve.   It  has  been  demonstrated  that  in  limit  cycle  condi- 
tions the  amplitude  of  the  nonlinearity  input  signal  can  adequately  be  relat 
ed  to  N    through  the  describing  function.   Consequently,  it  was  reasoned 
that  in  a  highly  oscillatory  condition  other  than  the  limit  cycle  condi- 
tion, the  system  can  be  thought  of  as  being  in  a  quasi-limit  cycle  at 
least  over  a  half  cycle  of  oscillation.   Thus,  it  appeared  logical  that 
the  describing  function  could  be  used  to  relate  the  average  M-point  to  the 
maximum  amplitude  over  the  half  cycle  in  question.   The  envelope  of  the 
oscillations  then  becomes  important  since  this  determines  amplitude  of  the 
oscillations.   It  was  then  theorized  that  the  total  envelope  could  be  ap- 
proximated by  a  series  of  exponential  segments  with  each  segment  dependent 
only  on  the  real  part  of  the  average  root  over  the  half  cycle  as  determined 
by  the  average  M-point  of  the  half  cycle. 

For  linear  M-points  in  the  unstable  region,  the  highly  oscillatory 
condition  is  assured  since  for  the  system  studied  a  stable  limit  cycle  al- 
ways results.  Thus  the  following  method  was  evolved.   For  a  chosen  linear 
M-point,  the  resultant  closed  loop  transfer  function  can  be  formed  and  its 
denominator  factored.   The  linear  differential  equation  is  easily  solved 
by  operational  methods  to  give  the  initial  linear  response.   An  impulse 
forcing  function  of  such  a  magnitude  was  assumed  such  that  the  linear  enve- 
lope of  signals  X.  and  X~  essentially  retained  symmetry  about  zero.   The 
amplitude  of  successive  peaks  of  both  signals  X1  and  X9  were  calculated 
until  a  peak  of  one  of  the  signals  was  greater  than  the  saturation  voltage. 
When  a  signal  peak  was  found  that  exceeded  the  saturation  voltage,  the 
average  M-point  was  calculated  for  that  half  cycle  by  using  the  peak  magni- 
tude in  the  describing  function  to  determine  the  N    for  the  saturated 

ave 

channel.  This  new  average  M-point  was  plotted  on  the  coefficient  plane. 
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In  order  to  obtain  the  information  necessary  to  calculate  the  next  exponent- 
ial segment  of  the  envelope,  a  curvilinear  grid  of  constant  d  and  UJn lines 
had  to  be  superimposed  on  the  coefficient  plane.   This  curvilinear  grid  was 
calculated  from  Mitrovic's  mapping  equations  of  Appendix  A  for  negative 
values  of   ^  •  From  the  position  of  the  average  M-point  with  respect  the 

(£  and  LDft    grid,  a  new  average  */,  and  (jOf\     was  read  which  were  assumed  to 
govern  the  response  over  the  next  half  cycle.   Thus  the  envelope  equa- 
tions for  signals  X  and  X~  for  the  next  half  cycle  can  be  written  as 

«t  -  V  e  *"** 

and 

x2  =  Xax-<2 

The  quantities  X, ,  and  X„ .  are  the  values  of  the  envelope  calculated  from 
li      2i 

the  linear  solution, or  the  end  points  of  each  preceding  exponential  seg- 
ment.  The  average  frequency  over  a  half  cycle  changes  as  the  governing 
4£  and  60ft  of  the  half  cycle  change  as 

From  this  relation  a  new  period,  T   ,  can  be  obtained.   Since  the  next 

r     '   ave ' 

peak  of  the  signal  under  consideration  occurs  a  half  cycle  later,  the 

magnitude  of  the  new  envelope  segment  at  the  end  of  the  next  half  cycle 

was  calculated  by  setting  t  =  T   /2  in  the  envelope  equations.   Thus  new 

values  of  X,    and  X0    were  determined  which  were  used  in  the  describ- 
Imax      zmax 

itig  function  to  determine  a  new  M-point  to  be  used  over  the  next  half  cycle. 
Thus  by  successive  calculations  the  average  M-point  locus  was  constructed. 

Several  average  M-point  loci  were  calculated  in  this  manner  but  none 
of  them  ever  crossed  the  stability  curve.   A  typical  locus  calculated  in 
this  manner  is  shown  in  Figure  IX- 1.   As  can  be  seen  the  locus  approaches 
but  never  reaches  the  stability  curve.   It  appears  to  continue  indefinitely 
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past  the  experimental  limit  cycle  while  never  quite  reaching  the  stability 
curve.   The  behavior  of  the  locus  calculated  in  the  manner  described  is 
reasonable  since  the  average  M-point  can  approach  the  stability  cuive  as 
close  as  one  pleases  but  the  z>     to  be  used  in  the  next  half  cycle  is 
always  negative  causing  the  next  peak  to  be  slightly  greater  than  the 
previous  peak.   The  method  as  outlined  does  not  take  into  account  the 
phase  difference  between  the  signals.   One  improvement  would  be  to  take 
the  phase  difference  into  account  possibly  by  using  quarter  cycle  incre- 
ments instead  of  half  cycle  increments. 

An  observation  made  in  the  application  of  the  method  for  the  several 
M-points  tried  is  that  the  manner  of  making  the  transition  from  the  linear 
M-point  to  the  first  average  M-point  does  not  affect  the  locus  appreciably 
for  succeeding  average  M-points.   That  is,  it  apparently  made  no  difference 
whether  the  exact  linear  M-point  equation  was  solved  completely  to  calcu- 
late the  first  peak  greater  than  the  saturation  voltage  or  that  the  first 
peak  greater  than  the  saturation  was  calculated  by  X  =  e  where 

the  linear  system  ^  and  UJ^   are  used.   Thus  an  initial  amplitude  of 
the  dominant  channel  signal  greater  than  saturation  could  be  assumed  as 
long  as  this  amplitude  does  not  cause  saturation  in  the  other  channel 
through  the  relationship 

"Ke 

Another  observation  to  be  noted  was  that  after  both  channels  saturate,, 

the  average  locus  departs  from  a  horizontal  or  vertical  line  as  the  case 

may  be.   If  a  straight  line  is  drawn  through  the  next  couple  of  calculated 

points,  this  straight  line  intersects  the  stability  curve  quite  close  to 

the  experimental  limit  cycle.   This  observatior   is  indicated  in  Figure  IX-1 

This  observation  is  probably  coincidental,  but  some  interpretation  may  be 

made  in  the  future. 
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APPENDIX  A 


General  Equations  of  Mitrovic's  Method  . 

If  the  overall  system  characteristic  equation  is, 

CLnS  -/-£tn-is  ~  +  <2.0_2  ^~Z+ +CL\  S  +  a0  =  O 

then  Mitrovic's  general  mapping  equations  required  in  this  thesis  are: 

1.   For  coefficients  a  and  a,  assumed  variable, 

o      1 


4-, 


and 


2.   For  coefficients  a  and  a  assumed  variable, 


and 


— ^  * 

The  0  functions  in  the  above  equations  are  various  functions  of  ^L 

only  which  occur  in  the  derivation* of  Mitrovic's  mapping  equations.   The  0 
functions  are  as  follows: 

0,1  W  =  -/ 

fa  IM)  * '2.  £ 

I 

I 

Thaler,  G.  J.  and  Ohta,  T.,  Mitrovic's  Method  -  Some  Fundamental 
Techniques,   United  States  Naval  Postgraduate  School,  Research  Paper,  No. 
39,  Jan  1964,  2 


From  the  above  listings  a  general  formula  evolves  from  which  next  0 
function  can  be  formed  by  knowing  the  two  preceding  0  functions.   This 

formula  is: 

For  convenience,  the  values  of  the  0  functions  for  U^     -   0  are  given 


b  a  1  cw  „ 


fcP)  =  *  1 


Thaler,  G„  J.  and  Brown,  R„  G.,  Analysis  and  Design  of  Feedback 
Control  Systems,  McGraw-Hill,  I960,  Ch  10. 
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APPENDIX  B 
Tabulated  Results  of  Analog  Computer  Study  with  E  .  =  Eg2  =  3  volts. 

sle  I  gives  the  results  of  the  first  set  of  computer  runs  travers- 


ing the  B  versus  B1  plane  with  a  linear  system  M-polnt. 


TABLE  I 


Linear 

System 

Resu 

Iting  Limit   Cycle 

_ , — , 

M- Point 

K 

K 

60 

X, 

N,     , 

x0 

N0    . 

e 

L 

Imax 

1mm 

2max 

2mm 

0.1 

Stable 

for  E.    up 

to 

200  vol 

ts 

0.5 

Stable 

for  E .    up 

to 

200   vol 

ts 

4.5 

3.65 

0.101 

1.0 

3.31 

0.907 

0,5 

5.0 

3.64 

0.106 

1.0 

3.97 

0.756 

5,5 

3.63 

0.111 

1.0 

4.6 

0.653 

6.0 

3.63 

0.118 

1.0 

5.0 

0.60 

8.0 

3,63 

0.126 

1.0 

7.2 

0.417 

10.0 

3.62 

0.132 

1.0 

9.22 

0.326 

0.1 

1.19 

3.62 

0.829 

0.435 

1.0 

0.2 

Stable 

for  E.    up 

to 

150  vol 

ts 

4.0 

3.51 

0.207 

1.0 

3.62 

0.807 

1.0 

4.5 

3.49 

0.242 

1.0 

3.8 

0.789 

5.0 

3.49 

0.266 

1.0 

4.45 

0.674 

6,0 

3.49 

0.283 

1.0 

5.65 

0.531 

8.0 

3.48 

0.302 

1.0 

8.0 

0.375 

10.0 

3.47 

0.314 

1.0 

10.3 

0.291 

0.2 

1.3 

5.08 

0.591 

0.895 

1.0 

1.5 

0.4 

1.51 

3.58 

0.838 

1.48 

1.0 

0.5 

1.6 

3.02 

0.994 

1.57 

1.0 

0.55 

Stable 

for  E .    up 
l 

to 

100  vol 

ts 

1.0 

Stable 

for  E,    up 

to 

100  vol 

ts 

3.8 

3.37 

0.384 

1.0 

3.35 

0.92 

4.5 

3.36 

0.45 

1.0 

4.5 

0.667 

1.5 

5.0 

3.36 

0.465 

1.0 

5.2 

0.577 

6.0 

3.36 

0.49 

1.0 

6.46 

0.464 

8.0 

3.35 

0.515 

1.0 

9.25 

0.324 

10.0 

3,34 

0.525 

1.0 

11.6 

0.269 
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Table    T    (Cont'd) 

Linear   System 

Resulting   Limit   Cycle 

M- Point 

K 

e 

Kt 

Ud 

X. 

Imax 

imm 

X0 
zmax 

2mrn 

0,2 

1,3 

7.28 

0.413 

0.95 

1.0 

0.5 

1,595 

4.9 

0.613 

1.94 

1.0 

0,7 

1,76 

3.89 

0.772 

2.46 

1.0 

0,9 

1,92 

3.19 

0.94 

2.79 

1.0 

0.97 

1.97 

3.04 

0.987 

2.82 

1.0 

2,0 

1,1 

Stable 

for  E .    up 

to  75  volt 

s 

2,0 

Stable 

for  E,    up 

to   75  volt 

5 

3,3 

3.15 

0.&55 

1.0 

3.32 

0.904 

3,5 

3.15 

0.705 

1.0 

3.72 

0.806 

4„0 

3,14 

0.755 

1.0 

4.63 

0,648 

6,0 

3,12 

0.86 

1.0 

7.72 

0.388 

8,0 

3.11 

0.9 

1.0 

10.7 

0.28 

10,0 

3,11 

0,925 

1.0 

13.7 

0.219 

0.7 

1,755 

5.25 

0.572 

2.72 

1.0 

0,8 

1,84 

4,88 

0.614 

3.0 

1.0 

0,9 

1,9 

4.55 

0.659 

3.28 

0.915 

1.5 

2.09 

3.82 

0.785 

5.0 

0.6 

1.7 

Stable 

for  E ,    up 

to  62.5  vo 

Its 

2,5 

2.84 

1,63 

1.0 

3.07 

0.977 

2,4 

3.0 

2.82 

1.28 

1.0 

4.5 

0.667 

4.0 

2.78 

1.53 

1.0 

6,95 

0.432 

5.0 

2.75 

1.6 

1.0 

9.5 

0.316 

5.3 

2,20 

3.75 

0.759 

17.8 

0.1686 

5.5 

2.18 

3.79 

0.772 

18.5 

0,1622 

6.0 

2,19 

3.7 

0.811 

19.8 

*    1515 

8.0 

2,17 

3.59 

O.B37 

27.2 

C-iloZ 

10.0 

2.18 

3.72 

0.317 

32.4 

O.092& 

0,2 

1.29 

9,2 

0.326 

0.96 

1.0 

0.5 

1.595 

6,25 

0.48 

1.98 

1.0 

0.7 

1.76 

5.25 

0.572 

2.6 

1.0 

0.8 

1.83 

4.85 

0.619 

2.88 

1.0 

0.9 

1.89 

4.6 

0.652 

3.14 

0.955 

1,0 

1,98 

4.4 

0,682 

3.4 

0.883 

1.5 

2,05 

3.96 

0.758 

4.85 

0,618 

2,5 

2.0 

2.1 

3.82 

0.785 

6.45 

0.465 

2.5 

2,13 

3.76 

0.798 

8.0 

0.375 

3.0 

2,14 

3.72 

0.807 

9.5 

0,316 

3.5 

2.145 

3.7 

0.811 

11.1 

0,27 

4,0 

2,16 

3.7 

0.811 

12.5 

0.24 

5.0 

2,17 

3.68 

0.815 

15.9 

0.1887 

6.0 

2.2 

3.7 

0.811 

19.3 

0.1555 

8.0 

2.15 

3.68 

0.815 

25.2 

0.119 

10,0 

2,14 

3.68 

0.815 

32.0 

0.0937 
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Table   I    ( 

[Cont'd) 

Linear 

System 

Re 

^suiting   Limit   Cycle 

M-Poi 

nt 

K 

e 

t 

(JO 

Imax 

1mm 

X 
2max 

2mm 

0.2 

1.295 

11.2 

0.268 

0.965 

1.0 

0.5 

1.595 

7.65 

0.392 

2.04 

1.0 

0.7 

1.755 

6.4 

0.469 

2.78 

1.0 

0.8 

1.83 

6.3 

0.477 

3.08 

0.974 

1.0 

1.94 

5.7 

0.527 

3.68 

0.815 

1.5 

2.0 

5.5 

0.545 

5.55 

0.54 

3.0 

2.0 

2.04 

5.45 

0.55 

7.25 

0.413 

2.5 

2.025 

5.25 

0.597 

9.0 

0.333 

3.0 

2.06 

5.27 

0.569 

10.8 

0.278 

3.5 

2.075 

5.25 

0.582 

12.5 

0.24 

4.0 

2.05 

5.35 

0.56 

14.2 

0.211 

5.0 

2.05 

5.35 

0.56 

18.3 

0.164 

2.065 

5.3 

0.567 

21.6 

0.139 

8.0 

2.07 

5.25 

0.582 

27.4 

0.1095 

10.0 

2.085 

j  «  ££ 

0.574 

35.4 

0.0847 

1.0 

1.935 

6.65 

0.451 

3.55 

0.845 

2.0 

2.02 

6.55 

0.458 

6.98 

0.43 

3.5 

4.0 

2,05 

6.15 

0.488 

13.9 

0.216 

6.0 

2.05 

6.2 

0.484 

?°.8 

0.1445 

8.0 

2.045 

6.3 

0.476 

28.2 

0.1065 

1.0 

1.92 

7.65 

0.392 

3.59 

0.836 

2.0 

2.03 

7.15 

0.387 

6.75 

0.444 

4.0 

4.0 

2.045 

7.15 

0.42 

13.9 

0.216 

8.0 

2.03 

7.2 

0.417 

28.4 

0.1055 

1.0 

1.925 

8.7 

0.345 

3.59 

0.836 

2.0 

2.01 

8.25 

0.364 

7.1 

0.423 

4.5 

4.0 

2.04 

8.15 

0.368 

14.1 

0.213 

8.0 

2.035 

8.35 

0.359 

21.4 

0.1405 

8.0 

2.035 

8.3 

0.361 

28.6 

0.105 

1.0 

1.915 

9.7 

0.309 

3.62 

0,829 

2.0 

2.00 

9.25 

0.324 

7. 15 

0.42 

5.0 

4.0 

2.035 

9.1 

0.33 

14.1 

0.213 

6.0 

2.025 

9.2 

0.326 

21.6 

0.139 

8.0 

2.025 

9.15 

0.328 

28.8 

0.1045 

.,  .  0 

1.915 

10.7 

0.28 

3.62 

0.829 

2.0 

2.00 

10.07 

0.298 

7.15 

0.42 

5.5 

4.0 

2.03 

10.05 

0.299 

14.2 

0.211 

6.0 

2.03 

10.0 

0.30 

21.5 

0.1395 

8.0 

2.025 

10.02 

0.294 

28.9 



0.104 
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Table   I    (Cont'd) 

Linear   System 

-■ 

Resulting   Limit   Cycle 

M-Point 

K 

e 

Kt 

UD 

Imax 

Imin 

2max 

2mm 

1.0 

1.915 

11.7 

0.256 

3.63 

0.827 

2.0 

2,00 

11.07 

0.271 

7.15 

0.42 

6.0 

4.0 

2,025 

11.05 

0.272 

14.2 

0.211 

6.0 

2.035 

11.0 

0.273 

21.4 

0.1405 

8.0 

2.025 

11.2 

0.268 

29 

0.1035 
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Table  II  gives  the  resulting  limit  cycle  for  linear  M-point  close 


to  dividing  line  No.  1, 


TABLE  II 


— — n 

Linear  System 

Resulting  Limit  Cycle 

M- point 

K 

e 

\ 

LO  ave 

lmax 

2max 

0.90 

1.92 

3.52 

2.98 

0.92 

1.93 

3.50 

3.0 

2.1 

0.96 

1.955 

3.42 

3.08 

1,0 

1.99 

3.34 

3.19 

0.86 

1.895 

4.2 

2.97 

0.88 

1.91 

4.1 

2.99 

2.3 

0.90 

1.915 

4.05 

3.03 

0.92 

1.925 

4.0 

3.07 

0.80 

1.85 

4.9 

2.89 

0.82 

1.87 

4.83 

2.95 

2.5 

0.84 

1.88 

4.75 

2.99 

0.86 

1.89 

4.7 

3.03 

0.88 

1.91 

4.65 

3.08 

0.90 

1.915 

4.6 

3.12 

0.80 

1.83 

6.0 

2.96 

3.0 

0.82 

1.85 

5.85 

2.98 

0.84 

1.865 

5.8 

3.02 

0.80 

1.83 

8.05 

2.96 

4.0 

0.82 

1.85 

8.0 

3.02 

0.79 

1  82 

10.35 

2.32 

5.0 

0.80 

1.83 

10.25 

3.0 

0.81 

1.84 

10.2 

3.04 

0.78 

1.815 

12.1 

2.98 

6o0 

0.79 

1.83 

12.4 

3.02 

0.80 

1.84 

12.5 

3.03 

. 
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Le  III  gives  the  resulting  limit  cycles  for  linear  M-points  close 

to  dividing  line  No.  2. 

TABLE  III 


Linear  System 
M-  point 

Resulting  Limit  Cycle 

Kt 

K 

e 

UJ^V^ 

X. 
Imax 

x0 

zmax 

3.0 

2.40  . 
2.42 
2.43 
2.45 

2.78 
2.73 
2.08 
2.08 

1.31 
1.49 
3.83 
3.88 

4.5 

4.68 

9.65 

9.75' 

4.0 

2.35 
2.38 
2.39 
2.41 

2.42 

2.83 

2.74 
2.71 
2.68 
2.09 

1.25 
1.42 
1.57 
1.58 
3.68 

6.1 
6.51 
6.7 
6.85 

12.8 

5.0 

2.40 
2.41 
2.42 

2.68 
2.11 
2.12 

1.51 

3.7 

3.68 

8.6 
16.0 
16.5 

6.0 

2.39 

2.40 

2.73 
2.11 

1.52 
3.66 

10.3 
19.0 

8.0 

2.38 

2.39 
2.40 

2.75 
?.13 
2.15 

1.60 
3.62 
3.65 

14.3 
27  4 
27  8 

■ 

10 

2.38 
2.39 

2.74 

M6 



1.62 
3.7 

18.7 
32.8 
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